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AB8TRACT 

Taylor  series  have  a  long  history  of  usefulness  in  numerical  analysis* 
especially  for  the  nuserical  solution  of  the  initial  value  problem  for  systems 
of  ordinary  differential  equations.  Since  recurrence  relations  for 
coefficients  of  Taylor  series  are  well  known*  it  is  possible  to  automate  the 
computation  of  arithmetic  operations  and  various  standard  functions  with 
arguments  which  are  themselves  series.  If  the  language  used  for  scientific 
computation  supports  user  defined  operators  and  data  types,  then  the 
facilities  built  into  the  language  compiler  itself  can  be  used  to  generate 
machine  code  for  the  evaluation  of  Taylor  coefficients.)  Examples  of  such  t 
languages  are  Pascal -SC*  Algol  68*  and  ADA  (a  trademark  of  the  U.  S.  ,_r 
Department  of  DefenseTT^tPascal-SC  (Pascal  for  Scientific  Computation)^ offers 
the  user  highly  accurate  floatingpoint  and  interval  arithmetic*  the  latter 
being  useful  for  automatic  computation  of  guaranteed  error  bounds  .*-~~Ia^this 
language,  series  with  real  coefficients  are  introduced  as  type  TAYLOR,  and^the 
corresponding  series  with  interval  coefficients  as  type  ITAYLOR.  Source  code 
is  given  for  the  operators  +/  *,  /,  **  and  the  functions  SQR,  SQRT*  EXP, 

SIN*  COS,  ARCTAN*  and  LN  with  arguments  of  these  types  and  Borne  other  useful 
functions  and  procedures.  Integer*  real,  and  interval  constants  are  also 
allowed  in  TAYLOR  or  ITAYLOR  expressions.  Suggestions  for  the  implementation 
of  additional  operators  or  functions  are  given.  An  application  of  Taylor 
series  and  the  methods  of  interval  analysis  t  ,  the  solution  of  the  initial 
value  problem  for  ordinary  differential  equations  is  made  using  types  TAYLOR 
and  ITAYLOR.  An  analysis  of  the  stability  of  this  method  is  made,  which  shows 
that  the  recurrence  relations  for  generation  of  the  Taylor  series  for  the 
solution  exhibit  a  mild  instability  which  has  no  significant  effect  on  the 
values  of  the  solution  computed  by  analytic  continuation. 
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SIGNIFICANCE  AND  EXPLANATION 


The  reliable  numerical  solution  of  the  initial  value  problem  for  systems  of 
ordinary  differential  equations  is  one  of  the  fundamental  problems  of  scientific 
computation.  Taylor  series  methods  for  this  purpose  have  long  been  recognized 
to  be  of  theoretical  importance,  but  their  use  in  practice  has  been  hampered  in 
the  past  by  the  need  to  differentiate  the  functions  defining  the  system  in  order 
to  obtain  coefficients  of  the  series.  Consequently,  methods  are  commonly  used 
which  take  linear  combinations  of  function  evaluations  to  be  "as  good  as"  Taylor 
polynomials  of  some  degree  as  approximate  solutions.  However,  since  recurrence 
relations  are  well  known  for  the  Taylor  series  coefficients  of  functions 
resulting  from  arithmetic  operations  and  standard  functions  involving  series 
arguments,  it  is  possible  to  generate  the  required  coefficients  automatically 
with  a  computer,  and  the  user  need  only  supply  code  for  the  functions  defining 
the  differential  equations.  Many  modem  compilers  allow  user  defined  data  types 
and  "overloading"  the  standard  operators  and  functions  so  that  these  series 
operations  can  be  applied  directly  to  vectors  consisting  of  the  coefficients  of 
the  Taylor  polynomials  involved.  Examples  of  such  languages  are  Pascal-SC, 

Algol  68,  and  ADA*  (a  trademark  of  the  Department  of  Defense).  Pascal-SC  offers 
many  advantages  for  numerical  computation,  since  it  is  based  on  accurate  real 
and  interval  arithmetic  for  vectors  as  well  as  scalars.  In  this  report,  source 
code  is  given  for  generation  of  both  real  and  interval  valued  Taylor  series,  and 
simple  programs  illustrate  the  application  of  each  to  the  numerical  solution  of 
ordinary  differential  equations.  In  these  programs,  both  the  order  of  the 
method  being  used  and  the  step  size  of  the  integration  are  under  the  control  of 
the  user.  Advantages  of  interval  computation  include  the  ability  to  study  the 
range  of  values  of  solutions  depending  on  ranges  of  initial  conditions  and 
parameters  in  the  equations,  and  to  obtain  rigorous  bounds  for  solutions  in 
applications  such  as  aerospace  in  which  reliability  is  important.  In  this 
report,  interval  calculations  are  used  to  investigate  the  stability  of  the 
process  of  generating  the  Taylor  coefficients.  For  rapidly  converging  Taylor 
series,  instability  in  the  generation  of  the  coefficients  has  little  effect  on 
the  computed  values  of  the  solution  because  as  the  relative  error  increases  for 
successive  terms,  the  terms  themselves  become  smaller.  Hence,  their 
contribution  to  the  error  of  the  sum  is  insignificant. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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AUTOMATIC  GENERATION  OP  TAYLOR  SCRIES  IM  PASCAL- SC : 

BASIC  OPSRATIONS  AMO  APPLICATIONS  TO  ORDINARY  DIFFERENTIAL  EQUATIONS 


George  Oorliss*  and  L.  B.  Rail** 


1.  Taylor  series.  polynomials.  and  forma.  A  fundamental  tool  of  nuawrical 
analyaia  la  tha  axpanalon  of  a  raal  function  f  of  a  raal  var labia  x  Into  a  Taylor  aarlaa  at 
x  -  x0,  which  gives  tha  axpraaalon 


f(x)  -  t  fU*1,(xn)(x  -  xn)<1"1)/(i-1)t, 
1-1  " 


valid  for  |x  -  xQ|  <  p,  whara  p  la  tha  radiua  of  convergence  of  tha  inflnlta  aarlaa  on  tha 
right-hand  alda  of  ( 1 . 1 )  •  Of  course,  In  actual  numerical  computation,  tha  Taylor 
polynomial 


f  (x)  -  Z  f,i"1,(x.)(x  -  xfl)U-1)/(i-1»l< 
B  1-1  00 


la  uaed  In  placa  of  tha  inflnlta  aarlaa,  Ihia  raaulta  In  the  truncation  error 


<1.3)  R  <f,jVxB)  -  f<x)  -  f  <x>  -  f(B,Ct)(x  -  xft)B/ni,  5  e  X, 

nun  v 


whara  X  denotes  tha  Interval  X  >  {•in{x,x0},max{x,x()}] ,  and  tha  remainder  term  Rn(f,x0»x) 
la  axprasaad  In  Lagrange  form.  This  approximation  of  fix)  by  fn(x)  gives  rise  to  a  problem 
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of  error  MtiMtiw  which  can  be  aolvtd  by  thi  methods  of  interval  analysis.  If  la 


interval  Inclusion  of  the  rsal  function  f'  ,  than 


f(x>  -  fn(x)  -  l^tf  ,x0/x)  a  F'n,U)(x  -  x^/nli 


this  allows  automatic  computation  of  guaranteed  arror  bounds  by  tha  usa  of  Interval 
arithmetic  [12],  [13].  - - 


In  order  for  Taylor  series  methods  to  be  useful  in  scientific  oomputatioo,  it 
be  possible  to  automats  the  calculation  of  the  no  mall sad  real  Taylor  coefficients 


c(i  ♦  1)  -  f^NxjXx  -  XglVil,  i  “  1,...,n-1, 


and  the  corresponding  interval  quantities 


C(1  ♦  1)  -  r(l)lX)(x  -  x*)1/!!,  1  -  1,...,n-l. 


These  calculations  can  be  carried  out  by  aeane  of  well-known  recurrence  relations  [1], 
[12],  [13],  (161  for  functions  defined  by  subroutines  or  expresslone  involving  arithmetic 
operations  and  a  variety  of  standard  functions  for  which  library  subroutines  are 
available.  A  very  important  application  of  automated  generation  of  Taylor  series  by 
recursion  Is  the  numerical  solution  of  the  initial-value  problem  for  ordinary  differential 
equations.  That  is,  it  is  required  to  find  y  “  y(x)  »  (y^Cx) ,y2(x) , . . . ,yK(x) ]  euch  that 


y]  “  fjlx.y),  TjU,,)  -  yJ0,  i  -  1,  •••»■» 


for  values  of  x  in  an  Interval  containing  Xq  [1],  [3],  [6],  [12],  [13]. 

Another  application  of  the  methods  in  this  paper  is  to  the  automatic  generation  of 
interval  inclusions  of  real  functions  by  means  of  their  interval  mean- value  and  Taylor 


forms  [12],  [13],  [19].  Suppose,  for  example,  that  f(x)  is  e  real  function,  such  as 


fix)  -  (x  4  3)/(x*  4  2). 


which  can  be  evaluated  by  tha  correapoedi ng  expression 


f  »-  (x  ♦  3)/(x**2  4  2)> 


In  a  Fascal-SC  program.  An  Interval  Inclusion  F  of  f  on  an  latarwaX  X,  for  which 


(1.10) 


f(X)  -  {fix)  |  X  a  x)  c  fix) 


can  be  obtained  alnply  by  declaring  tha  variables  F  and  X  to  ba  of  type  IHTEKVAL,  and  than 
evaluating  the  axpraaalon  corresponding  to  (1.9), 


(1.11) 


F  i-  (X  4  3>/<X**2  4  2) 


using  interval  arithmetic,  a  standard  feature  of  Fascal-SC  [23] .  An  inclusion  obtained  in 
this  way  nay  ba  too  coarse  In  tha  Sanaa  that  FIX)  la  a  such  larger  Interval  than  needed  to 
contain  f(X).  In  this  case,  an  interval  inclusion  provided  by  the  naan-value  fora 


(1.12) 


F, (X)  -  fix)  4  F’ IX) IX  -  x),  X  «  X, 


can  ba  batter,  particularly  if  tha  width  of  X  is  not  large  [121,  1131 ,  [19].  In  11.12),  F* 
denotes  an  Interval  inclusion  of  the  derivative  f'  of  fi  F*(X)  Is  obtained  automatically  by 
evaluating  (1.11)  with  F  and  X  of  type  XTAYLOR,  as  will  ba  explained  below.  Interval 
inclusions  of  f  are  also  provided  by  ISylor  forms  of  higher  order  (19),  in  general. 


(1.13) 


F  IX)  -  I  fU,(x)(X  -  x)1/!!  ♦  F<B>(X)(X  -  x)B/nl,  x  «  X. 
"  1-0 


■3' 


Am  torn  c«n  be  fm*r*t«4  automatically  from  the  nprmloni  (1.9)  and  (1.11)  by  tha  iaa 
of  types  TAYLOR  and  ITAYLOR,  respectively.  Recursive  generation  of  real  and  Interval 
Taylor  coafficionta  makes  poaaibla  an  adaptive  Method  for  calculation  of  interval 
lnclueiona  of  real  functional  in  which  n  ia  increaaed  until  ?„(*)  includea 
r  ,(X).  Zt  ia  alao  poaaibla  to  reduce  tha  width  of  confuted  lnclueiona  by  making  uae  of 
the  fact  that  the  interaectlon  of  Interval  lncluaiona  ia  likewise  an  Interval  lnclualon. 

Previous  implementations  of  automatic  generation  of  Taylor  coefficients  in  computer 
languages  such  aa  FORTRAN  have  used  interpretation  [20]  or  pre-compilation  [10]  to  activate 
the  necessary  subroutines  [Id].  In  more  aodarn  languages,  the  caapilsr  itself  can  be  used 
to  produce  the  necessary  routines,  leading  to  a  saving  of  programming  effort  and  an 
Increase  in  clarity  of  the  source  code.  The  uae  of  Pascal-SC,  a  language  of  this  type, 
will  be  explained  in  the  next  section. 

2.  Pascal-SC.  The  aethod  for  autoaatic  generation  of  Taylor  aeries  given  in  this 
report  ia  based  on  coaputatlon  with  the  coefficients  of  Taylor  polynomials  of  arbitrary 
length,  considered  as  specific  aatheaatical  entitles.  This  requires  that  the  language 
support  1)  user  defined  data  types,  as  do  descendants  of  ALGOL-60  such  as  Pascal  and  ADA 
(ADA  is  a  trademark  of  the  0.8.  Ospartswnt  of  Defense);  and  11)  user  defined  operators,  as 
do  ALGOL-68  and  ADA. 

Pascal-SC  [2]  is  an  extension  of  Pascal  irtilch  provides  both  user-defined  data  types 
and  user-defined  operators.  This  paper  assumes  a  modest  familiarity  with  standard  Pascal 
[8] .  For  the  reminder  of  this  Section,  we  outline  some  of  the  extensions  which  aaks 
Pascal-SC  well  suited  to  the  applications  in  this  paper.  The  reader  who  wishes  to  omit  the 
discussion  of  programming  language  issues  my  proceed  directly  to  the  definition  of  the 
data  types  TAYLOR  and  ITAYLOR  in  Section  3. 

Pascal-SC  was  developed  with  the  needs  of  scientific  computation  in  Bind.  It  is  an 
implementation  of  Jensen  and  Wlrth  Pascal  [S]  which  also  provide*  intervals,  complex 
numbers,  complex  Intervals,  as  built-in  elamantary  scalar  data  types  [23] .  A  full  range  of 
standard  operators  is  provided  to  manipulate  the  elementary  scalar  data  types,  as  well  ee 
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vectors  and  Matrices  built  of  thaaa  typos  (23) . 

Standard  Pascal  supports  uaar— defined  data  typos  built  from  elementary  data  typos* 
this  feature  will  bo  used  to  dsfiao  variables  of  typo  TAYLOR  and  XTAYLOR  ( interval  Taylor) 
In  Section  3. 

Pascel-SC  allows  the  user  to  doflno  operators.  Most  computer  lan qua 90s  allow 
programmers  to  define  functions,  subroutines,  or  procedures,  but  except  for  APL,  the 
lan9US9es  nost  often  used  for  scientific  computation  require  that  such  user- defined 
functions  be  called  using  a  prefix  notation  (eg.  8XM  (X)),  while  built-in  operators  are 
called  using  an  infix  notation  (eg.  A  ♦  B).  Programmers  can  define  operators  in  Pascal-SC 
to  extend  the  language  in  a  uniform  way,  retaining  the  familiar  infix  notation  for 
operators  whose  operands  are  variables  of  ueer-deflnsd  types  (eg.  A  +  B,  where  A  and  B  are 
variables  of  type  TAYLOR). 

Operators,  functions,  and  procedures  in  Paaeal-8C  can  be  overloaded,  that  is,  the 
name  of  an  operator,  a  function,  or  a  procedure  can  have  different  meanings,  depending  on 
the  type  or  number  of  its  operands.  For  example,  the  standard  Pascal  or  FORTRAN  operator 
is  said  to  be  overloaded  because  aA  ♦  B"  for  Integer  variables  A  and  B  has  a  different 
meaning  from  "A  ♦  B"  for  real  variables  A  and  B.  The  support  of  Pascal-SC  for  overloading 
of  user  defined  operators  is  essential  to  the  uniform  extension  of  the  language  because  we 
wish  to  define  the  Marling  of  "a  ♦  B“  for  variables  which  represent  Taylor  series  with  real 
or  with  interval  coefficients. 

The  support  of  Pascal-SC  for  user-defined  operators  and  for  overloading  is  very 
similar  to  that  provided  by  ADA.  ADA's  PACKA6B  concept  would  allow  a  more  secure 
implementation  of  data  abstractions  (9)  for  real  and  interval  valued  Taylor  series.  Pie 
operations  on  intervals,  however,  also  require  support  for  directed  rounding  of  floating¬ 
point  respite  in  order  to  guarantee  that  the  desired  answer  is  contained  in  the  interval 
computed.  The  early  implementations  of  ADA  do  not  provide  an  accuracy  of  floating-point 
computations  which  can  ccmpets  with  Pascal-SC. 

Pascal-SC  features  a  highly  accurate  arithmetic  based  on  a  general  theory  (111  for 
real  and  complex  maabers,  real  and  complex  intervals,  and  vectors  and  Mtrlces  over  these 
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typM.  Operations  on  floating-point  numbers  are  rounded  to  tha  closest  floating-point 
nabtr  to  the  true  result,  or  upward  or  downward  to  the  closest  neighboring  floating-point 
number  under  the  oontrol  of  the  user.  this  accuracy  meets  the  proposed  ins  standard  for 
floating-point  arithmetic  [14].  In  addition,  scalar  products  of  sectors 

N 

(2.1)  SCALP! A,B,ROOMD)  -  I  A  *B. 

i-1  1  l 

are  calculated  with  the  sees  accuracy  (to  the  closest  floating-point  number),  and  with  the 

tint  options  for  rounding  [23] .  A  sufficiently  long  accumulator^ is  ueed-to  store —  - 

intermediate  results  in  the  evaluation  of  the  scalar  (or  inner)  product,  this  capability 
can  also  be  used  to  obtain  results  of  the  same  high  accuracy  in  evaluation  of  a  given 
arithmetic  expression  so  that  1.0B+99  ♦  1. OB-99  -  1.0B+99  yields  1. OB-99. 

3.  Types  TATLOR  and  ITATLOR.  Me  wish  to  provide  the  developer  of  scientific 
software  with  a  set  of  tools  with  which  Taylor  series  methods  can  be  implemented  easily  for 
a  variety  of  noser leal  problems.  The  ability  of  the  computer  to  perform  formula 
translation  is  used.  Compilers  sinew  the  first  rORTRAJf  compiler  have  produced  machine  code 
for  the  evaluation  of  an  expression  such  as 

(3.1)  T  f  (**T  ♦  SI MIX)  *  4.0)  •  (3.0  •  (T«*2)  +  6.0). 

This  is  done  by  analysis  of  the  expression  and  application  of  the  rules  for  evaluation  of 
formulas.  If  the  rules  for  differentiation  or  recursive  generation  of  Taylor  coefficients 
are  applied  in  the  same  way,  then  code  for  the  evaluation  of  the  corresponding  quantities 
results  [16] .  In  this  way,  fast  and  inexpensive  operations  performed  by  the  compiler  avoid 
the  overhead  Involved  in  invoking  symbolic  differentiation  software.  This  leads  to  a  more 
efficient  implementation  of  Taylor  series  generation  all  the  way  from  initial  coding 
through  program  execution. 

The  normalised  Taylor  coefficients  of  a  function  f(x)  expanded  at  x  ■  Xg  are 


dmtinmd  bf 


e.rci*  *  «  -  t  •  *  -  x  -  o,i, a,. 
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whara  DIN  la  tka  lan  9th  of  tha  truncated  sarlaa  which  is  actually  a  tor  ad.  this  raal  oar 
lntarval  vactor  of  noraallsad  Taylor  ooaffloiants  is  tha  basis  tor  tha  data  typos  nnOR 
and  iraTLOR.  for  tha  raaalndar  of  this  papar,  tha  tarn  *  serin*"  la  naad  to  rafar  to  tha 
Taylor  polynomial  given  by  aquation  (3.3)  or  its  lntarval  anal 09. 

in  what  follows,  tha  general  rula  will  ba  adoptad  that  all  varlahlas  or  expressions 
of  tha  scalar  typaa  MTMSR.  MU,  or  INTERVAL  ara  traatad  as  constants  for  tha  purposas  of 
diffarantiation. 

To  font  tha  raal  data  typa  TAYLOR,  tha  DZN  nornaliaad  Taylor  ooaffielanta  in  (3.3) 
ara  storad  as  a  vactor  of  floating-point  n unbars,  tha  appropriate  declarations  in  Pascal- 
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nnn  hWn  an  *-• -  - 4m  flAI  anwnh  Iav  ika  nM  n  snail 
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F. LENGTH  <  DIM  if  F  is  being  built  up  recursively!  if  F  has  been  defined  by  term-by-term 
differentiation  of  another  series,  or  if  F  has  been  defined  as  a  quotient  of  two  series 
both  of  whose  leading  tens  are  zero  (see  Section  4.2).  This  field  has  been  added  to  the 
record  for  type  TAYLOR  given  in  [18]  for  internal  documentation  and  so  that  i)  only  series 
terms  actually  uaed  need  to  be  processadi  and  ii)  1' Hospital's  rule  can  be  applied  to 
certain  indeterminant  forms  0/0  which  may  appear. 

The  normalised  Taylor  coefficients  themselves  are  stored  in  the  array  of  real 
nuabers  named  TC,  that  la. 


F.TCtK]  -  F<K-1,(X0)(X  -  X0)(K"1)/KI,  X  -  1,...,DXM. 


The  size  of  the  stap  being  used  for  expansion  is  F.T  -  X  -  Xfl.  Series  are 
generated  uaing  a  fixed  stepslze  for  which  the  series  Slight  even  be  divergent.  The  series 
for  F  at  a  different  point  Z  la  readily  computed  at  a  cost  proportional  to  DIMs 


F.TCtK]  i-  F.TC[K]*((Z  -  X0)/F.T))s*(K-1)»  K-2,...,DXM, 


while  the  cost  of  series  generation  is  usually  proportional  to  DIM2.  The  presence  of  the 
stepslze  in  the  record  also  mates  it  possible  to  check  that  an  operation  is  not  being 
performed  on  two  series  with  different  stepsises. 

One  of  the  Important  problems  to  which  interval  analysis  has  been  applied  since  its 
beginnings  is  the  problem  of  controlling  the  truncation  error  of  Taylor  series  methods 
[12].  Hence  it  is  natural  to  support  Taylor  series  whose  normalized  coefficients  are 
intervals.  The  appropriate  declarations  in  Pascal-9C  are 


CONST  DIN  -  nr 


{  Oser  supplies  n  } 


KR  DIMTYPE  -  1..DXNI 

INTERVAL  -  RECORD  INF,  SUP  I  RIAL  END, 

(3.7)  I VICTOR  »  ARRAY [ DIMTYPE}  Of  INTERVAL! 

XTAYLOR  “  RECORD  LENGTH  I  DINTYPEi 

T  t  REAL! 

TC  t  I VECTOR  END* 

Aa  types  ITAYLOR  ahd  TAYLOR  are  As  mum,  except  that  Ae  nonuliud  coefficients  of  tha 
former  ara  intervals.  Tha  aaaa  racarranca  ralationa  are  used  to  generate  aarlaa  of  aach 
typa. 

Tha  atapalza  T  remains  raal.  Ala  oorraaponda  to  bounding  tha  range  of  values  of  a 
function  f  at  ona  raal  n unbar  x.  Thara  ara  scam  applications  for  which  it  la  nacaaaary  to 
bound  tha  rakja  of  f  on  an  lntarval ,  aa  in  (1.13).  In  Ala  caaa,  ona  can  taka  T  «  1  and 
for*  tha  norma  11 sad  coafflcianta  by  computing  tha  naadad  powers  of  (X  -  x)  by  lntarval 
arithmetic,  or  else  Introduce  a  new  data  typa  in  Alch  T  la  of  typo  INTERVAL,  and  a  sat  of 
operators  corresponding  to  those  given  hare. 

Aa  declarations  (3.6)  and  (3.7)  of  types  TAYLOR  and  ITAYLOR,  respectively,  are 
basic  to  the  discussion  of  operators  in  tha  next  section. 

4.  Implementation  of  operators  and  functions  for  typos  TAYLOR  and  ITAYLOR.  As 
indicated  above,  the  ability  of  a  compiler  to  perform  formula  translation  can  also  be  used 
to  produce  machine  coda  for  Aa  evaluation  of  the  normalised  Taylor  coefficients  [1],  [3], 
(10),  (12),  [13],  (15),  (16),  (18).  If  tha  value  of  function  f  la  obtained  by  a 
composition 

(4.1)  f  ■  f j  o  fj  o  ...  o  fB 

of  a  finite  number  of  elementary  functions,  then  derivatives  of  f  can  be  computed  by  tha 
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chain  rul«  fro*  tha  derivatives  of  f  ^ , . . .  f a .  This  la  a  tedious  and  error-prone  calculation 
to  do  by  hand,  but  tha  computer  doaa  it  not  only  rapidly,  but  alao  accurately. 

Racurranca  ralationa  for  calculating  tha  normalised  Taylor  ooeffidente  for  tha  baalc 
arithmetic  oparatlona  and  for  tha  elementary  functions  are  wall  known  (aaa  [16],  for 
example) •  Hence  machine  coda  can  be  generated  to  expand  the  Taylor  aarlaa  for  f  at  any 
point  x  ■  x0  at  which  f  la  analytic,  These  recurrence  ralationa  are  both  more  efficient 
and  aore  accurate  than  maaerical  differentiation  [1?].  Recursive  generation  of  the  aarlaa 
nay  be  aildly  unstable  {6),  but  the  interval-valued  Taylor  aerlea  introduced  in  Section  3 
can  give  guaranteed  bounds  for  the  effect  of  any  such  Instability.  In  Section  6,  we  show 
that  any  Instability  in  the  aeriea  generation  has  no  significant  effect  on  the  aeries 
sum. 

Rail  [18]  outlines  an  approach  to  abstract  data  types  for  real  and  interval-valued 
Taylor  series.  Our  implementation  generally  follows  that  outline.  This  report  discusses 
extensions  and  soaa  of  the  iapleaentation  details.  The  code  is  included  as  Appendix  C  of 
this  report,  first  of  all,  in  order  for  expressions  to  be  evaluated  correctly  whan  they 
include  variables  of  type  TAYLOR  or  ITAYLOR,  the  arithmetic  operations  and  the  standard 
functions  must  be  defined  in  a  manner  which  Incorporates  the  appropriate  recurrence 
relatione  for  the  generation  of  the  normalized  Taylor  coefficients.  Our  implementation  in 
Pascal -SC  attempts  to  follow  the  principles  of  uniformity,  compactness ,  locality,  and 
linearity  for  a  good  programming  language  design  [22]  .  Next  we  attempt  to  justify 
significant  departures  from  two  of  these  principles. 

The  principle  of  uniformity  in  programming  language  design  says  that  the  mm  things 
should  be  done  the  same  way  whenever  they  occur.  Thus  "A  *  B*  means  "add",  regardless  of 
the  types  of  the  variables  A  and  B.  The  other  arithmetic  operators  enjoy  the  same 
uniformity,  but  the  standard  functions  do  not.  For  example,  exp(x)  is  EXP(X)  if  X  is  RIAL, 
ItXP(X)  if  X  is  INTERVAL,  TEXP(X)  if  X  is  TAYLOR,  and  ITXXP(X)  if  X  is  ITAYLOR.  EXP  and 
I EXP  are  built-in  functions  which  wars  designed  to  suggest  the  type  of  their  operand  and 
result  as  an  aid  to  reading  the  code.  That  is  especially  useful  since  Pascal  tends  to 
violate  the  principle  of  locality  by  placing  the  declaration  of  a  variable  far  from  its 


us«<  Ns  chose  to  maintain  uniformity  of  our  extensions  with  the  built-in  functions'  It  is 
Important  to  be  able  to  determine  the  type  of  a  variable,  sad  It  mould  be  quit*  non-uniform 
if  IEXP  were  the  only  function  in  this  family  which  requires  a  prefix' 

The  principle  of  locality  suggests  that  all  relaxant  parts  of  the  program  are  found  in 
the  same  place.  He  attempt  to  follow  this  principle  in  each  of  our  program  unit*,  but  the 
use  of  the  global  constant  DIM  and  the  global  types  RVRCTOR,  I  VECTOR,  TAYLOR,  and  ITAYLOR 
is  a  xiolation.  The  us*  of  such  global  types  needed  in  the  headings  of  the  operators  and 
functions  is  very  difficult  to  avoid.  Their  us*  has  the  advantage  that  all  of  the 
information  about  the  length  of  the  series  to  be  used  is  located  in  only  one  place, 
const  DIN  -  n,  so  it  la  easy  to  change. 

In  roughly  their  order  of  importance,  the  goals  of  this  lnpleawntatlon  ares 

o  Consistent  set  of  software  tools, 

o  Correct  answer  whenever  possible. 

o  Useful  error  messages  when  no  correct  answer  is  possible, 
o  Readable  cod*  for  future  edaptatlone. 
o  Ifflcient  execution, 

o  Ooagsact  code. 

For  example,  this  implies  that  although  efficient,  compact  code  is  sought,  efficiency  and 
compactness  are  aanetlaws  sacrificed  for  higher  goals.  In  particular,  it  is  important  that 
other  programmers  be  able  to  reed  the  code,  perhaps  in  order  to  improve  its  efficiency. 

Binary  oparations  with  on*  operand  of  type  TAYLOR  may  appear  with  the  other  operand  of  type 
INTBGIR,  REAL,  or  TAYLOR)  and  tha  two  oparands  nay  appear  in  althsr  ordar.  Similarly, 
binary  oparations  with  on*  operand  of  type  I TAYLOR  may  have  a  second  operand  of  type 
ItmtGIR,  interval,  or  ITAYLOR.  The  operators  built  into  Faacal-SC  do  not  support  the 
nixing  of  RIAL  and  INTERVAL  operands  because  real  nutbers  are  viewed  es  being  potentially 
inexact  (23] .  Our  extensions  of  the  arithmetic  operator*  to  interval  valued  Taylor  series 


Maintain  uniformity  with  this  convention.  Uhls  Is  recognised,  but  not  explicitly  stated  in 
tIB] .  U  t  proij i  isnr  Is  certain  that  a  real  lumber  X  is  sxact  so  that  it  eay  safely  be 
mixed  with  an  interval,  IMTST  (X)  converts  X  into  the  interval  [X,  X] • 

The  library  of  subroutines  to  support  computations  with  types  TAYLOR  and  ITAYLOR 
includes  operators  (♦,-,*,/,•*),  special  power  functions  (aqc,  eqrt,  exp),  standard 
functions  (sin,  cos.  In,  arctan),  and  additional  functions  (tan  and  the  Rungs  function 
f(x)  “  1/(1  +  x2) ,  to  which  the  user  can  add  more  functions  and  procedures  as  desired.  The 
analytic  operations  of  term-by-term  differentiation  of  real  and  interval  series,  as  well  as 
term-by-term  differentiation  of  interval  series  are  also  provided  by  means  of  functions  for 
the  given  purpose.  There  is  also  a  set  of  utility  functions  and  procedures  to  perform 
frequently  needed  tasks,  such  as  reading  and  writing  real  and  interval  series,  taking  the 
midpoints  of  the  coefficients  of  an  interval  series  to  obtain  a  real  series,  and  so  on. 

The  following  abbreviations  are  used  in  the  cods  to  make  it  as  easy  as  possible  to 
locate  a  desired  operation  with  any  text  editors 

K  INTEGER 

R  RIAL 

I  INTERVAL 

T  TAYLOR 

IT  ITAYLOR. 

Using  these  abbreviations  to  distinguish  between  Instances  of  overloading,  the  operators 
which  are  needed  to  support  variables  of  type  TAYLOR  and  ITAYLOR  ares 

Addition  (Section  4.1) s 

♦  T,  K  ♦  T,  T  +  X,  R  ♦  T,  T  ♦  R,  T  ♦  T 

♦  IT,  K  ♦  IT,  IT  ♦  X,  I  ♦  IT,  IT  ♦  I,  IT  ♦  IT 

Subtraction  (Section  4.1) t 

-  T,  K  -  T,  T  -  X,  R  -  T,  T  -  R,  T  -  T 
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-  it,  r  -  it, 

IT  -  X, 

R  - 

IT, 

IT  - 

R, 

IT 

-  IT 

Multiplication  (Section  4.2) i 

R  •  T, 

T  *  X, 

R  * 

T, 

T  * 

R, 

T 

*  T 

K  *  IT, 

IT  *  X, 

I  • 

IT, 

IT  • 

I. 

IT 

*  IT 

Oivision  (Section  4.2) i 

K  /  T, 

T  /  X, 

R  / 

T, 

T  / 

R. 

T 

/  T 

K  /  IT, 

IT  /  X, 

I  / 

IT, 

IT  / 

I. 

IT 

/  IT 

Power  (Section  4.3) t 

K  ••  K,  R  ** 

X,  X  *• 

R. 

R 

**  R 

I  •• 

X,  X  ** 

I. 

I 

*«  I 

I  »•  T,  T  *• 

X,  R  ** 

T, 

T 

**  R. 

T 

•* 

T 

K  **  IT,  IT  ** 

X.  I  ** 

IT, 

IT 

**  I# 

IT 

•* 

IT 

Iaplamutlon  total  la  of  «ach  operator  art  discussed  in  the  Sections  shown.  Pascal-SC 
provides  no  power  operator,  so  “  east  be  defined  for  the  scaler  types  before  It  can  be 
extended  to  types  TAYLOR  and  ITinOR*  lbs  discussion  of  •*  Is  postponed  to  follow  the 
Introduction  in  Section  4.3.1  of  special  cases  of  exponentiation t  aqr,  sqrt,  and  exp. 
The  priorities  of  the  operators  9iven  in  this  Section  eret 
Highest)  Unary  addition  and  subtraction,  functions* 

Multiplication,  division,  and  powers!  *,  /,  *• 

bowest<  Binary  addition  and  subtraction!  ♦,  - 

In  particular,  note  that  the  priority  of  ••  relative  to  *  and  /  is  different  than  in 
PORTRAM.  ' 

Tor  types  TAYLOR  and  I T ATI/O It ,  laplesontatlon  has  been  provided  for  the  standard 
functions  which  are  supported  in  Pascal-SC  for  types  INTMER,  RIAL,  and  INTERVAL.  They 
are  i 
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Special  powers  (Section  4.3.1)  i 

TSCX  (  T),  T8QRT  (  T),  TUP  (  T) 

XTSQR  (XT),  ITSQRT  (XT),  ITUP  (XT) 

Standard  functions  ( Section  4.4) i 

T8M  (  T),  TCOS  (  T),  TM  (  T),  TARCTAlf  (  T) 

ITS IS  (IT),  ITCOS  (IT),  XTLS  (IT),  ITARCTAM  (IT) 

Additional  functions  (Section  4.5) » 

TRDRGK  (  T),  TTAH(  T) 

ITROMGS  (IT),  XTTAN(XT) 

Differentiation  and  integration  (Section  4.6) t 

TDIPT(  T),  TXSTGRLI  T) 

XTDXFF(IT),  ITIMTSRL(IT) 

Miscellaneous  utilities  (Section  4.7)i 

VRMQLL,  T_I DWT_tKRO (  T),  T_IO«ITCOM8TAST(  T),  ITMIDPT(IT) , 

IVJMDIX.  IT_IDEKT_ZKRO( IT) ,  IT_ID«MT_COMSTAHT( IT) ,  WRITK_SBRIES(  T) , 
RBADJtMTKRVAL  SIRIIS(IT),  ¥RIT*_IMTEFVAI.(  I),  WRITB_IMTWWAL_SJrRIS8(  IT) 

A  brief  description  of  the  Method  for  Introduction  of  use r-de fined  functions  will  be 
given  in  Section  4.5.  Oo«»e  implementation  details  of  the  operators  and  functions  will  now 
be  discussed.  Die  recurrence  relations  are  taken  from  [16].  In  following  the  conventions 
of  Pascal-SC,  minor  differences  froa  the  indices  found  thsrs  are  due  to  our  starting  the 
series  indices  at  1  instead  of  starting  at  0.  In  each  Section,  operations  involving  the 
scalar  types  are  discussed  before  turning  to  types  TAYLOR  and  I TAYLOR.  The  Pascal-SC 
source  code  for  each  of  these  operators  and  functions,  as  well  as  for  utilities  which  they 
require,  is  included  as  Appendix  C  to  this  report. 


quit*  straightforward 


Addition i 

♦  T,  A  ♦  T,  T  ♦  A,  *  ♦  T, 

♦  it,  a  ♦  it,  it  ♦  a,  i  *  n. 


r  ♦  *. 

1*  ♦  X, 


*  ♦  * 
it  ♦  it 


Subtraction  > 


-  T,  A  -  T.  f  -  It,  * 
-IT,  It  -  IT,  IT  -  A,  X 


T,  T  -  A,  T 
IT,  IT  -  I,  IT 


T 

IT 


Addition  and  subtract  ion  of  a  oonatant  altara  only  tho  nlM  of  a  variable,  not  tba  vetoes 
of  any  of  ita  derivatives.  Interval  constant*  only  require  that  the  appropriate  built-ia 
interval  operator  be  used.  Otherwise,  addition  or  subtraction  of  aeries  is  dsns  tore- by¬ 
te  re. 

If  U  i*  I  ♦  S,  then 

(4.1.1)  O.TCtK]  i-  P.TC(A)  ♦  e.TClK)  ,  It  -  1, . .  .,DXM. 

4.2.  Multiplication  and  Division. 

Multiplication! 

K  •  T,  T  •  A,  A  *  T,  T  •  A,  T  •  T 
A  •  IT,  IT  *  A,  I  •  IT,  IT  •  I,  IT  *  IT 

Multiplication  and  division  of  two  Taylor  aortas  is  dons  by  the  well-known  Laihais 
rule  for  the  My  lor  coefficients  of  a  product  [16). 


no  -  r  •  G,  than 

K 

(4.3.1)  O.TCW  -  £  P.TCH)  •O.TClK-I+l)  ,  X  -  1, . . .  .DIM. 

1-1 

lha  aealar  product  of  two  vectors  la  evaluated  in  Paaeal-SC  fay  tha  standard  function  acALP 
to  tha  cloaaat  floating  point  n unbar.  Pa at  aarlaa  multiplication  techniques  vara  not  uaad 
harm  faaeaaaa 

o  In  many  applications  of  *,  tha  aarlaa  for  0  la  faalng  ganaratad 
recursively.  That  la.  tha  variables  P  or  0  involve  Q  ltaalf. 
o  lha  accuracy  of  SCALP  mould  not  ha  aval labia, 

o  tha  apaad  of  SCALP,  aapacially  whan  ran  tana  are  aero,  aakea 

theea  tachniguea  laaa  attractive. 

Maltlpllcatioa  or  dlvlaion  of  a  aarlaa  by  a  conatant  la  dona  tem-by-ten.  Diviaion 
of  a  conatant  fay  a  aarlaa  la  dona  fay  generation  of  the  aarlaa  for  C/r(x). 

Diviaion i 

*  /  T,  »  /  X,  X  /  T,  T  /  A,  T  /  T 
X  /  It,  It  I  X,  I  /  It,  It  /  X,  It  /  IT 

If  0  i>  p  /  o,  than  0  *  o  -  P,  and  Lalfanla*  rule  applies i 

O.TCtll  -  P.TC[1J  /  G.TC[Hi 
for  K  -  2,..., DIM 

(4.3.3) 


If  etxg)  •  G.TCl  1]  -  0,  than  we  attempt  return  thn  oorrnot  «a»Nr  Awtvtr  possible. 
If  r( xQ)  •  r.TCl 1]  Is  slso  0,  thsn  m  can  apply  1' Hospital's  ruin  because  the  series  for 
both  F  snd  G  are  known.  0.1C(1]  “  r* (xQ)  /  G’(Xg)  »  F.TC(2}  /  G.TC[2],  If  this  quotient 
exists,  but  0  f  r'  /G’  as  functions. 

if  B  i«  F  /  6,  and  F(x0)  ■  G(xq)  "  0,  then  let 

(4.2.3)  V.ycCKl  «-  r.TCtr+lli  H.TCtK)  «-  C.TC(*H)|  K  -  1... ..DIM-1. 

Then, 

(4.2.4)  0  i-  V  /  W. 

Thus,  l'Boapltal's  rule  Is  implemented  as  a  recursive  call  to  the  division  operator  with 
operands  whose  series  length  has  been  reduced  by  one.  This  approach  would  not  be  possible 
in  a  language  which  does  not  support  recursion.  Further,  cases  in  which  the  series  for 
both  f  and  g  have  several  leading  seros  are  handled  automatically  by  the  language. 

L'Hoepltal's  rule  is  applied  in  a  similar  manner  when  a  constant  quotient  or  li visor 
is  equal  to  sero. 

4.3.  Power  Operators,  the  power  operator  **  defined  by  F  **  G  •  F®  is  not 
standard  in  Pascal  or  Pascal-SC,  but  can  be  implemented  in  the  latter  for  data  types  for 
which  it  is  meaningful  by  the  use  of  the  operator  concept.  Oodlng  of  •*  is  simplified  by 
the  Introduction  of  a  set  of  basic  power  functions.  These  are  implemented  separately 
o  for  uniformity  with  Pascal-SC  which  provides  these  functions 
for  standard  data  types, 

o  to  provide  tighter  bounds  for  interval  operands,  and 
o  for  efficiency. 
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4.3.1.  Special  Knot  Function*  This  Mt  of  functions  consists  of  ths 
square,  square  root,  and  natural  exponential  function  of  variables  of  types  TAYLOR  and 
ITAYLORt 


T SfiR  (  T),  n0KT  (  T).  TEXT  (  T) 

IT8QR  (IT),  IT8QRT  (XT),  ITIZP  (IT) 

These  functions  are  called  by  the  operator  **  when  appropriate.  For  exaaple,  if  X  is  of 
type  ITAYLOR,  then  both  X  **  2  and  X  **  INTPT(2.0)  are  actually  performed  by  a  call  to 
XT8QR  (X).  The  use  of  this  function  rather  than  X  *  X  is  important  in  interval 
computations,  since,  for  exaaple,  f— 1 ,— 1 1 a  "  (0,1)  while  t-1,-1)  •  1-1,-IJ  ■  [-1,+1J. 
Further,  the  squaring  functions  T8QR  and  IT  SQR  are  twice  as  fast  as  the  multiplication  Y  * 
Y  for  variables  of  the  corresponding  types. 

nie  recurrence  relations  to  generate  the  series  terms  for  these  functions  can  be 
derived  easily  using  Lsibnls*  rule.  The  squares  of  real  and  Interval  Taylor  series  are 
cosputed  as  follows. 

If  U  ■  sqr(F),  then  Leibnls'  rule  for  a  product  can  be  shortened  toi 

For  K  -  1,  ...,  DIM, 

X  DIV  2 

(4.3.1)  O.TCtXj  -  £  F.TC(Il*F.TC(X-I-MJ» 

I  ■  1 

if  X  is  odd,  then  O.TCIR]  -  O.TC[K)  +  SQR  (F.TC((XVl)/2] ) . 

The  inner  product  contains  only  TR0MC(X/2)  terms.  If  F  is  of  type  ITAYLOR  and  Includes 
negative  numbers,  then  XTSQR  (F)  provides  tighter  bounds  than  does  F  *  F.  The  SQR 
functions  are  nested  T8QR  and  IT8QR  to  indicate  the  type  of  operand  accepted  and  value 
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returned 


A  similar  function  mi  written  for  COBS.  Its  MMitiou  had  length  YMMC(X/3)  bat 
thsy  war*  nested  to  yisld  •  cost  proportional  to  DIM*.  COM  la  not  included  In  the  library 
because  r  •  SQM  (F)  is  faster. 

Tbs  functions  in  the  next  set  calculate  square  roots  of  real  and  interval  Taylor 
variables. 

If  a  i“  sqrt  (r),  then  0  *  0  -  P.  The  algorithm  runs  as  fellows  i 

O.TCll]  i*  SQRT  <F.TC[1])» 

O.TC(2)  i-  P.TC[2]  /  <2  *  O.TC(1])l 
for  K  -  3,  ....  DIM, 

(4.3.2) 

K  DXV  2 

PROD  t  0.TC[I1*0.TC[X-I-M] » 

1-2 

if  X  is  Odd,  then  noo  t-  PNOO  *  SQM  (U.TC[(XVl)/2))» 

O.TCIX]  i-  (F.TC(X)  -  PMDD)  /  (2  *  O.TCtl]). 

If  F(x„)  -  0,  and  T  is  not  a  constant  saries,  then  SQMT  (F)  cannot  be  oaapotad  unless  P  is 
the  constant  0,  becauee  F'trg)  is  not  ds fined.  The  SQMT  functions  are  nsaied  TSQMT  and 
IT8QRT  to  indicate  the  type  of  operand  accepted  and  value  returned. 

The  natural  expontlal  functions  (base  -  e)  are  now  defined  for  types  TAYLOR  and 
ITAXLOR. 

If  O  m  (F),  then  0*  -  0  *  F*.  This  gives  the  algorithsi 

O.TC(1]  l-  OF  (F.TC(IJ); 
for  X  -  2,..., DIN, 

(4.3.3) 

K-1 

O.TCIX)  i-  (  I  O.TC(I)»F.TC[X-I-MI*(X-I)  )  /  (X-1); 

1-1 

Although  this  forsnila  would  appear  slightly  slnplsr  with  the  change  of  index  J  -  X  *  X,  it 


mi  Implemented  in  this  way  so  that  ths  O.TC  tons  rowain  stationary  in  tho  inner  product 
as  X  increases.  Thus,  only  ths  vector  F.TC  weeds  to  be  "reversed",  the  tfXP  functions  are 
named  TtXV  and  XTSXP  to  indicate  the  type  of  operand  aeoepted  and  value  jreturaed. 

4.3.2.  The  Operator  •*.  The  f sally  of  power  operators  **  seems  to  be  the 
most  difficult  to  implament  as  suggested  by  (5] .  Hone  of  the  operators  are  especially 
difficult,  but  there  are  aany  einor  details  to  be  considered.  Ths  implementation  of  **  for 
types  TAYLOR  and  ITAYLOR  is  based  on  ths  standard  power  functions  above,  and  the  power 
operators  **  for  the  scalar  types  INTEGER  (X),  RIAL  (R),  and  INTERVAL  (X). 

Scalar  Powers: 

X  **  X,  R  **  X,  X  **  R,  R  ••  R 

I  ••  X,  X  •*  X,  I  **  I 

Integer  powers  are  implemented  using  repeated  squaring.  Real  and  interval  powers 
which  fit  no  special  case  are  computed  by  P  **  G  “  EXP  (G  *  LN( F) ) .  He  have  not  attempted 
optimal  implementations  of  the  scalar  power  operators  because  it  is  hoped  that  they  will  be 
provided  as  standard  operators  in  a  later  release  of  Pascal-SC,  an  approach  that  is 
especially  attractive  for  Interval  operands  because  ths  Interpreter  hides  information  from 
programmers  which  can  be  used  for  correctly  directed  roundings. 

Real  and  Interval  Taylor  Powers: 

X  ##  t,  T  **  X,  R  **  T,  T  **  R,  x  *•  T 

X  «*  XT,  IT  •*  X,  I  •*  IT,  IT  •»  I,  IT  *•  IT 

The  power  operators  for  a  constant  to  a  variable  power  follow  the  pattern  of  TEXP  or 
ITKXP,  as  appropriste.  A  series  which  represents  a  constant  (only  its  first  term  is  non¬ 
zero)  is  handled  as  a  special  case  for  accuracy  (especially  for  interval  series)  and  for 


efficiency 


Ih«  operators  T  *•  X  and  XT  ••  X  taka  oars  to  return  the  correct  answer  ebon sear  that 
is  possible  and  to  produce  an  appropriate  error  nos sage  when  it  is  not  possible,  the 
resolution  of  various  eases  is  shown  in  Table  4.1. 

Exponent! 

Base.TC  •  0 

Base. TCP]  -  0 

Base. TCP]  <>  0 


Consider  a  aeries  whose  first  tern  is  sero,  but  which  has  other  terns  which  are  non- 
sero.  Raising  such  a  aeries  to  a  negative  power  is  undefined  becauee  it  is  equivalent  to 
dividing  by  sero.  but  raising  such  a  series  to  the  power  0  defines  a  function  which  is 
identically  equal  to  1,  except  for  a  raswvahle  singularity  at  t  •  Xq.  hence  it  is 
appropriate  to  give  1  as  the  answer.  Raising  the  series  to  a  positive  Integer  power  is 
lnplenented  by  repeated  squaring  because  the  recurrence  relation  which  is  nost  often  used 
[20]  requires  division  by  BASt.TCP] ,  which  is  sero. 

The  special  cases  of  an  exponent  equal  to  1  cr  2  are  singled  out  for  Individual 
treetnent  in  order  to  achieve  the  eaxleo*  possible  accuracy  (especially  when  the  base  is  an 
Interval  aeries)  and  for  efficiency. 

Except  in  the  special  cases  shown  in  the  table,  if  0  ■  P  •*  X,  where  ■  is  of  type 


>  2 


<  0 


Ohdef . 


■  0 


Ondef. 


I 

8QR  (BAH)  I  By  suit.  Ondef. 

I _ 

I 

8QR  (BABB)  I  By  recurrence 

I 


Table  4.1.  Resolution  of  Cases  for  **■ 


XRTB8BR,  then  r  *  O'  •  B  •  0  *  P*.  This  gives  the  algorithm 


O.TC  I"  0 t 


D.TC[  1 J  i-  P.TCMJ  **  I; 

(4.3.4)  fftr  X  i*  2  to  DIM, 

K-1 

O.TC[K]  i-  (  Z  (E*(K-I>  -  I  ♦  1)  «  U.TC[I)  *  P.YC[K-X-H])  /  ((X-1)«F.TC[I) >. 

1-1 

The  integer  exponent  appears  in  tha  racurranea  only  as  a  aultipliar.  Hence  tha  speed  of 
this  algorithm  is  nearly  Independent  of  the  else  of  tha  exponent.  That  is  why  this 
algorithm  is  preferred  to  repeated  squaring. 

The  operators  T  **  R  and  IT  •*  I  are  similar  to  T  *•  X  and  IT  *•  X,  respectively, 
except  that  the  additional  special  cases  of  an  exponent  equal  to  1/2  or  to  an  integer  are 
handled. 

The  operators  T  •*  T  and  IT  **  IT  are  Included  primarily  for  completeness;  the  authors 
have  never  seen  a  differential  equation  containing  a  variable  to  a  variable  power,  for 
example.  Perhaps  any  such  problems  which  arise  are  at  once  simplified  by  logarithmic 
differentiation.  With  the  tools  described  here,  it  may  be  advantaglous  to  attack  the 
problem  in  its  original  form. 

Within  the  oparatore  T  •*  T  and  IT  •*  IT,  the  cases  in  which  either  the  bane  or  the 
exponent  aeries  represent  a  constant  function  are  treated  as  special  for  reasons  of 
accuracy  and  efficiency.  Otherwise,  P  ••  G  -  EXP  (G  *  LM  (P)),  using  TEXP  and  TUI  or  ITEXP 
and  ITU>,  as  appropriate. 

4.4.  Standard  Functions.  There  are  many  useful  library  functions  which  could  be 
provided,  we  have  chosen  to  implement  the  functions  which  are  built  into  Pascal-9C  for  the 
standard  scalar  data  types,  and  a  few  others.  Additional  functions  can  be  added  as  they 
are  needed  by  following  the  models  provided  by  this  report.  In  addition  to  tha  standard 
power  functions  of  Section  4.3.1  (which  Include  EXP  and  XEXP),  other  standard  functions 
implemented  for  types  TAYLOR  and  ITAYLOR  are 
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TSIN  (  T),  TC08  <  T),  TLH  (  T) ,  TMtCTXN  (  T) 

ITS  IN  (IV),  ITCOS  (IT),  ITLN  (IT),  ITANCTMI  (IT) 

If  0  -  sin  (F)  tai  T  «  cot  (F),  than  0*  -  V  *  F'  and  V*  -  -  0  *  F*. 

0.TCI1)  l-  8 in  (F.TC(ll)l  V.TC(I)  I-  COS  (F.TC(IJ)l 
for  X  i-  2,.,.,OIM, 

X 

(4.4.1)  O.TCtX]  t-  (  I  V.TClIl  •  F.TCEX-X+1]  *  (X-I)  )  /  (X-1)» 

1-2 

X 

V.TCtX]  »-  -  (  I  O.TC(I)  *  F.TClX-I+1)  *  (X-I)  )  /  (X-1). 

1-2 

Tha  SIN  and  COS  function*  ara  naaad  TSIN  and  TCOS  or  ITS IN  and  ITCOS  to  lndicata  tha 
typo  of  oparand  which  thay  accept  and  valua  they  r a turn,  since  the  aerie*  for  SIN  and  COS 
are  always  computed  together,  the  library  also  contains  procedures  T_8lN_cos  and  IT_8IN_C06 
which  return  both  the  SIN  and  COS  of  variables  of  type  TXtu»  and  rrxnOR,  respectively,  in 
the  mum  call. 

If  U  «-  In  (F),  then  O'  •  F  -  F’. 

O.TCll]  i-  LN  (F.TCl  1) ) t 

(4.4.2)  for  K  t-  2 . DIN, 

X-1 

O.TCtX]  i-  (F.TCtK)  -  (  I  O.TCtl)  •  F.TClX-Itl]  *  (1-1)  )  /  (K-1))  /  F.TC[1). 

1-2 

There  Is  a  eisprint  In  this  rscurrencs  relation  in  ([16],  p.  42),  but  its  implementation  is 
straightforward. 


If  O  »-  Cretan  (F)  and  V  i-  1  /  ( 1  ♦  F2) ,  than  0*  -  V  *  F*. 

V  »«  1  /  {  1  ♦  8QR  (F)  )i 
0.TCI1]  ARCTAN  (  F.TCI1]  )> 

(4.4.3)  for  X  «-  3,..., DIM, 

K 

U.TC (X)  «-  (  I  V.TC(I)  *  F.TC(X-I+1)  *  (X-I)  )  /  (X-1). 

1-2 

Sinca  tha  aarlaa  for  the  Rungs  function  V(F)  ((7),  p.  78)  la  required  to  compute  the 
aerlea  for  arctan(F),  functiona  TRUNGE  and  ITRUNCT  are  included  in  the  library  along  with 
the  functiona  TARCTAN  and  ITARCTAN . 

4.5.  Dear  Defined  Functiona.  If  a  programmer  require*  an  operation  or  a  function 
which  la  not  included  in  thle  report,  tha  requirement  can  be  net  either  by  a  composition  of 
operators  and  functions  which  are  already  provided,  or  by  a  careful  derivation  of  the 
necessary  recurrence  relations  following  the  models  in  this  report.  For  example,  the 
tangent  functiona  TTAN(  T)  and  ITTAN(IT)  are  implemented  essentially  by 

(4.5.1)  TTAN(  T)  1-  TSIN(  T)  /  TCOS(  T),  ITT  AN  ( IT )  »-  IT8IN(  IT)  /  ITCOS(IT), 

respectively  in  the  aet  of  additional  functiona  provided  in  the  library.  The  tangent 
functions  can  aleo  be  implemented  directly  by  recurrence  relations,  using  the  fact  that  y  - 
tan(x)  satisfies  the  differential  aquation 

(4.5.2)  y’  -  1  ♦  y2.  y(*0>  -  tan(x0), 

(12),  (13).  Thus,  for  U  t-  TTANI  T),  for  example, 

(4.5.3)  U.TC ( 1 1  «-  TAN(T.TC(1] ),  O.TC(2)  i-  (1  ♦  SQR(T.TC[ 1 ) ) )  *  T.T. 
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The  succeeding  coefficients  can  be  obtained  in  s  slnple  way  froe  the  recurrence  relation 
<4.3.1)  for  HQk(  t;,  and  ITCAN(tT)  is  computed  analogously. 

the  tangs  function  f(x)  -  1  /  <1  ♦  x2),  which  la  an  auxiliary  function  for  the  series 
expansion  of  the  arctangent,  la  leplsesnted  in  the  library  by 

0  1-  T8QR(IT) I 

(4.5.4)  O.TC[t]  i-  1  ♦  0.9C{1) > 

TRONGE  i-1/O. 

ITRUUGK  is  coaputed  slellarly. 

4.6.  Differentiation  and  Integration.  Functions  which  return  the  results  of  tsre- 
by-terw  differentiation  and  integration  of  TAYLOR  and  I TAYLOR  series  are  also  provided. 

For  series  with  1  <  LENGTH  <_  DIM,  differentiation  decreases  the  length  of  the  series  by 
one: 

(4.6.1)  O.TCtK]  i-  T.TCIK  ♦  1]  *  RATIO  /  (X  ♦  1),  X  -  1, . . . ,T. LENGTH  -  1, 

where  RATIO  -  1  /  T.T  if  0  -  TOIFF  <T>,  and  RATIO  -  1  /  INTPT  (T.T)  if  0  -  ITDIFF  (T). 
Integration  results  in  a  aeries  with  its  first  coefficient  set  to  0  and  its  length 
increased  by  onei 

(4.6.2)  U.TC(K]  t—  T.TCIX  -  1]  *  RATIO  /  (K  -  1),  X  -  2, ... ,T. LENGTH, 

with  O.TC[1)  -  o  end  RATIO  -  T.T  for  0  -  TINTGRL(T),  while  O.TC(1)  -  INTPT(O)  and  RATIO  - 
INTPT(T.T)  for  0  -  ITINTGRL(T).  The  result  of  Integration  of  a  series  of  length  DIN  will 
be  truncated  to  length  DIM. 
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4.7.  W^egUgMow  Ptilltiee.  Some  useful  function*  and  procedural  axe  provided 
tor  convenience.  These  ere  the  transfer  function  TMIDPT  (IT),  the  special  functions 
VMOLL,  XvmOLL,  the  comparison  functions  T__IDENT_ZERO,  T_I DENT_CON8TAWT ,  rT_IDUrr_**Jto, 
IT_ID*NT_COH8TART,  and  the  input/output  procedures  WRIT«_I»T*RVAL  (I),  WRITEJBIRIE8  (T), 
READ_INTEAVAL_8ERISS  (IT),  OTUTI_INT«RVAL_8ERIES  (IT).  The  purpose*  Of  nost  Of  the** 
utilities  are  indicated  by  their  names. 

The  transfer  function  TMIDFT  (IT)  forms  a  TAYLOR  series  free  a  series  of  type 
I TAYLOR.  The  coefficients  of  the  result  series  are  the  eidpolnts  of  the  corresponding 
coefficients  of  the  interval  series. 

The  parameterless  functions  VRNULL,  rvRHULL  yield  sero  real  and  interval  vectors, 
respectively,  of  length  DIM.  They  are  standard  Paacal-SC  functions  123] . 

The  comparison  functions  ylsld  the  BOOLEAN  value  TRUE  if  their  argiseent  satisfies  the 
stated  condition  (the  series  la  identically  equal  to  sero  or  a  constant),  otherwise,  FALSE. 

The  input/output  procedures  are  also  self-explanatory.  The  procedure  mXTE_XMTEKVAL 
la  included,  since  the  standard  Fascal-SC  procedure  IWRITE  only  prints  the  digits  of  the 
lower  and  upper  endpoints  of  intervals  which  agree  up  to  the  last  [23] .  WRITE_INTEKVAL, 
however,  prints  all  digits  of  each  endpoint. 

5.  The  initial-value  problem  for  ordinary  differential  aquations.  Taylor  series 
■ethods  for  the  mmerlcal  solution  of  initial-value  problem*  for  systems  of  ordinary 
differential  equations  have  been  studied  by  aany  authors  (see  [6]  or  [13]  for  summaries), 
and  have  been  used  for  applications  such  as  dynamics  and  parameter  identification.  Each 
component  of  the  solution  of 

(3.1)  y4'  «  fjfx.y),  yt(x„)  -  yi0,  i  -  1,...,m, 

is  expressed  as  a  Taylor  aeries  expanded  at  i  •  xg  using  recurrence  relations  derived  from 
the  differential  equation.  Various  error  control  strategies  have  been  employed.  The 
strategy  of  analysing  the  radius  of  convergence  of  each  component  series  has  ths  desirable 
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•Ida  effect  of  producing  such  analytic  information  aa  tha  location  and  orders  of  tha 
alngularitiaa  in  tha  solution.  Onca  tha  radiua  of  oonvsrgsnoe  ia  known,  a  a ta paisa  can  bo 
choaon  which  ia  as  largo  as  possible  subject  to  error  control  and  stability  constraints. 
Than  each  component  of  the  solution  ia  extended  by  analytic  continuation  and  tha  process  is 
repeated  at  the  next  Integration  step,  this  algorithm  la  diaeusaed  in  greater  detail  la 
[6]. 

A  program  RDIQ_SOLV  for  solving  aquation  (S.1)  is  given  as  Appendix  A  of  this 
report,  the  program  is  written  for  the  case  m  -  1,  but  can  be  modified  easily  to  handle 
systems  of  several  equations,  tha  variables  f  and  YPRXKS  are  declared  to  be  of  type 
TAYLOR,  and  the  equation  ia  written  in  a  natural  way.  To  solve  a  different  equation,  it  is 
only  necessary 

o  to  change  the  line  in  RMQ^SOLV  which  contains  the  differential  aquation, 
o  to  copy  from  the  library  into  the  source  program  any  operators  or  functions 
required  by  the  new  differential  equation. 

Because  the  differential  equation  is  written  using  the  types  and  operators  discussed  in  the 
preceding  Sections,  the  needed  recurrence  relations  are  implemented  by  the  Pascal-SC 
compiler  and  need  not  be  derived  explicitly  by  the  user. 

The  program  prints  the  series  terms,  extends  the  solution  by  analytic  continuation  to 
compute  the  initial  condition  at  the  next  step,  and  repeats  tha  process.  The  program 
RD*Q__SOI/V  in  Appendix  a  is  intentionally  slepls  to  illustrsts  tha  use  of  tha  Taylor 
operators  and  to  explore  the  stability  of  tha  ssrlss  gansratlon.  It  would  require  an  error 
control  mechanism  in  ordar  to  bs  of  practical  valua  for  the  solution  of  initial  veins 
problems.  Rlthsr  scalar  (6]  or  lntsrval  [12]  error  control  techniques  can  be  used. 

The  program  XDBQ^SOLVS  Hated  aa  Appendix  B  of  this  report  ccaputes  interval-valued 
approximate  solutions  to  aquation  (5.1)  for  tha  casa  m  ■  1,  but  can  be  modified  for  systems 
of  aavaral  aquations.  It  dlffera  from  the  program  RD*^_SOLV  only  in  that 

1)  the  variables  Y  and  YPRXNS  are  of  type  XTAYLOR  instead  of  type  TAYLOR,  and 
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ii)  additional  coda  ha*  baan  ad  dad  to  Monitor  tha  relative  arror 
introduced  by  lnatability  in  tha  aarla*  generation  proceae. 

Thaaa  two  program*  ara  designed  to  aerv*  aa  example*  of  on*  way  in  which  tha  tool*  of 
thia  report  can  b*  uaed.  Th*y  ara  simple,  menu-driven  profraaa  which  allow  direct  uaar 
intervention  at  aach  integration  atep.  By  observation  of  the  outcoam  of  each  atop,  the 
uaar  can  experiment  with  arror  control  strategies.  A  Dear  Manual  containing  nor*  detailed 
lnat ruction*  tor  uaing  thaaa  program*  ia  included  aa  Appendix  D  of  thia  report. 

Tha  bound*  coaputad  by  ID*{J_80LV  are  for  tha  interval-valued  Taylor  polynoadal 
(3.3).  They  are  not  global  arror  bound*  for  tha  aolutlon  of  tha  differential  aquation. 
Global  arror  bound*  are  readily  coaiputabl*  uaing  interval  remainder  tana  (aa*  [12]),  but. 
for  simplicity,  tha  prograaa  given  hare  contain  no  arror  bounding  or  control  atratagy. 

6.  An  application i  Stability  of  aerie*  generation,  in  thia  auction,  w*  praaant  an 
example  which  uaaa  tha  Taylor  and  interval  Taylor  operator*.  Thia  example  waa  choaan 
bacauaa  it  illuatrataa  the  uaaa  of  theaa  operator*  and  becauae  it  addreaaaa  the  laaua  of 
atabillty  in  tha  generation  of  the  aerie*,  tti*  latter  laaua  la  central  in  ahowlng  that 
Taylor  aarlea  nethoda  ara  reliable  for  practical  computation*. 

A  numerical  computation  ia  aaid  to  be  unstable  if  it*  relative  arror  grow*  without 
bound  aa  the  computation  proceed*.  It  ia  poaaibl*  that  the  recurrence  relation*  being  uaed 
might  be  unatable,  although  lnatability  haa  never  baen  obaervad  in  practice.  Thia  example 
uaaa  the  Taylor  and  interval  Taylor  operator*  to  explore  the  atabillty  of  the  recurrence 
relation*  in  one  application.  In  thia  example,  there  la  lnatability  in  the  generation  of 
the  term*  of  the  aarlea,  but  that  doe*  not  aerloualy  affect  the  accuracy  of  the  aerie* 
auamatlon.  Tha  atabillty  of  the  recurrence  relatione  in  other  application*  can  be  handled 
almllarly. 
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Consider  th«  initial  value  problem 


(6.1)  y'  -  y2,  y(0)  -  t, 

whose  solution  Is  y(x)  «•  1  /  (1  -  x).  A  program  ! RDCQ_80LV )  for  solving  equation  (6.1) 

using  the  Taylor  function  T8QR  is  given  as  Appendix  A  of  this  report.  The  effect  of 

program  RD*Q_SOLV  la  to  generate  the  normalised  Taylor  coefficients  (4.3.1)  of  the  solution 
recursively.  This  recurrence  is  accomplished  automatically  by  the  Taylor  function  T8QR  in 
the  statement  TPRIMX  :■  T8QR  (y).  In  this  case,  the  same  solution  Is  obtained  if  T  **  2  or 
r  ‘  r  is  used  Instead  of  TSQK  (Y)i  however,  the  use  of  *  •*  2  requires  the  compilation  of 
much  sure  cose,  while  *  *  »  Is  not  as  fast  as  T8QR  (Y). 

Ns  wish  to  explore  the  stability  of  the  recurrence  relation  (6.2).  This  Issue  is 
separate  from  the  issue  of  the  stability  of  Taylor  series  methods  for  solving  Initial  value 
problems.  If  an  infinite  Taylor  series  were  used,  the  method  would  be  A- stable,  but  when  a 

truncated  series  is  used,  the  region  of  stability  is  bounded.  Stetter  [21]  showed  that  the 

region  of  stability  for  truncated  Taylor  series  methods  is  the  same  as  that  for  related 
Runge-Kutta  methods.  The  real  interval  Of  stability  is  relatively  large  here  because  long 
series  are  used.  Tor  example,  the  real  intervals  of  stability  are  [>8.85,  0]  and 
(-16.29,  0],  respectively,  if  DIM  "  20  and  40  terns  of  the  series  are  used. 

Ns  will  outline  the  theoretical  analysis  of  the  stability  or  recurrence  (6.2).  A  more 
complete  discussion  appears  in  [6].  Let  U(K)  denote  the  actual  and  Y(JC)  denote  the 
computed  normalised  Teylor  coefficients,  bat  Y(1)  ■  0(1)  (1  +  e)  »  (1  ♦  c)  from  (6.1). 

Then  0(K)  -  h*”1,  and 

(6.3)  Y(K)  -  0(K)  (1  ♦  e)K, 

so  the  series  generation  la  unstable.  However,  the  summation  of  the  series  is  unaffected 
by  this  Instability  since 
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DIM  DIM  .  . 

y(X,)  -  E  Y(K)  -  (1  ♦  e)  E  (h(1  t>) 

K-1  KfI 

(6.3) 

-  (1  ♦  c)  y(h<  1  +  e>)  ♦  Olh01*1) 

-(It  e)(y(h)  ♦  y*  (C)h)  ♦  0(hDU,)f  h  <  C  <  h(1  ♦  e). 

That  is,  ths  Instability  in  the  ssrlss  generation  is  equivalent  to  s  small  srror  in  ths 
point  at  which  ths  sariss  is  evaluated.  This  is  because  EgY(K)  is  a  convergent  ssrlss,  so 
ths  tons  for  which  instability  causas  ths  relative  srror  to  ba  largest  ara  themselves  wary 
small. 

This  suggests  using  interval  arithmetic  to  kaap  track  of  tha  potantlal  growth  in  ths 
sariss.  Tha  program  1DBQJKJLV  11 a tad  as  Appendix  B  of  this  raport  doss  so. 

By  dsclaring  T  and  TMIIM8  to  bs  of  typa  itaylor,  tha  statamant  YPRXMK  i-  XTBQR(T) 
lnvokas  tha  function  ITSQR  for  Interval  valued  sariss  to  ganarata  Interval  norma 11  rad 
Taylor  coofflclants  according  tha  racurranca  relation  (4.3. t).  Tha  lengths  of  successive 
coaff iclanta  maasura  tha  stability  of  tha  racurranca.  Table  6.1  shows  tha  Interval  valued 
series  solution  of  agnation  (6.1)  for  DIM  -  15,  y(x„)  -  y(0)  -  [0.99  ,  1.01] 

<c  -  0.01),  and  h  -  0.5.  Tha  computed  instability  is  equal  to 

length  (T.TC(K]  ) 

<6'5>  ‘computad  ’  aidpolnt  (T>TC(K, )  ' 

a  measure  of  tha  relative  error  in  Y.TC(K)  which  appears  to  grow  as  K  increases.  Ths 
theoretical  Instability  is  equal  to 

*«•«>  ^Theoretical  “  (1  ♦  e)** 

Tabla  6.1  shows  that  these  two  values  are  vary  close,  and  that  ths  theoretical  bound  is 
slightly  larger  than  tha  actual  bound,  as  it  should  be.  Tha  interval  estimate  for  y(0.5) 
agrees  wall  with  tha  interval  [y(0.46),  y( 0 . 54 ) ]  *[1.8518,  2.1739]. 
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MAD  nmtKVXL  INITIAL  CONDITIONS  XO,  y(XO)  I 


INITIAL  CONDITIONS  AT  XO  -  [  0.00000X400,  0.00000X4-00], 

TO  -  {  9.90000X-01,  1.01000X400]. 

INTXR  STXMI8X  X  -  XOi  O.S 
Computing  MTlM  t*nu  ... 


step 

Left 

Right 

Coaputed 

Theoretical 

Xndpoint 

Xndpoint 

Instability 

Instability 

1 

l 

9.900008-01, 

1.010008400] 

1.0108400 

1.0108400 

2 

t 

4.900508-01, 

5.100508*01} 

1.0208400 

1.0208400 

3 

l 

2. 425758-01, 

2.575758-01] 

1.0308400 

1.0308400 

4 

t 

1.200758-01, 

1.300768-01} 

1.040X400 

1.041X400 

5 

f 

5. 943698-02, 

6.568818-02] 

1.0508400 

1.0518400 

6 

l 

2.942138-02, 

3.317258-02] 

1.060X400 

1.0628400 

7 

t 

1.456358-02, 

1.675218-02] 

1.0708400 

1.072X400 

8 

t 

7.20S948-03, 

8.459828-03} 

1.0808400 

1 .0838400 

9 

{ 

3.568438-03, 

4.272218-03] 

1.0908400 

1.094X400 

10 

I 

1.766378-03, 

2.157478-03] 

1.1008400 

1.1058400 

11 

[ 

8. 743548-04, 

1.089528-03] 

1.110X400 

1.116X400 

12 

l 

4.328058-04, 

5. 502088-04] 

1.1198400 

1.1278400 

13 

t 

2. 142398-04, 

2.778558-04] 

1.1298400 

1.138X400 

14 

I 

1.060488-04, 

1.403178-04] 

1.139X400 

1.1498400 

15 

t 

5.249388-05, 

7.085998-05] 

1.1498400 

1.1618400 

TUX  TALDX  AT  X  -  (  5.000001-01,  5.00000X-01] 

IS  T  -  l  1.96034X400,  2.04033X400], 

Table  6.1.  Interval  bounds  for  instability. 


7.  Inol— isntatlon  details.  Tbs  software  described  in  this  report  was  created  and 
tested  using  the  Raacal-sC  coepller  developed  at  the  University  of  Karlsruhe  for  the  Zilog 
HCZ-1  coaputer  with  the  MO  2.06  operating  ays  tee.  No  other  clalaus  of  correctness  or 
usability  are  made. 
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APPENDIX  A 


PROGRAM  RDMQ_SOLVE  (INPUT,  DATA,  OUTPUT )  I 

(*  SOLVE  A  FIRST  ORDER  DIFFERENTIAL  EQUATION)  T'  -  SO*  (V)  *> 

COMET  DIM  -  30| 

TYPE  DIMTTPE  -  1 . .DIM; 

R VECTOR  -  ARRAY  [DIMTTPE]  OP  REAL  / 

TAYLOR  -  RECORD  LENGTH  t  DINTYPI) 

T  t  REAL; 

1C  i  R  VECTOR  EMDl 

CHOICE  *  1..3) 


FLAG 
I,  INI 

X,  Y,  TP  RIME 
OATA 


i  CHOICE! 

I  DIMTTPE  | 
t  TAYLOR) 
t  TEXT » 


FUNCTION  VRNULL  s  EVBCTOR; 

VAR  It  DIMTTPE  I  Ut  RVECTOR) 

BEGIN 

FOR  I  I-  1  TO  DIM  DO  DflJ  t-  0.0/ 

VRNULL  1-  U 

END)  <*  FUNCTION  VRNULL  *) 

FUNCTION  TSQR  (Tl  TAYLOR)  «  TAYLOR) 

(•  taqulrMl  VRNULL,  SCALP,  SQR  *) 

VAR  I,  J,  K,  HALF)  DIMTTPE) 

X,  Yt  RVECTOR) 

U  i  TAYLOR) 

BEGIN 

X  t-  VRNULL)  Y  I-  VRNULL) 

O. LENGTH  i«  T. LENGTH) 

U.T  t-  T.T) 

U.TC  i-  VRNULL) 

U.TCtl)  i-  SQR  (T.ICH1  )) 
lit]  x-  T.TCtl]) 

FOR  I  I*  I  TO  U. LENGTH  DO 
BEGIN 

X[K]  t-T.TC(KJ) 

HALF  »•  X  DIV  2) 

FOR  J  l-  1  TO  HALF  DO 
BEGIN 

I  t-  K  -  J  ♦  1) 

YtJj  «-  T.TC[I] 

END)  (•  FOR  J  *) 

U.TC tH)  «-  2.0  •  SCALP  (  X,  Y,  0>) 

IF  K  MOO  2  -  1  THEN 
BEGIN 

HALF  )•  HALF  4  1) 

U.TCtKli-  U.TC (X]  +  SQR  (T.TClHALF] ) 
END  (•  IF  *) 

END)  (•  FOR  E  ») 

TSQR  t-  U 

END)  (•  FUNCTION  TSQR  (TAYLOR)  *) 


(*  TSQR(T)  •) 


•“'v.-y  ■  ;  ■ 
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FUNCTION  MENU  CHOICE  :  CHOICE  ) 


VAR  I:  INTEGER) 
BEGIN 

WRITELN) 

WRITELN  CENTER: 

1 

-  GIVE  NEW 

INITIAL  CONDITIONS')) 

WRITELN  (' 

2 

-  CONTINUE 

EXTENDING  THE  SOLUTION '  ) ) 

WRITELN  ( ' 

3 

-  STOP')) 

READ  (I); 

IF  {(I  >-  3)  OR 

(I  <- 

0))  THEN  I 

3; 

MENU_CHOICE  :-  I 

END)  (*  FUNCTION  HENU_CHOICE  •) 

PROCEDURE  PRNT_TAY_COEF  (Y:  TAYLOR;  INDEX:  DIMTYPB); 

BEGIN 

WRITELN  CYC,  INDEX:  5,  ')  «  ',  Y.TC  [INDEX]) 

END)  (*  PROCEDURE  PRNT_TAY_COEF  *) 

FUNCTION  SUM  (VAR  A:  RVECTOR?  DIM,  ROUND:  INTEGER)  :  REAL; 

EXTERNAL  480) 

BEGIN  <*MAIN  PROGRAM  RDEQjSOLVE*  ) 

(* .  INITIALIZE  ♦) 

FLAG  2; 

X.  LENGTH  DIM) 

Y.  LENGTH  :»  DIM) 

RESET  (DATA); 

WHILE  FLAG  <-  2  DO  (*  LOOP  FOR  NEW  INITIAL  CONDITIONS  •) 

BEGIN 
FLAG  :«  2) 

X. TC  :-  VRNULL) 

Y. TC  :-  VRNULL) 

WRITELN  ('READ  REAL  INITIAL  CONDITIONS  X0,  Y(X0);’)> 

READ  (DATA,  X.TC[1]))  READ  (DATA,  Y.TC[1))> 

WRITELN;  WRITELN) 

WRITELN  ('INITIAL  CONDITIONS  AT  X0  -  ',  X.TCfl], 

WRITELN  ('  Y0  «  Y.TC[1),  ’.')) 

WHILE  FLAG  -  2  DO  (*  LOOP  FOR  ANALYTIC  CONTINUATION  *) 

BEGIN 

(* . .  READ  STEP  SIZE  •> 

WRITELN  ('ENTER  STEPSIZE  X  -  X0:  ')) 

READ  (X.T) ) 

Y.T  :«  X.T) 

WRITELN  ('Confuting  series  terms  ...'); 

FOR  I  :-  2  TO  DIM  DO  (•  LOOP  FOR  SERIES  GENERATION  *) 

BEGIN 

(*  YOUR  FIRST  ORDER  DIFFERENTIAL  EQUATION  GOES  HERE:  *) 

YPRIME  TSQR  (Y); 

IMl  :-  I  -  1) 

Y.TClI]  :»  YPRIME. TC [IMl]  *  Y.T  /  IMl) 

END)  ( *FOR* ) 

(* .  PRINT  TABLE  *) 

WRITELN)  WRITELN) 

WRITELN  ('THE  TAYLOR  COEFFICIENTS  OF  Y  ARE:')) 

FOR  I  1  TO  DIM  DO  PRNT_TAY  COEF  (Y,  I); 
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(* .  PERFORM  THE  ANALYTIC  CONTINUATION  *) 

Y.TC [ 1 ]  SUM  (Y.TC,  DIM,  0)» 

FOR  Is-  2  TO  DIM  DO  Y.TC[I]  0.0) 

X.TC  [1]  :«  X.TC[1]  +  X.Tl 
WRITELN; 

WRITELN  ('THE  VALUE  AT  X  -  X.TC[1])I 

WRITELN  ('  IS  Y  »  Y.TC[1), 

FLAG  M£NU_CHOICE 
END  ( *WHILE* ) 

END  (*WHILE*) 

END.  (*  MAIN  PROGRAM  RDEfi_SOLVE  *) 


* 


« 


-3 

i 


I 


I 
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APPENDIX  a 


PROGRAM  ID*Q_SOLVX  (INPUT,  DATA,  OUTPUT) I 

<•  SOLVE  A  FIRST  ORDER  DIFFERENTIAL  EQUATION  •) 

{*  *’  -  SQR  <*)  *) 

(*  SOLUTION  IS  IN  INTERVAL  FORM  *) 

CONST  DIM  -  ISf 

TYPE  DIMTYPE  -  1..DIM| 

INTERVAL  -  RECORD  INF,  SUP  t  REAL  END) 

I VECTOR  -  ARRAY [DIMTYPE]  OF  INTERVAL} 

ITAYLOR  -  RECORD  LENGTH  s  DIMTYPE} 

T  i  REAL} 

TC  I  I VECTOR  END} 

CHOICE  “  1 . .3} 

VAR  FLAG  I  CHOICE} 

I,  INI  1  DIMTYPE} 

X,  Y,  YPRIMB  t  ITAYLOR} 

DATA  :  TEXT  1 

EPSILON, 

COMPOUND  I  REAL} 

(#  Transfer  Functions  *) 

FUNCTION  INTPT  (  RA  >  REAL  )  I  INTERVAL} 

EXTERNAL  41} 

FUNCTION  IMTVAL  (  RA, RB 1  REAL  )  1  INTERVAL} 

EXTERNAL  42} 

FUNCTION  IINF  (  Al  INTERVAL)  !  REAL} 

EXTERNAL  43} 

FUNCTION  I8UP  (  Al  INTERVAL)  1  REAL} 

EXTERNAL  44} 

(*  Comparison*  •) 

OPERATOR  <«  ( A,Bl  INTERVAL  )  RESt  BOOLEAN  } 

EXTERNAL  48} 

OPERATOR  >-  (A, Si  INTERVAL  )  RESi  BOOLEAN  I 
EXTERNAL  50} 

OPERATOR  IN  (RAiRXAL}  Bi  INTERVAL)  RESi  BOOLEAN} 
EXTERNAL  47} 

OPERATOR  IN  (XAi  INTEGER}  Bl  INTERVAL)  RESi  BOOLEAN} 
EXTERNAL  46} 

OPERATOR  X  (A.Bi  INTERVAL  )  RESi  BOOLEAN  } 

EXTERNAL  52} 

(*  Lsttlcs  Operators  *) 

OPERATOR  ♦+  (A, 81  INTERVAL)  RE81  INTERVAL} 

EXTERNAL  63} 

OPERATOR  ••  (A.Bi  INTERVAL)  RESt  INTERVAL} 

EXTERNAL  60} 
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(*  Arithmetic  Operators  *) 


OPERATOR  +  (  At  INTERVAL  )  RXSl  INTERVAL) 

EXTERNAL  67 j 

OPERATOR  -  (  As  INTERVAL  )  RXSl  INTERVAL; 

EXTERNAL  66; 

OPERATOR  +  (  A, Bl  INTERVAL  )  RESl  INTERVAL) 

EXTERNAL  68) 

OPERATOR  +  (  XAt  INTEGER)  Bl  INTERVAL  )  RES)  INTERVAL) 

EXTERNAL  69) 

OPERATOR  +  (  At  INTERVAL)  KB  l  INTEGER  )  RXSl  INTERVAL) 

EXTERNAL  70) 

OPERATOR  -  (  A , B 1  INTERVAL  )  RXSl  INTERVAL) 

EXTERNAL  73) 

OPERATOR  -  (  XAt  INTEGER)  Bt  INTERVAL  )  RX8t  INTERVAL) 

EXTERNAL  75) 

OPERATOR  -  (  At  INTERVAL)  KB)  INTEGER  )  RXSl  INTERVAL) 

EXTERNAL  74) 

OPERATOR  •  (  A , B l  INTERVAL  )  RXSl  INTERVAL) 

EXTERNAL  76) 

OPERATOR  *  (  XAt  INTEGER)  Bt  INTERVAL  )  RXSt  INTERVAL) 

EXTERNAL  79) 

OPERATOR  •  (  At  INTERVAL)  KB  I  INTEGER  )  RESl  INTERVAL) 

EXTERNAL  SO) 

OPERATOR  /  (  A,  Bt  INTERVAL  )  RXSl  INTERVAL) 

EXTERNAL  85) 

OPERATOR  /  (  XAt  INTEGER)  Bt  INTERVAL  )  RXSt  INTERVAL) 

EXTERNAL  S3) 

OPERATOR  /  (  At  INTERVAL)  KB t  INTEGER  >  RXSt  INTERVAL) 

EXTERNAL  86) 

FUNCTION  ISC  ALP  (VAR  A,  Bt  I  VECTOR)  AKDIM  t  INTEGER)  t  INTERVAL) 
EXTERNAL  88) 

(*  Standard  Functions  *) 

FUNCTION  IABS  (  'it  INTERVAL  )  s  REAL) 

EXTERNAL  101) 

FUNCTION  I8QR  (  Yl  INTERVAL  )  t  INTERVAL) 

EXTERNAL  102) 

FUNCTION  I8QRT  (  Yt  INTERVAL  )  t  INTERVAL) 

EXTERNAL  105) 

FUNCTION  IXXP  (  Yt  INTERVAL  )  I  INTERVAL) 

EXTERNAL  106) 

FUNCTION  ILK  (  Yt  INTERVAL  )  s  INTERVAL) 

EXTERNAL  107) 

FUNCTION  IARCTAN  (  Yt  INTERVAL  )  t  INTERVAL) 

EXTERNAL  108) 

FUNCTION  ISIN  (  Yt  IN1.RVAL  )  i  INTERVAL) 

EXTERNAL  109) 

FUNCTION  IC06  (  Yt  INTERVAL  )  t  INTERVAL) 

EXTERNAL  110) 
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(*  Input  /  Output  •) 


PROCEDURE  IRKAO  (  VAR  FiTKXT;  VAR  At  INTERVAL  ) I 
EXTERNAL  92  j 

PROCEDURE  WRITE  (  VAR  Pi  TOT ;  At  INTERVAL  )» 

EXTERNAL  91; 

FUNCTION  X8UM  (At  IVECTOR;  DIM i  DIKTYPI)  t  INTERVAL; 

VAR  Bt  I VECTOR; 

It  DIMTYPE; 

BEGIN 

FOR  I  i-  1  TO  DIM  DO  B[X]  i-  XNTPT  (1.0); 

I SUM  I-  X8CALP  (A,  B,  DIM) 

END;  (*  FUNCTION  ISDN  •) 

FUNCTION  IVRHULL  t  I VECTOR; 

VAR  It  DIMTYPE;  Ut  IVECTOR; 

BEGIN 

FOR  I  i-  1  TO  DIM  DO  0(1)  I-  XNTPT  (0.0); 

IVRHULL  t-  V 

END;  (•  FUNCTION  IVRHULL  *) 

FUNCTION  ITSQR  (Tt  ITATLOR )  t  ITATLOR;  (*  ITSQR(IT)  •) 

(*  Raquireat  IVRHULL,  I SCALP,  XBfiR  *) 

VAR  I,  J,  K,  HALF t  DIMTYPE; 

X,  Yt  IVECTOR; 

U  t  ITAYLOR; 

BEGIN 

X  i”  IVRHULL;  Y  i"  IVRHULL; 

U.LENGTR  t-  T.  LENGTH; 

U.T  t-  T.T; 

U.TC  t-  IVRHULL; 

U.TC(I)  I-  I8QR  (T.TC(I)); 

X(1]  t-  T.TC(1); 

FOR  X  t-  2  TO  U. LENGTH  DO 
BEGIN 

X(Xl  t-  T.TC(K); 

HALF  I-  X  DIV  2; 

FOR  0  t»  1  TO  HALF  DO 
BEGIN 

I  c-  X  -  J  ♦  1; 

Y(J)  t-  T.TCtU 
END;  (•  FOR  J  *) 

U.TC (X)  I"  2  •  I SCALP  (  X,  Y,  HALF); 

IF  X  MOO  2  “  1  THEN 
BEGIN 

HALF  t-  HALF  ♦  1; 

U.TC(X)  t"  U.TCtX)  ♦  I8QR  (T.TC(HALF)  ) 

END  (*  IF  •) 

END;  (•  FOR  X  *) 

ITSQR  t-  V 

END;  (•  FUNCTION  XT8QR  (ITATLOR)  •) 
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FUNCTION  MENU_CHOICE  I  CHOICE; 

VAR  I:  INTEGER; 

BEGIN 
WRITELN ; 

WRITELN  ('ENTER:  1  -  GIVE  NEW  INITIAL  CONDITIONS'); 

WRITELN  ('  2  -  CONTINUE  EXTENDING  THE  SOLUTION'); 

WRITELN  ('  3  -  STOP'); 

READ  (I); 

IF  ((I  >-  3)  OR  (I  <-  0))  THEN  I  :«  3; 

MENU_CHOICE  I 

END;  (*  FUNCTION  MENU  CHOICE  *) 


PROCEDURE  WRITE_INTERVAL  (INT:  INTERVAL); 

BEGIN 

WRITE  ('[',  INT. INF: 12,  ',  *,  INT.SUP:12,  ')'); 
END;  {•  PROCEDURE  WRITE  INTERVAL  *) 


PROCEDURE  PRNT_ITAY_COEF  (Y:  ITAYLOR;  INDEX:  DIMTYPE)  I 
BEGIN 

WRITE  ('!(',  INDEX: 5,  ')  -  '); 

WRITE_INTERVAL  (Y.TC  [INDEX]  ) ; 

WRITELN 

END;  (*  PROCEDURE  PRNT  ITAY_COEF  *) 


FUNCTION  INTERVAL  LENGTH  (INT:  INTERVAL)  :  REAL; 
BEGIN 

INTERVAL_LENGTH  :»  INT. SUP  -  INT. INF 
END;  (*  FUNCTION  INTERVAL  LENGTH  •) 


FUNCTION  RELATIVE ^LENGTH  (INT:  INTERVAL)  :  REAL; 
BEGIN 

RELATIVELENGTH  :»  2.0  •  (INT. SUP  -  INT. INF) 
/  (INT. SUP  +  INT. INF) 
END;  (•  FUNCTION  RELATIVE  LENGTH  ») 


FUNCTION  RELATIVEJSRROR  (INT:  INTERVAL)  :  REAL: 

BEGIN 

RELATIVE_ERROR  :-  2.0  *  INT. SUP  /  (INT. SUP  +  INT. INF) 
END;  (*  FUNCTION  RELATIVE  ERROR  •) 


BEGIN  (*  MAIN  PROGRAM  IDEQ_SOLVE  *) 

(* .  INITIALIZE  *) 

FLAG  :»  2; 

X.  LENGTH  :-  DIM; 

Y.  LENGTH  DIM; 

RESET  (DATA); 


WHXLI  rLAG  <-  2  DO  (*  LOOP  FOR  NEW  INITIAL  CONDITIONS  *) 

BEGIN 
FLAG  2| 

X. TC  s-  IVRNULLl 

Y. TC  i-  IVRNDLLl 

WRITELN  (‘READ  INTERVAL  INITIAL  CONDITIONS  XO,  Y(X0):')| 

IREAD  (DATA,  X.TC(1))|  IREAD  (DATA,  Y.TC[1])l 
WRITELM I  WRITELNl 

WRITE  ('INITIAL  CONDITIONS  AT  XO  -  '  )  I 


WRITE_INTERVAL  (X.TC(1])l  WRITELN  (',*»* 

WRITE  ( '  YO  »  '  )  I 

WRITE  INTERVAL  (Y.TC  Mill  WRITELN  ('♦*)» 

FIAG  m  2  DO  ( *  LOOP  FOR  ANALYTIC  CONTINUATION  *  ) 

BEGIN 

(• .  READ  STEP  SIZE  *) 


WRITELN  CENTER  STSPSIZE  X  -  XO:  ')l 
READ  (X.TC 
Y.T  X.Ti 

WRITELN  ( 'Computing  series  terns 

FOR  I  !-  2  TO  DIM  DO  (*  LOOP  FOR  SERIES  GENERATION  •) 

BEGIN 

(•  YOUR  FIRST  ORDER  DIFFERENTIAL  EQUATION  GOES  HEPEi  > 

YPRIME  s-  ITSQR  (Y); 

IMl  i-  I  -  1| 

Y.TC (I)  :•  YPRIME. TC  [IM1]  *  INTPT  (Y.T  /  IM1)| 

END)  ( *FOR* ) 

(• .  PRINT  TABLE  *) 

EPSILON  0.5  •  RELATIVE_LBNGTH  (Y.TC[1])l 
COMPOUND  «-  1.0; 

WRITELNl  WRITELNl 

WRITELN  ('Step  Left  Right  Computed 

Theoretical’)! 

WRITELN  ( •  Endpoint  Endpoint  Instability 

Instability' )i 

WRITELN i 

FOR  I  1  TO  DIM  DO  (•  LOOP  FOR  ERROR  MEASUREMENT  *) 

BEGIN 

COMPOUND  : -  COMPOUND  *  (1.0  *  EPSILON ) | 

WRITE  (I:3)|  WRITE  ('  * ) I 

WRITE_INTERVAL  ( Y  .TC  U  1  )  I 

WRITE  ('  ',  RELATIVE_ERROR  (Y.TC  [Ij >« 10); 

WRITE  ('  ',  COMPOUND:  10)1  WRITELN 

END!  ( *FOR* ) 

(* .  PERFORM  THE  ANALYTIC  CONTINUATION  *) 

Y.TC (1]  I SUM  (Y.TC,  DIM)| 

FOR  I  :«2  TO  DIM  DO  Y.TC(I)  INTPT  (0.0)| 

X.TC  (1 )  X.TC  1 1 1  +  INTPT  (X.T): 

WRITELNl 

WRITE  ('THE  VALUE  AT  X  -  *)|  WRITE_INTERVAL  (X.TC 1 1 1 ) I  WRITELNl 
WRITE  ('  IS  Y  -  *>l  WRITE_INTERVAL  (Y.TC(1)  )l  WRITELN  (' 

FLAG  MENU_CHOICE 
END  ( *WHILE*) 

END  (*WHILE* ) 

END.  (*MAIN  PROGRAM  IDEQ_SOLVE* ) 
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APPENDIX  C 


Source  Cod*  for  TAYLOR  and  I TAYLOR  Operators,  Functions,  and  Procedures 


Ths  Pascal-SC  source  cod*  for  the  operators,  function*,  and  procedure* 
described  In  this  report  is  contained  in  seven  files t 

1 .  Real  and  Interval  Taylor  addition  and  subtraction  operators 
( RIT_ADD. LIB ) ; 

2.  Reel  and  interval  Taylor  aultipllcatlon  operators  including  T9QR 
and  ITSQR  (RIT_HUL.LIB)» 

3.  Real  and  interval  Taylor  division  operators  (RIT_DIV.LIB)i 

4.  Real  Taylor  power  operators  and  funetlons  (RT_POW.LIB) » 

5.  interval  Taylor  power  operators  and  functions  (ITJPOW.LIB)j 

6.  Real  and  interval  Taylor  functions  (RIT_FNS.LIB) » 

7.  Utility  functions  and  procedures  (UTIL. LIB). 

Bach  file  is  headed  by  a  table  of  contents  to  assist  in  the  location  of  the 
needed  routines.  In  addition,  a  complete  table  of  contents  for  all  files  of 
TAYLOR  and  ITAYLOR  operators,  functions,  and  procedures  is  given  in  the  file 
contents. LIB.  The  first  line  of  cod*  for  each  subroutine  contains  a  short 
identification  in  coaMents  brackets,  for  exa*ple,  {•  K  +  T  •),  to  assist  in 
locating  it  with  the  help  of  a  text  editor. 
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C.  1.  a t«X  and  Inf  rv«l  Taylor  addition  and  Subtraction 


K\-z  w  i'-u  n 


(•  UT  ADD. LIB  -  RIAL  AND  I  NTS  AVAL  TAYLOR  ADO  AM)  ROBTRACT  «<<« 


Oentwtii 

♦  T  OPERATOR  4  (Ti  TAYLOR)  RM  t  TAYLOR/ 

X  4  T  OPERATOR  4  (Ki  INTEGER)  Ti  TAYLOR)  RM  I  TAYLOR) 

T  4  X  OPERATOR  4  (Ti  TAYLOR)  Ki  INTEGER)  RM  t  TAYLOR) 

X  4  T  OPERATOR  4  (Ri  RIAL)  Ti  TAYLOR)  RM  i  TAYLOR) 

T  4  R  OPERATOR  4  (Ti  TAYLOR)  At  RIAL)  RM  t  TAYLOR) 

T  4  T  OPERATOR  4  (TA,  TBi  TAYLOR)  RM  I  TAYLOR) 

4  XT  OPERATOR  4  (Tl  XTAYLOR)  RM  t  1TAYLOR) 

X  4  it  OPERATOR  4  (Xi  INTEGER  I  Tl  XTAYLOR)  RM  I  XTAYLOR) 

XT  4  X  OPERATOR  4  (Tl  XTAYLOR)  Kl  IDT  BOER)  RM  I  XTAYLOR) 

I  4  IT  OPERATOR  4  (X)  INTERVAL)  Tl  XTAYLOR)  RM  t  XTAYLOR) 

IT  4  I  OPERATOR  4  (Tt  XTAYLOR)  Xl  INTERVAL)  RM  I  XTAYLOR) 

IT  4  IT  OPERATOR  4  (TA,  TBs  ITAYLOR)  RM  I  XTAYLOR) 

-  T  OPERATOR  -  (Tl  TAYLOR)  RM  i  TAYLOR) 

X  -  T  OPERATOR  -  (El  INTEGER)  Tl  TAYLOR)  RM  I  TAYLOR) 

T  -  X  OPERATOR  -  (Tl  TAYLOR)  Xl  INTEGER)  RM  I  TAYLOR) 

R  -  T  OPERATOR  -  (Rl  RIAL)  Tt  TAYLOR)  RM  I  TAYLOR) 

T  -  R  OPERATOR  -  (Tl  TAYLOR)  Rl  RIAL)  RM  I  TAYLOR) 

T  -  T  OPERATOR  -  (TA,  TBi  TAYLOR)  RM  I  TAYLOR) 

-  XT  OPERATOR  -  (Tl  XTAYLOR)  RM  I  XTAYLOR) 

E  -  XT  OPERATOR  -  (X)  INTEGER)  Tl  XTAYLOR)  RM  I  ITAYLOR) 

IT  -  K  OPERATOR  -  (T)  XTAYLOR)  Xl  INTEGER)  DM  I  XTAYLOR) 

X  -  IT  OPERATOR  -  (El  INTERVAL)  Tl  XTAYLOR)  RM  I  XTAYLOR) 

IT  -  I  OPERATOR  -  (Tl  XTAYLOR)  El  INTERVAL)  RM  i  ITAYLOR) 

XT  -  IT  OPERATOR  -  (TA,  TRi  ITAYLOR)  RM  I  XTAYLOR) 


*) 


OPERATOR  4  (Tl  TAYLOR)  RM  i  TAYLOR) 

BEGIN  RM  I-  T  END)  (•  4  TAYLOR  •) 


(•  4  T  •) 


OPERATOR  4  (Ki  INTEGER)  Tl  TAYLOR)  RH  I  TAYLOR)  (•  X  4  T  •) 

VAR  Ol  TAYLOR) 

BEGIN 
0  l«  T) 

O.TCIU  «•  K  4  T.TC(U) 

RH  i-  0 

END)  (•  INTEGER  4  TAYLOR  •) 


OPERATOR  4  (Tl  TAYLOR)  Xl  INTEGER)  RH  t  TAYLOR)  (*  T  4  X  •> 

VAR  Vi  TAYLOR) 

BEGIN 
0  I"  T) 

O.TCIU  i“  T.TCtl]  ♦  X) 

RH  i-  0 

END)  (•  TAYLOR  4  INTEGER  •) 
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OPERAYOR  ♦  (Rt  REAL)  Ti  TAYLOR)  RES  I  TAYLOR |  (*  R  +  T  •) 

VAR  Ut  TAYLOR) 

BEGIN 
0  l-  T| 

O.TCH1  i-  R  +  U.TC (1)» 

RES  i-  0 

END)  <*  RIAL  +  TAYLOR  *) 


OPERATOR  +  (Tt  TAYLOR )  Rt  RIAL)  RES  t  TAYLOR)  (•  T  +  R  *) 

VAR  Ot  TAYLOR) 

BEGIN 
0  I-  T) 

O.TCJ1]  i-  U.TC111  ♦  R) 
i  RES  t-  V 

END)  (*  TAYLOR  ♦  REAL  •) 


OPERATOR  +  (TA,  TBt  TAYLOR)  RES  t  TAYLOR)  (*  T  +  T  *) 

VAR  Ui  TAYLOR)  It  DIMTYPE) 

BEGIN 

IP  TA.T  <>  TB.T 
THEN  BEGIN 

WRITE  (  'ERROR:  ADDITION  OP  TAYLOR  VARIABLES'  )| 

WRITE LN  ( '  WITH  UNEQUAL  SCALE  FACTORS' )) 

SVR  (0)  END  (•  RETURN  TO  OPERATING  SYSTEM*) 

ELSE  BEGIN 

U.LENGTH  I-  T A. LENGTH) 

IP  U.LENGTH  >  TB. LENGTH  THEN  U.LENGTH  t-  TB. LENGTH) 

U.T  t-  TA.T) 

FOR  I  I-  1  TO  U.LENGTH  DO  U.TCtI)  t-  TA.TCfll  ♦  TB.TCCII) 

RES  l-  U 

END  (•  ELSE  *) 

END)  (*  TAYLOR  ♦  TAYLOR  *) 


OPERATOR  <Tt  ITAYLOR)  RES  t  ITAYLOR)  <*  ♦  IT  *) 

BEGIN  RES  t-  T  END)  (*  ♦  ITAYLOR  *) 


OPERATOR  *  (Kt  INTEGER)  Tt  ITAYLOR)  RES  t  ITAYLOR)  (*  K  *  IT  *) 

VAR  Ut  ITAYLOR) 

BEGIN 
U  t-  T) 

U.TC[1]  t-  K  ♦  T.TCll]) 

RES  l-  U 

END)  (*  INTEGER  *  ITAYLOR  •) 


OPERATOR  *  (Tt  ITAYLOR)  Kt  INTEGER)  RES  t  ITAYLOR)  (*  IT  *  K  *) 

VAR  Ut  ITAYLOR) 

BEGIN 
U  t“  T) 

U.TCP1  t-  T.TC(1)  ♦  X) 

RES  t*  U 

END)  (*  ITAYLOR  +  INTEGER  *) 
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(•  I  4  IT  *) 


OPERATOR  4  (Xl  INTERVAL;  Tt  ITAYLOR)  Mt  I  ITAYLOR; 
VXR  Ot  ITAYLOR; 

begin 

0  I"  T; 

O.TCd)  t-  X  4  T.TCdll 

RX8  i-  0 

END;  («  INTERVAL  +  ITAYLOR  *) 


OPERATOR  4  (Ti  ITAYLOR;  Kt  INTERVAL)  RI8  t  ITAYLOR;  (•  IT  4  I  •) 

VAR  Oi  ITAYLOR; 

BEGIN 
U  »■  T; 

U.TCd)  i-  T.TCdl  4  Kf 
RES  t-  0 

BHD;  (•  ITAYLOR  4  INTERVAL  *) 


OPERATOR  4  (TA,  TBi  ITAYLOR)  RE8  t  rTAYLOR;  (•  IT  4  IT  *) 

VAR  Ot  ITAYLOR;  It  DIMTYPS; 

BEGIN 

IF  TA.T  <>  TB.T 

THEN  BEGIN 

WRITE  (‘ERROR!  ADDITION  OF  ITAYLOR  VARIABLES'); 

WRITXXJI  (‘  WITS  UNEQUAL  SCALE  FACTORS'); 

8VR  (0)  END  (•  RETURN  TO  OPERATING  8Y8TEM* ) 

ELSE  BEGIN 

0. LENGTH  I-  TA. LENGTH; 

IF  U. LENGTH  >  TB. LENGTH  THEN  U. LENGTH  »•  TB. LENGTH; 

U.T  t-  TA.T; 

FOR  I  l«  1  W  U. LENGTH  DO  O.TCU)  I-  TA.TClI)  4  TB.TClI] ; 

RES  i-  U 

END  (•  ELSE  •) 

END;  <*  ITAYLOR  4  I TAYLOR  *) 


OPERATOR  -  (Tl  TAYLOR)  RES  I  TAYLOR;  (*  -  T  •) 

VAR  Ui  TAYLOR;  It  DZHTYPE; 

BEGIN 

U. LENGTH  i-  T. LENGTH; 

U.T  i-  T.T; 

FOR  I  i»  1  TO  o. LENGTH  DO  O.TC[I]  t-  -  T.TC(I) ; 

RES  i-  U 

END;  (•  -  TAYLOR  *) 


OPERATOR  -  (Kt  INTEGER;  Tt  TAYLOR)  RES  t  TAYLOR;  (•  X  -  T  *) 

VAR  OtTAYLOR;  It  DIMTYPE; 

BEGIN 

U. LENGTH  t-  T. LENGTH; 

U.T  t"  T.T; 

U.TCI1]  t-  X  -  T.TCdl » 

FOR  It-  2  TO  U. LENGTH  DO  O.TCtI]  t-  -  T.TCUJ ; 

RES  t-  U 

END;  (•  INTEGER  -  TAYLOR  •) 
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(*  I  -  X  *> 


* 


OPERATOR  -  (Ti  TAYLOR;  Ki  INTEGER)  AES  t  TRYLOR; 
TAR  Ut  TAYLOR; 

SWIM 
0  1“  T; 

O.TC[1]  i-  T.TCll]  -  X; 

res  i-  o 

xm>;  t*  TAYLOR  -  INTEGER  *) 


OPERATOR  -  (Ri  REAL;  Ti  TAYLOR)  RES  t  TAYLOR;  (*  R  -  T  •) 

VAR  Ut  TAYLOR;  It  DINTYPE; 

BEGIN 

0. LENGTH  t-  T. LENGTH; 

|  U.T  i-  T.T; 

O.TCID  t-  R  -  T.TCll]; 

TOR  1  1*210  0. LENGTH  DO  O.TC[I]  I-  -  T.TC[I] I 
RES  i-  0 

END;  (*  REAL  -  TAYLOR  •) 


OPERATOR  -  (Tl  TAYLOR;  Rt  REAL)  RE8  t  TAYLOR;  (*  T  -  R  *) 

VAR  Ut  TAYLOR; 

BEGIN 
O  <*■  T; 

U.TCM]  1-T.TCH]  -  R; 

RES  t-  O 

END;  (•  TAYLOR  -  REAL  •) 


1 


OPERATOR  -  (TA.  TBt  TAYLOR)  RES  t  TAYLOR;  (*  T  -  T  *) 

VAR  0i  TAYLOR;  It  DDfTYPE; 

BEGIN 

IP  TA.T  <>  TB.T 
THEM  BEGIN 

WRITE  ( ‘ERROR!  SUBTRACTION  OP  TAYLOR  VARIABLES*)  I 
WRITELN  ( '  WITH  UNEQUAL  SCALE  FACTORS' )l 

SVR  (0)  END  (•  RETURN  TO  OPERATING  SYSTEM*) 

ELSE  BEGIN 

U. LENGTH  i«  TA. LENGTH; 

IP  U. LENGTH  >  TB. LENGTH  THEM  U. LENGTH  I-  TB. LENGTH; 

O.T  t"  TA.T; 

POR  I  t-  1  TO  U. LENGTH  DO  U.TCtI)  t-  TA . TC [ I ]  -  TB.TC[I]; 

RES  t-  U 
END  (*  ELSE  •) 

END;  (•  TAYLOR  -  TAYLOR  *) 


OPERATOR  -  (Tt  ITAYLOR)  RX8  I  ITAYLOR;  (*  -  IT  *) 

VAR  Ut  ITAYLOR;  It  DINTYPE; 

BEGIN 

U. LENGTH  t*  T. LENGTH; 

O.T  t“  T.T; 

POR  I  t-  1  TO  U. LENGTH  DO  U.TCfl)  t-  -  T.TCll]; 

RES  Is  U 

END;  (•  -  ITAYLOR  •) 
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{•  E  -  IT  *) 


OPERATOR  -  (Xt  INTEGER)  Tl  itatlor)  res  I  ITAYLOR) 

VAR  Ot ITAYLOR)  Zt  DWTYPtl 

inn 

O.LENQTH  i-  T. LENGTH) 

O.T  T.T» 

U.TCHJ  l-  R  -  T.TC( t]  I 

POR  It-  2  TO  0. LENGTH  DO  O.TCtI)  I-  -  T.TCtl] I 

RRfl  «-  0 

KMDl  (•  INTEGER  -  ITAYLOR  •) 

OPERATOR  -  (Tl  ITAYLORl  Ki  INTEGER)  REE  I  ITAYLORl  (•  IT  -  K  •) 

VAR  Ol  I TAYLOR) 

REGIN 

0  f-  T I 

O.TC[1]  i-  T.TCtl]  -  El 

RES  i-  0 

END)  (*  ITAYLOR  -  INTEGER  *) 

OPERATOR  -  (Kt  INTERVAL I  Tt  ITAYLOR)  RES  l  ITAYLOR)  (*  I  -  IT  *) 

VAR  0 1  ITAYLOR |  It  OURTPEi 
BEGIN 

U. LENGTH  i-  T. LENGTH) 

O.T  i-  T.Tl 

O.TC(t]  l"  E  -  T.TClHt 

POR  Ii*  2  TO  0. LENGTH  DO  O.TC(I]  l«  -  T.TCtl]) 

RES  i-  0 

END)  (•  INTERVAL  -  ITAYLOR  *> 

OPERATOR  -  (Tt  ITAYLOR)  Et  INTERVAL)  RES  I  ITAYLOR)  (•  IT  -  I  •) 

VAR  Ot  ITAYLOR » 

BEGIN 
0  l»  T) 

O.TC[1)  i-  T.TCtl]  -  El 

RES  t-  0 

END)  (*  ITAYLOR  -  INTERVAL  •) 

OPERATOR  -  (TA.  TEt  ITAYLOR)  RES  t  ITAYLOR)  t«  IT  -  IT  •) 

VAR  Ot  ITAYLORt  It  DINTYPE) 

BEGIN 

IP  TA.T  <>  TB.T 
THEN  BEGIN 

WRITE  ( 'ERROR)  SUBTRACTION  OP  ITAYLOR  VARIABLES' ) I 
NRIT1LM  ('  WITH  UNEQUAL  SCALE  PACTORS')) 

SVR  (0)  end  (*  RETURN  TO  OPERATING  SYSTEM* ) 

ELSE  BEGIN 

0. LENGTH  t-  TA. LENGTH) 

IP  0. LENGTH  >  TB. LENGTH  THEM  0. LENGTH  I-  TB. LENGTH) 

O.T  1“  TA.T) 

POR  I  t-  1  TO  0. LENGTH  DO  O.TCtI)  t-  TA.TCtI]  -  YB.TCtI] l 
RES  i-  O 
END  (*  ELSE  •) 

END)  (*  ITAYLOR  -  ITAYLOR  •) 


(•  END  OP  RIT  ADD. LIB 


>»»>»»»»»»»»»»»»»»»»»»  •) 


C.2.  Rgfl  and  Interval  Taylor  Multiplication  Op«r«tori.  Including  TSOR 
and  IT80R. 


\ 

i 


I 

! 


(*  RIT  MOL. LIB  -  RIAL  AMD  INTERVAL  TAYLOR  MULTIPLY  <<<<<<<<<<<<< 


Contents: 


X  *  T 
T  *  K 
R  *  T 
T  *  R 
T  *  T 
TSQR(T) 

X  •  IT 
IT  *  K 
I  •  IT 
IT  •  I 
IT  •  IT 
ITSQR(  IT) 


OPERATOR  *  (Xi  INTEGER!  Tt  TAYLOR)  RES  s  TAYLOR I 
OPERATOR  *  (Ti  TAYLOR;  Kt  INTEGER)  RES  t  TAYLOR; 
OPERATOR  *  (Ri  REAL;  Ti  TAYLOR)  RX8  I  TAYLOR; 
OPERATOR  *  (T>  TAYLOR;  R;  REAL)  RES  t  TAYLOR; 
OPERATOR  *  (TA,  TBi  TAYLOR)  RES  t  TAYLOR; 

FUNCTION  TSQR  (Ti  TAYLOR)  t  TAYLOR; 

OPERATOR  *  (X;  INTEGER;  Tt  ITAYLOR)  RES  t  ITAYLOR; 
OPERATOR  *  (Tt  ITAYLOR;  Xt  INTEGER)  RES  s  ITAYLOR; 
OPERATOR  •  (Xt  INTERVAL;  Tt  ITAYLOR)  RES  t  ITAYLOR; 
OPERATOR  •  (Tt  ITAYLOR;  Kt  INTERVAL)  RES  t  ITAYLOR; 
OPERATOR  *  (TA,  TBi  ITAYLOR)  RES  t  ITAYLOR; 

FUNCTION  ITSQR  (Tt  ITAYLOR)  t  ITAYLOR; 


OPERATOR  *  (Kt  INTEGER;  Tt  TAYLOR)  RE8  t  TAYLOR;  (*  X  •  T  *) 

VAR  ItDIHTYPE;  Ut  TAYLOR; 

BEGIN 

U. LENGTH  i-  T. LENGTH; 

U.T  t-  T.T; 

FOR  I  t-  1  TO  U. LENGTH  DO  U.TCtI)  t-  X  *  T.TCU]  ; 

RES  :*  U 

END;  (»  INTEGER  •  TAYLOR  •) 


OPERATOR  •  (Tt  TAYLOR;  Kt  INTEGER)  RES  t  TAYLOR;  (•  T  •  X  •) 

VAR  ItDIHTYPE;  Ut  TAYLOR; 

BEGIN 

U. LENGTH  t-  T. LENGTH; 

U.T  t»  T.T; 

FOR  I  1-  1  TO  U. LENGTH  DO  U.TC(l)  I-  T.TCtI)  *  X; 

RES  I*  U 

END;  (•  TAYLOR  *  INTEGER  *) 


OPERATOR  *  (Rt  REAL;  Tt  TAYLOR)  RES  t  TAYLOR;  (*  R  *  T  •) 

VAR  It  DIMTYPE;  Ut  TAYLOR; 

BEGIN 

U. LENGTH  i-  T. LENGTH; 

U.T  t“  T.T; 

FOR  I  t-  1  TO  U. LENGTH  DO  U.TC(I]  t-  R  •  T.TC(I); 

RES  t-  0 

END;  (*  REAL  •  TAYLOR  •) 
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I 


OPERATOR  •  (Tl  TAYLOR)  Ri  REAL)  RES  ■  TAYLOR) 
VAR  I l  DIRT TP* |  Oi  TAYLOR) 


(*  T  *  R  *) 


RHGIH 

0. LENGTH  I-  T. LENGTH) 

U.T  «-  T.T) 

FOR  I  1  TO  U. LENGTH  DO  D.TCtX]  I-  T.TC(I)  *  R) 
RES  t-  0 

EDO)  <*  TAYLOR  *  REAL  *) 


OPERATOR  *  (TA,  TBi  TAYLOR)  RES  «  TAYLOR)  (*  T  *  T  •) 

(*  REQUIRES)  VRNULL,  SCALP  •> 

VAR  I,  J,  K  t  DIMTYPE) 

X,  Y  t  RVECTOR) 

U  )  TAYLOR) 

BEGIN 

IP  TA.T  <>  TB.T 
THEN  BEGIN 

WRITE  (’ERROR)  MULTIPLICATION  OP  TAYLOR  VARIABLES')) 

WRITE LM  ( '  WITH  UNEQUAL  SCALE  FACTORS ' ) ) 

SVR  (0)  END  (*  RETURN  TO  OPERATING  SYSTEM  *) 

ELSE  BEGIN 

U. LENGTH  I-  TA. LENGTH) 

IP  U. LENGTH  >  TB. LENGTH  THEM  U. LENGTH  )■  TB. LENGTH) 

U.T  t-  TA.T) 

X  )■  VRNULL)  Y  I-  VRNULL) 

POR  I  i-  1  TO  U. LENGTH  DO  BEGIN 
XIIJ  )»  TA.TC(Il) 

POR  J  I"  1  TO  I  DO  BEGIN 
*  l»  I  -  J  ♦  t| 

YfJJ)-  TB. TC [X] 

END)  <•  POR  J  LOOP  •) 

U.TCtI]  )-  SCALP  (X,  Y,  0) 

END  <•  FOR  I  LOOP  *) 

END)  (*  ELSE  *) 


RE8  )-  U 

END)  (•  TAYLOR  «  TAYLOR 


•) 


-49- 


FUNCTION  TSQR  (Ti  TAYLOR)  i  TAYLOR; 

<«  taquirMi  VRMULL,  SCALP,  SQR  *) 

VAR  I,  J,  K,  HALT;  DIMTYPC; 

X,  Yt  RVECTOR; 

0  i  TAYLOR; 

BEGIN 

X  l-  VRMULL;  Y  l”  VRMULL; 

U.LENGTH  i-  T. LENGTH; 

U.T  t-  T.T; 

U.TC  l-  VRMULL; 

U.TC(1]  «•  SQR  (T.TCtI] ); 

X(11  «-  T.TC[1]; 

FOR  K  I-  2  TO  U.LENGTH  DO 
BEGIN 

X[K]  «-  T.TC  [X]  I 
HALF  l-  K  DIV  2; 

FOR  J  1-  1  TO  HALF  DO 
BEGIN 

I  »-  K  -  J  ♦  1» 

Y[J]  «-  T.TCtI] 

END;  (•  FOR  J  •) 

U.TCtX]  «-  2.0  •  SCALP  (  X,  Y,  0); 

IF  X  HOD  2  -  1  THEN 
BEGIN 

HALF  I-  HALF  +  1; 

U.TCtX]  «-  U.TCtX]  ♦  SQR  (T.TCtHALF]  ) 
END  (•  IF  *) 

END;  (•  FOR  X  •) 

TSQR  :«  U 

END;  (•  FUNCTION  TSQR  (TAYLOR)  •) 


OPERATOR  •  (X;  INTEGER;  T:  ITAYLOR)  RES  t  ITAYLOR; 

VAR  I iDIMTYPE;  Ui  ITAYLOR; 

BEGIN 

U. LENGTH  j-  T. LENGTH; 

U.T  !■  T.T; 

FOR  I  I-  1  TO  U. LENGTH  DO  U.TCtU  »-  X  •  T.TCtI); 
RES  i-  U 

END;  (*  INTEGER  *  ITAYLOR  •) 


OPERATOR  •  (Ti  ITAYLOR;  Xi  INTEGER)  RES  <  ITAYLOR; 

VAR  IiDIMTYPE;  Ui  ITAYLOR; 

BEGIN 

U. LENGTH  1-  T. LENGTH; 

U.T  T.T; 

FOR  I  1  TO  U. LENGTH  DO  U.TClI]  »-  T.TCtI)  *  X; 
RES  «-  U 

END;  (*  ITAYLOR  *  INTEGER  •) 
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(*  Z  *  IT  •) 


OPERATOR  *  (Kt  INTERVAL;  Tl  ITATLOR)  RES  I  ITATLOR; 
VAR  IiDIWTYPE;  Ut  ITATLOR; 

BEGIN 

□.LENGTH  1“  T. LENGTH; 

D.T  J“  T.T; 

FOR  Z  t-  1  TO  U. LENGTH  DO  O.TC(I]  t-  K  •  T.TCUll 
RES  l-  O 

END;  <•  INTERVAL  •  ITATLOR  •) 


OPERATOR  *  (Tt  ITATLOR;  Kt  INTERVAL)  RES  t  ITATLOR;  (•  IT  *  I  *) 

VAR  I s  DIMTYPS ;  Ot  ITATLOR; 

BEGIN 

U. LENGTH  I-  T. LENGTH; 

U.T  t-  T.T; 

FOR  I  t-  1  TO  0. LENGTH  DO  D.TCII]  «“  T.TC(I)  *  K; 

RES  a 

END;  (*  ITATLOR  *  INTERVAL  *) 


OPERATOR  *  (TA,  TBl  ITATLOR)  RES  t  ITATLOR;  (•  IT  *  IT  *) 

(*  Rtqninii  IVRNULL,  ZSCALP  *) 

VAR  I,  X  i  DINTTPB; 

X,  T  i  I VECTOR; 

O  >  ITATLOR; 

BEGIN 

IF  TA.T  <>  TB.T 
THEN  BEGIN 

WRITE  ( 'ERROR;  MULTIPLICATION  OP  ITATLOR  VARIABLES’)! 

WRITE!*  (•  WITH  UNEQDAL  SCALE  FACTORS’); 

SVR  (0)  END  (•  RETURN  TO  OPERATING  8T8TEN  •> 

ELSE  BEGIN 

X  t-  IVRNULL;  T  »•  IVRNULL; 

U. LENGTH  i-  TA. LENGTH; 

IF  U. LENGTH  >  TB. LENGTH  THEN  U. LENGTH  t-  TB. LENGTH; 

U.T  «■  TA.T; 

FOR  I  l-  1  TO  U. LENGTH  DO  BEGIN 
XtU  l-  TA.TCtll; 

FOR  J  lr  1  TO  I  DO  BEGIN 
K  t«  I  -  J  ♦  1; 

TtJll-  TB.TC[KJ 
END;  (•  FOR  J  LOOP  •) 

U.TCIIJ  l-  I SCALP  (X,  T,  I) 

END  (*  FOR  I  LOOP  *) 

END;  (•  ELSE  *) 


RES  t-  U 

END;  (•  ITATLOR  •  ITATLOR  •> 
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(•  ITSQR(IT)  •) 


FUNCTION  ITSQR  (Ti  ITAYLOR)  t  ITAYLOR) 

(•  Raquir*S!  IVRHOLL,  1 8CALP,  ISQR  *) 

VAR  I,  J,  K,  HALF)  DIMTYPE) 

X,  »l  I  VECTOR  J 
U  I  ITAYLOR I 

BBC  It) 

X  IVWTOLLl  Y  «-  IVRNULL) 

U. LENGTH  I-  T.  LENGTH) 

O.T  «“  T.T) 

U.TC  I-  IVRNOLLf 

O.TCtl]  t-  ISQR  (T.TC(ll)) 

X[1]  J-  T.TC[1]) 

FOR  K  «-  2  TO  0.  LENGTH  00 
BEGIN 

X[K]  i«T.TClK]| 

HALF  :«  X  DIV  2) 

FOR  J  i-  1  TO  HALF  00 
BEGIN 

I  K  -  J  ♦  1| 

YIJJ  i«  T.TCII] 

END)  (•  FOR  J  *) 

O.TCfX)  »-  2  •  I  SCALP  (  X,  Y,  HALF)) 

IF  K  HOD  2  -  1  THEN 
BEGIN 

HALF  I-  HALF  ♦  1) 

U.TC(X))-  U.TC(K)  ♦  ISQR  (T.TClHALF)  ) 
END  (*  IF  •) 

END)  (•  FOR  X  *) 

ITSQR  i-  U 

END)  (*  FUNCTION  ITSQR  (ITAYLOR)  •) 


(*  END  OF  RIT  NUL.LIB 


>>>»»»»»»»»»»»»>»»»»»»»  •) 
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C.3.  B»«l  and  Interval  T«vlor  Division  Op«r«tor!. 


(•  HIT  DIV.LIB  -  MU  AMO  INTERVAL  TAYLOR  DIVIDE  <<<<<<<<<<<<<<< 


Content* i 

T  /  K 
K  /  T 
R  /  T 
T  /  R 
T  /  T 
K  /  IT 
IT  /  K 
I  /  IT 
IT  /  I 
IT  /  IT 


OPERATOR  / 
OPERATOR  / 
OPERATOR  / 
OPERATOR  / 
OPERATOR  / 
OPERATOR  / 
OPERATOR  / 
OPERATOR  / 
OPERATOR  / 
OPERATOR  / 


{El  INTEGER)  Tt  TAYLOR)  RES  t  TAILOR) 

(Tt  TAYLOR)  Kt  INTEGER)  RES  I  TAYLOR) 

(Si  REAL)  Ti  TAYLOR)  RES  i  TAYLOR; 

(Ti  TAYLOR)  Ri  REAL)  RES  I  TAYLOR) 

(TA,  TO l  TAYLOR)  RES  i  TAYLOR) 

(El  INTEGER)  Tt  ITAYLOR)  RES  I  ITAYLOR) 
(Ti  ITAYLOR)  El  INTEGER)  RES  t  ITAYLOR) 
(El  INTERVAL)  Tt  ITAYLOR)  RES  I  ITAYLOR) 
(Tl  ITAYLOR)  Et  INTERVAL)  RES  1  ITAYLOR) 
(TA,  TBi  ITAYLOR)  RE8  I  ITAYLOR) 


OPERATOR  /  (Ti  TAYLOR)  El  INTEGER)  RES  I  TAYLOR)  (•  T  /  E  *) 

VAR  Oi  TAYLOR)  II  DIHTYPE) 

BEGIN 

O. LENGTH  I-  T. LENGTH) 

O.T  I-  T.T) 

IT  E  -  0 

THEN  BEGIN 

miTELN  ('ERROR!  DIVISION  OP  TAYLOR  VARIABLE  BY  ZERO')) 

SVR  (0)  END  (•  RETURN  TO  OPERATING  SYSTEM  *) 

ELSE  BEGIN 

FOR  I  I-  t  TO  U. LENGTH  DO  O.TCtI)  I-  T.TC(I)  /  E 
END)  (•  ELSE  •) 

RES  I*  U 

END)  (•  TAYLOR  /  INTEGER  •) 
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(*  X  /  T  •) 


I 


OPERATOR  /  (Xt  INTEGER)  Ti  TAYLOR)  RES  I  TAYLOR ) 
<•  Mqolraai  VRNOLL  •) 

VAR  0,  Nt  TAYLOR) 

I,  J,  IH,  RXVi  DINTYP1) 

RATIO!  REAL) 

BIGIN 

0. LENGTH  t-  T. LENGTH) 

U.Ti-  T.T) 


U.TC  I-  VRNOLL)  N.TC  I-  VRNOLL)  (*  ZERO  RESULT  AND  NONE  VECTORS  *) 


1 

1 


IP  T.TCMJ  <>  0.0 

THEN  BEGIN  {*  NORMAL  DIVIDE  •) 

RATIO  i-  1.0  /  T.TCl 1 ] ) 

O.TCI1]  »•  X  *  RATIO) 

POR  I  I-  2  TO  0. LENGTH  DO 
BEGIN 

POR  J  «-  1  TO  (1-1)  DO 

BEGIN 

REV  i-  I  -  J  ♦  1) 

H.TC[J]  I-  T.TCfREVJ 
END)  (•  POR  J  *) 

U.TC [I]  1“  -  SCALP  (U.TC,  N.TC,  0)  *  RATIO 
END  (*  FOR  I  •) 

HMD  (*  THEN  •) 

ELSE  IP  X  <>  0 
THEN  BEGIN 

WRITE  ('ERROR!  DIVISION  BY  TAYLOR  VARIABLE’)) 

HRITELN  (  '  EQUAL  TO  ZERO')) 

SVR  (0)  (•  return  to  operating  system  •) 

END  (•  THEN  •) 

ELSE  BEGIN  (•  APPLY  L' HOSPITAL'S  ROLE  *) 

IP  U. LENGTH  -  1 
THEN  BEGIN 

WRITE  ('ERROR!  DIVISION  WITH  TAYLOR  VARIABLE' ) ) 

WRITELN  ('  OP  LENGTH  ONE.')) 

SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  •) 

END)  (*  IP  •) 


U. LENGTH  l-  T. LENGTH  -  1) 

POR  I  1-1  TO  0. LENGTH  DO 
BEGIN 

IP1  i-  I  ♦  1) 

D.TC(I)  !-  T.TC[IP1] 

END)  (•  POR  I  *) 

0  I-  0  /  0  (•  THIS  18  A  RECURSIVE  CALL  •) 

END)  (*  IP  •) 

RES  l-  U 

END)  (•  INTEGER  /  TAYLOR  *) 
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(*  R/T  •> 


OPERATOR  /  (R»  REAL;  Ti  TAYLOR)  REE  l  TAYLOR) 

(*  Raqulrcai  vrmull  •) 

VAR  0,  «•  TAYLOR) 

X,  J,  XVI,  REV i  DXMmSl 
RATIO)  REAL) 

BROIM 

U. LENGTH  i-  T. LENGTH) 

U.Ti-  T.T) 


O.TC  t-  VRMOLL)  R.TC  »-  VRMOLL)  {*  EBRO  REVOLT  AMD  MORE  VECTORS  *) 
IV  T.TCH1  <>  0.0 

THEM  REGIE  (*  NORMAL  DIVIDE  *) 

RATIO  i-  1.0  /  T.TC[1J) 

O.TCtll  l*  R  •  RATIO) 

VOR  I  I-  a  TO  0. LENGTH  DO 
BEGIN 

FOR  J  )■  1  TO  (X-t)  DO 

BEGIN 

REV  )•  I  -  J  ♦  t» 

W.TC[<J]  i«  T.TCtREV] 

END)  (*  FOR  J  »> 

U.TCtX]  )•  -  ECALF  (O.TC,  N.TC,  0)  *  RATIO 
END  (*  FOR  I  *) 

END  (*  THEN  *) 

ELSE  IF  R  <>  0 
THEN  BEGIN 

man  cbrrori  division  by  taylor  variable*  )> 

WRITELM  ('  BQOAL  TO  IBBO* )/ 

SVR  (0)  (•  RETURN  TO  OFERATING  SYSTEM  *) 

END  (•  THEM  •) 

ELSE  BEGIN  (•  APPLY  L'NOSFITAL'S  ROLE  *> 

IF  0. LENGTH  -  1 
THEM  BEGIN 

WRITE  ('ERROR)  DIVISION  WITH  TAYLOR  VARIABLE* )> 

WRITE LN  ( *  OF  LENGTH  OWE.* )) 

SVR  (0)  (•  RETURN  TO  OFERATING  SYSTEM  •) 

END)  (*  IF  *) 

0. LENGTH  I-  T. LENGTH  -  I) 

FOR  X  i-l  TO  U.  LENGTH  DO 
BEGIN 

IF1  «-  I  +  1) 

O.TCtll  i-  T.TCIXP1J 
END)  (•  FOR  •) 

0  l-  0.0  /  0  <*  THIS  IS  A  RECURSIVE  CALL  •) 

BHD)  (•  IF  «) 

REE  )■  0 

END I  (*  DEAL  /  TAYLOR  *) 
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t 


OPERATOR  /  (Tl  TAYLOR*  Rl  RIAL)  RES  t  TAYLOR*  (*  T  /  R  *) 

VAR  Oi  TAYLOR*  It  DIMTYPE* 

BEGIN 

0. LENGTH  I"  T. LENGTH* 

U.T  »■  T.T* 
ir  R  -  0 

THEM  BEGIN 

WRITE LN  ('ERROR*  DIVISION  OT  TAYLOR  VARIABLE  BY  SERO')> 

BVR  (0)  END  (*  RETURN  TO  OPERATING  SYSTEM  •) 

ELBE  BEGIN 

FOR  I  *-  1  TO  U. LENGTH  DO  D.TCtI)  I-  T.TC(I]  /  R 

END*  (•  ELBE  •) 

RE8  «-  U 

END*  (•  TAYLOR  /  REAL  *) 


OPERATOR  /  (TA,  TBl  TAYLOR)  RE8  t  TAYLOR* 
(*  Raquiraat  VRNULL,  SCALP  *) 


(*  T  /  T  *) 


(*  THIS  OPERATOR  ATTEMPTS  TO  HANDLE  0/0  SITUATIONS  USING  L' HOSPITAL' 8  RULE 
TA.TCM1  -  0.0  <>  0.0 


TB.TCt  1] 


CASE  1 
TA*  /  TB' 

CASE  3 
TA  /  TB 


CASE  2 
ERROR 

CASE  3 
TA  /  TB 


CASE  1  18  HANDLED  BY  RECURSION  UNLESS  THE  SERIES  LENGTH  IS  ONLY  1  •) 

VAR  U.  N  .  TAYLOR* 

I,  J,  IP1,  REV*  DIMTYPE* 

RATIO*  REAL* 


BEGIN 

IP  TA.T  <>  TB.T  (*  UNEQUAL  SCALE  FACTORS  *) 

THEN  BEGIN 

WRITE  ('ERROR*  DIVISION  OF  TAYLOR  VARIABLES')* 
WRITE LM  ('  WITH  UNEQUAL  SCALE  FACTORS')* 

BVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  •) 

END*  (*  IF  •) 

U. LENGTH  *-  TA. LENGTH* 

IF  U. LENGTH  >  TB. LENGTH  THEN  U. LENGTH  **  TB. LENGTH* 

U.T  *■  TA.T* 


U.TC  i-  VRNULL*  W.TC  *-  VRNULL* 


(*  EERO  RESULT  AMD  WORE  VECTORS  •) 


I 


IF  TB.TC(I)  <>  0.0 

ran  begin  (*  cask  3>  usual  divide  *) 

RATIO  I-  1.0  /  TB.TCtl] I 
U.TCtl]  :«  TA.TCtU  *  RATIOl 
FOR  I  l-  2  TO  U. LENGTH  DO 
BEGIN 

FOR  J  «-  1  TO  (1-1)  DO 
BEGIN 

REV  !■  I  -  J  ♦  1 1 
W.TC  [J]  I-  TB.TC(RBV) 

END I  (*  FOR  J  •) 

U.TCtl)  I-  (TA.TCtU  -  SCALP  (U.TC,  N.TC,  0))  *  RATIO 
END  (•  FOR  I  •) 

END  (*  THEN  «> 

ELSE  IF  (TA.TCtU  <>  0.0)  OR  (U. LENGTH  -  1) 

THEN  BEGIN  (*  CASE  It  BOTH  SERIES  IDENTICALLY  ZERO, 

OR  CASE  2t  TA  /  0  *) 

WUTE  (  ‘ERROR i  DIVISION  BY  TAYLOR  VARIABLE* ) ; 

WRXTELN  ('  EQUAL  TO  ZERO')! 

SVR  (0)  {*  RETURN  TO  OPERATING  SYSTEM  *> 

END  (*  THEM  *) 

ELSE  BEGIN  (*  APPLY  L'HOSPITAL'S  RULE  •) 

IF  U. LENGTH  -  1 
THEN  BEGIN 

WRITE  (  *  ERROR  I  DIVISION  WITH  TAYLOR  VARIABLE')! 

WRXTELN  ('  OF  LENGTH  ONE.')! 

SVR  (0)  {*  RETURN  TO  OPERATING  SYSTEM  *) 

END!  (*  IF  •) 

U. LENGTH  t»  TA. LENGTH  -  1/ 

W. LENGTH  t«  U. LENGTH! 

W.T  «-  U.T! 

FOR  I  x-  1  TO  U. LENGTH  DO 
BEGIN 

IP1  !-  I  ♦  l! 

U.TCtl]  I-  TA.TC tIPI] f 
W.TCtl]  t-  TB.TCIIP1] 

END!  (*  FOR  I  •) 

U  l-  U  /  W  (*  THIS  IS  A  RECURSIVE  CALL  •) 

END!  (•  IF  •) 

RES  !•  U 

END!  («  TAYLOR  /  TAYLOR  *) 


OPERATOR  /  (Xs  INTEGER!  Tl  XTAYLOR)  RES  I  ITAYLOR! 
(*  Htquirtu  IVRNULL,  I  SCALP  *) 

VAR  U,  Wt  ITAYLOR! 

I,  J,  XP1,  RHVt  DIMTYPEl 
RATIO!  INTERVAL! 


(*  X  /  IT  •) 


BEGIN 

U. LENGTH  t»  T. LENGTH; 

U.Ti*  T.Ti 

O.TC  i-  IVRNULL!  W.TC  I-  IVRNULL;  (•  ZERO  RESULT  AND  WORK  VECTORS  •) 


(*  NORMAL  DIVIO!  *) 


IP  NOT  (0.0  IN  T.TC[1] ) 

THEM  BIGIN 

RATIO  1  /  T.TCI1] : 

U.TC[1]  «-  K  •  RATIO l 
FOR  I  i-  2  TO  0. LENGTH  DO 
BEGIN 

FOR  J  i-  1  TO  (1-1)  DO 

BEGIN 

REV  I  -  J  ♦  If 
W.TCtJ]  l-  T.TCIREV] 

END)  <*  FOR  J  *) 

U.TC(I]  I-  -  I SCALP  (B.TC,  N.YC,  I)  •  RATIO 
END  (•  FOR  I  *} 

END  (•  THIN  *) 

ELSE  IF  (K  <>  0)  OR  (T.TC[1]  <>  WTPT  (0.0)) 

THEN  BEGIN 

NRITX  ('ERROR:  DIVISION  BY  ITAYLOR  VARIABLE'): 

NRITELM  ('  EQUAL  TO  ZERO'): 

SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  *) 

END  (•  THEM  •) 

ELSE  BEGIN  (*  APPLY  L'HOSPITAL'S  ROLE  •> 

IF  U. LENGTH  -  1 
THEN  BEGIN 

WRITE  ( 'ERROR:  DIVISION  WITH  ITAYLOR  VARIABLE'): 

NRITELM  ('  OF  LENGTH  OMI.‘): 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END:  (*  IF  *) 


U. LENGTH  1-  T. LENGTH  -  1: 

FOR  I  i-1  TO  U. LENGTH  DO 
BEGIN 

IP1  :-  I  ♦  It 
U.TC(I)  :-  T.TC [IP1] 

END:  (•  FOR  I  •) 

U  :-  0  /  U  (•  THIS  IS  A  RECURSIVE  CALL  •) 

END:  (•  IF  *) 

RES  l-  U 

END:  <*  INTEGER  /  ITAYLOR  •) 


OPERATOR  /  (T:  ITAYLOR:  K:  INTEGER)  RES  :  ITAYLOR:  (*  IT  /  K  •) 

VAR  0:  ITAYLOR:  I:  DIMTYPE: 

BEGIN 

U. LENGTH  :-  T. LENGTH: 

U.T  :“  T.T: 

IF  I  •  0 

THEM  BEGIN 

NRITELM  ('ERROR:  DIVISION  OF  ITAYLOR  VARIABLE  BY  EERO'): 

SVR  (0)  END  (•  RETURN  TO  OPERATING  8YSTEN  •) 

ELSE  FOR  I  1  TO  D. LENGTH  DO  U.TCtl]  T.TClI]  /  K| 

RES  :-  U 

END:  (•  ITAYLOR  /  INTEGER  •) 
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<*  I  /  IT  *) 


OPERATOR  /  (Xi  INTERVAL  I  Tl  ITAYLOR)  RX8  l  ITAYLOR; 
{•  Mqulrasi  XVRNULL,  I8CALP  •) 

VAR  a,  Wi  ITAYLOR; 

I,  J,  XU,  REV i  DIMTYPE; 

RATIO!  INTERVAL; 

BROW 

0.  LENGTH  «-  T. LENGTH; 

U.Tl*  T.T; 


O.TC  i-  IVRNULL;  N.TC  !-  IVRWULL;  (•  ZERO  RESULT  AMD  WORK  VECTORS  *) 
ir  NOT  (0.0  IN  T.TC(I)) 

THE*  BEGIN  (•  NORMAL  DIVIDE  *) 

RATIO  !-  1  /  T.TCMIJ 
0.TCI1]  «»  X  •  RATIO; 

FOR  X  l-  2  TO  0. LENGTH  DO 
BEGIN 

FOR  J  !•  1  TO  (1-1?  DO 
BEGIN 

REV  t-  I  -  J  ♦  1; 

W.TCUJ  «-  T. TC [REV] 

END;  (•  FOR  J  •> 

O.TCdl  I-  -  I  SCALP  (O.TC,  N.TC,  I)  *  RATIO 
END  (•  FOR  I  •) 

END  (*  THEN  *) 

ELBE  IP  (X  <>  XNTPT  (0.0))  OR  (T.TC[1)  <>  XNTPT  (0.0)) 

THEN  BEGIN 

WRITE  ('ERROR!  DIVISION  BT  XTATLOR  VARIABLE'); 

NRITELN  ( *  EQUAL  TO  EEROM; 

8VR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  •) 

END  <*  THEN  *) 

ELSE  BEGIN  (•  APPLY  L’ HOSPITAL' 8  RULE  *) 

ir  U. LENGTH  -  1 
THEN  BEGIN 

WRITE  ('ERROR!  DIVISION  WITH  ITAYLOR  VARIABLE'); 

NRITELN  ('  OP  LENGTH  OMR.’); 

8VR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  *) 

END;  (•  IP  *) 


U. LENGTH  l-  T. LENGTH  -  1; 

POM  I  i-1  TO  U. LENGTH  DO 
BEGIN 

IP1  i-  I  ♦  1» 

U.TCIX)  l-  T.TCUP1J 
END;  (*  FOR  •) 

U  i-  0  /  U  (•  THIS  IS  A  RECURSIVE  CALL  •) 

END;  (•  IP  •) 

RES  !-  0 

END;  (•  INTERVAL  /  ITAYLOR  •) 


59- 


(*  IT  /  I  *) 


OPERATOR  /  (Ti  ITAYLOR;  Kt  INTERVAL)  RIS  !  ITAYLOR; 

VAR  Oi  I TAYLOR I  It  DIKTYPE; 

BEGIN 

0. LENGTH  tm  T. LENGTH; 

U.T  t-  T.T; 

IP  0.0  IN  K 
THEN  BEGIN 

WRITE LN  ( ’ERROR)  DIVISION  OP  ITAYLOR  VARIABLE  BY  ZERO'); 

SVR  (0)  END  (•  RETURN  TO  OPERATING  8YSTEM  •> 

ELSE  FOR  I  »-  1  TO  U. LENGTH  DO  D.TCH}  »«  T.TCtI)  /  K> 

RES  I-  U 

END;  (*  ITAYLOR  /  INTERVAL  ») 


OPERATOR  /  (TA,  TBi  ITAYLOR)  RE8  I  ITAYLOR;  (*  IT  /  IT  •) 

(*  Requires i  IVRNULL,  ISCALP  *) 

(*  THIS  OPERATOR  ATTEMPTS  TO  HANDLE  0/0  SITUATIONS  U8ING  L'HOSPITAL'S  RULE 


TA.TCI1]  -  0.0  <>  0.0 


|  CASE  1  CASE  2 

-  0.0  |  TA*  /  TB*  ERROR 

TB.TCt  1]  | 

<>  0.0  j  CASE  3  CASE  3 

I  TA  /  TB  TA  /  TB 


CASE  1  IS  HANDLED  BY  RECURSION  UNLESS  THE  SERIES  LENGTH  IS  ONLY  1  *) 

VAR  U,  W  t  ITAYLOR; 

I,  J,  IP1,  REVt  DIMTYPE; 

RATIO;  INTERVAL; 

BEGIN 

IP  TA.T  <>  TB.T  (•  UNEQUAL  SCALE  FACTORS  *) 

THEN  BEGIN 

WRITE  CERRORt  DIVISION  OP  ITAYLOR  VARIABLES'); 

WRITE LN  (•  WITH  UNEQUAL  SCALE  FACTORS'); 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END;  (*  IP  •) 

U. LENGTH  l-  TA. LENGTH; 

IP  U. LENGTH  >  TB. LENGTH  THEN  U. LENGTH  I-  TB. LENGTH; 

U.T  i"  TA.T; 


U.TC  i-  IVRNULL;  W.TC  I-  IVRNULL;  (*  ZERO  RE8ULT  AND  WORK  VECTORS  *) 


IF  NOT  (0.0  IN  TB.TC(I) ) 

THEM  BEGIN  (•  CAM  3l  USUAL  DIVIDE  *) 

RATIO  I-  1  /  TB.TC[1]| 

D.TCI 1 ]  l-  TA.TCtll  *  RATIO l 
FOR  I  >-  2  TO  U. LENGTH  00 
BBGIM 

FOR  J  «-  1  TO  (1-1)  DO 

BBGIM 

REV  «-  I  -  J  ♦  1» 

M.TC(J)  «-  TB.TC(REV) 

BMDl  (»  FOR  J  •) 

U.TCtlJ  1-  (TA.TCII)  -  ISC ALP  (O.TC,  M.TC,  I)>  *  RATIO 
BHD  (*  FOR  I  «) 

BHD  (•  THEM  *) 

EL SB  IF  (D. LENGTH  -  1)  OR 

(MOT  (0.0  IM  TA.TCflJ))  OR 
(TA.TCdl  <>  IMTPT  (0.0))  OR 
(TB.TC(I)  <>  IMTPT  (0.0)) 

THEM  BBGIM  (•  CASE  It  BOTH  SBRIBS  IDENTICALLY  ZERO, 

OR  CASE  2i  TA  /  0  *) 

WHITE  ( 'ERRORt  DIVISION  BY  I TAILOR  VARIABLE* ) } 

MRITSLM  ('  EQUAL  TO  ZERO’ ) I 

SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  •) 

END  (•  THEM  *) 

ELSE  BEGIN  («  APPLY  L'HOSPZTAL'S  ROLE  *) 

(*  ONLY  IF  IMF  -  SOT  -  0  •) 

IF  a. LENGTH  -  1 
THEM  BEGIN 

WRITE  ('ERROR!  DIVISION  WITH  ITAYLOR  VARIABLE* ) » 

WRITELN  ('  OF  LENGTH  ONE.')} 

SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  •) 

END}  (*  IF  *) 

U. LENGTH  I-  TA. LENGTH  -  1} 

W. LENGTH  U. LENGTH} 

M.T  I"  U.T / 

FOR  I  1-  1  TO  0. LENGTH  DO 
BEGIN 

IP1  !-  I  ♦  1| 

U.TC[Il  »-  TA.TC [IP1] | 

M.TC(I)  I-  TB.TClIPI) 

EMOl  (*  FOR  I  *) 

U  !-  U  /  M  (•  THIS  IS  A  RECURSIVE  CALL  •) 

END;  (*  IF  •) 

RES  i*  U 

END}  (*  ITAYLOR  /  ITAYLOR  •) 


(*  END  OF  RIT_DIV.LIB  »»»»»»»»»»»»»»»»»»»»»  •) 


(*  RT  PON. LIB 


UAL  TAYLOR  F0MXR8 


<<<<<<<<<<<<<<<<<<<<<<<<< 


Crmtantat 

TSQR(T)  FUNCTION  TSQR  (Ti  TAILOR)  t  TAYLOR) 
TSQRT(T)  POMCTIOM  T8QRT  (Tt  TAYLOR)  I  TAILOR) 
TEXP(T)  (UNCTION  TUP  (Tt  TAYLOR)  t  TAILOR; 


K  **  K 
R  **  K 
R  **  R 
K  •*  R 
T  •*  K 
T  •*  R 
R  •*  T 
K  ••  T 
T  **  T 


OPERATOR  ••  (BASE,  EXPONENT  t  INTEGER)  RES  t  INTEGER; 
OPERATOR  **  (BASE:  REAL;  EXPONENT;  INTEGER)  RES  t  REAL: 
OPERATOR  •*  (BASE,  EXPONENT!  REAL)  RES  I  REAL; 

OPERATOR  ••  (BASE:  INTEGER;  EXPONENT)  REAL)  RE8  i  REAL; 
OPERATOR  *«  (BASEi  TAYLOR;  EXPONENT)  INTEGER)  RES  t  TAYLOR; 
OPERATOR  **  (BASE!  TAYLOR;  EXPONENT!  REAL)  RES  t  TAYLOR; 
OPERATOR  “  (BASE!  REAL;  EXPONENT!  TAYLOR)  RES  t  TAYLOR; 
OPERATOR  ••  (BASE!  INTEGER;  EXPONENT!  TAYLOR)  RE8  I  TAILOR; 
OPERATOR  **  (BASE,  EXPONENT!  TAYLOR)  RES  t  TAYLOR; 


■ . . . . — — — . ——————— — —  *) 

FUNCTION  TSQR  (Ti  TAILOR)  I  TAYLOR;  (*  TSQR(T)  •) 

{•  RaqulraS)  VRNULL,  SCALP,  SQR  *) 

VAR  I,  J,  K.  HALF)  D INTYPE; 

X,  It  RVBCTOR; 

O  i  TAYLOR; 

BEGIN 

X  !•  VRNULL;  I  )-  VRNULL; 

U. LENGTH  ]-  T. LENGTH; 

U.T  i-  T.T; 

O.TC  )-  VRNULL; 

U.TCll]  i-  SQR  (T.TClU); 

X(1]  i-  T.TCtlJ; 

FOR  E  I-  2  TO  U. LENGTH  00 
BEGIN 

X(K)  !■  T.TC(X) ; 

SALT  i-  X  DXV  it 
PON  J  1  TO  HALF  SO 

BEGIN 

I  ;•  E  -  S  ♦  1; 

Y(J)  T.TCII] 

END;  (*  FOR  J  *) 

U.TCtXj  i-  2.0  *  SCALP  (  X,  Y,  0); 

IF  X  NOD  2  -  1  THEN 

BEGIN 

HALF  i-  HALF  ♦  t; 

U.TClXjt-  O.TC [X]  ♦  SQR  (T.TCIRALP]  > 

END  (•  IF  *) 

END;  (•  FOR  X  •) 

TEQR  )-  U 

END;  <*  FUNCTION  TSQR  (TAYLOR)  •> 
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(*  TSQRT(T) 


POMCTION  TSQRT  (Tl  TAILOR )  t  TAYLOR) 

(•  R*qulr«si  VRMOLL,  SQRT,  8QR,  SCALP  •) 

VAR  X,  X,  INDEX,  INDX2l  DIKTYPX) 

RATIO I  REAL) 

X,  Y  s  RVXCTORl 
0  >  TAYLOR) 

BEGIN 

IP  (T.TC(1)  <  0.0)  OR 

((T.TCM1  -  0.0)  AND  (T. LENGTH  >-  2)) 

THEN  BEGIN 
NRXTSLN) 

WRXTELN  (’ERROR)  8QOARE  ROOT  Or  TAYLOR  VARIABLE  <-  ZERO.')) 
SVR  (0)  END)  (*  RETORM  TO  OPERATING  SYSTEM  •) 

(•  ELSE  *> 

0. LENGTH  t-  T. LENGTH) 

O.T  »-  T.T; 

X  i-  VRMOLL)  Y  VRNULL) 

O.TC  I-  VRMOLL) 

O.TCtl]  i“  SQRT  (T.TC(1])|  (*  REAL  •) 

IP  O. LENGTH  >-  2  THEN  BEGIN 
RATIO  i-  1.0  /  (2.0  •  O.TCM]  )) 

O.TC(2]  t-  T.TC[2]  *  RATIO) 

IP  U. LENGTH  >-  3  THEM  BEGIN 

O.TC[3]  i-  (T.TC(3]  -  8QR  (O.TC[2]))  *  RATIO) 

POR  K  I-  4  TO  O. LENGTH  DO 
BEGIN 

IP  K  NOD  2-0 
THEN  BEGIN 

INDEX  K  DIV  2) 

INOX 2  i-  INDEX  -  1) 

X  [XNDX2]  )-  O.TC  (INDEX] 

END  (•  THEN  •) 

ELSE  BEGIN 

INDEX  i-  iK  ♦  1)  DIV  2) 

U.TCIK]  t-  8QR  (O.TC(XNDEX) ) 

END)  (•  IP  •) 

POR  I  i-  1  TO  INDX2  DO 
BEGIN 

INDEX  t-  E  -  I) 

Y(X)  >-  O.TC (INDEX)) 

END)  (•  POR  I  *) 

O.TCtK)  i-  (T.TC(K)  -  O.TCIE)  -  2.0  *  SCALP  (X,  Y,  0))  •  RATIO 
END  (*  POR  X  •) 

END  (•  IP  U. LENGTH  >»  3  •) 

END)  (•  IP  0. LENGTH  >-  2  *) 

TSfiRT  x-  0 

END)  (*  PONCTION  TSQRT  (TAYLOR)  •) 


I 


FUNCTION  TEXP  (T«  TAYLOR )  i  TAYLOR;  (•  TEXP(T)  *) 


(♦  Requires i  VRN DLL,  SCALP  •) 
VAR  J,  X,  INDEX  l  DIMTYPE: 

RATIO:  REAL: 

X  :  RVECTOR; 

U  :  TAYLOR: 

BEGIN 

X  j-  VRNULL; 

U. LENGTH  T. LENGTH: 

U.T  :«  T.T: 

U.TC  VRNULL: 

U.TC [  1  ]  i-  EXP  (T.TC[1]>: 


'4 

i 


FOR  K  :-  2  TO  U. LENGTH  DO 
BEGIN 

RATIO  :«  1.0  /  (X  -  1): 

FOR  J  j-  2  TO  X  DO 
BEGIN 

INDEX  :»  X  -  J  ♦  1: 

X [INDEX]  T.TC[J]  *  (J  -  1)  *  RATIO 
END:  (•  FOR  J  *) 

U.TC[E]  SCALP  (U.TC,  X,  0) 

END:  (*  FOR  X  •) 

TEXP  I-  U 

END:  (•  FUNCTION  TEXP  (TAYLOR)  •) 


OPERATOR  **  (BASE,  EXPONENT  :  INTEGER)  RES  :  INTEGER: 
VAR  Ul  INTEGER; 


(*  X  **  X  •) 


BEGIN 

IF  EXPONENT  <  0 
THEN  BEGIN 

WRITELN  ('ERROR:  INTEGER  *•  NEGATIVE  INTEGER'  ); 

SVR  (0)  (•  RETURN  to  operating  system  *) 

END:  (*  IF  *) 

IF  BASE  -  0 

THEN  IF  EXPONENT  -  0 
THEN  BEGIN 

HRITELN  ('ERROR:  ZERO  ••  ZERO’): 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END  (*  THEN  •) 

ELSE  U  :“  0 


(•  NORMAL  BRANCH  *> 


ELSE  BEGIN 
U  :"  1; 

WHILE  EXPONENT  >  0  DO 
BEGZN 

IP  EXPONENT  MOD  2  -  1  THEN  U  :-  U  •  BASE: 
EXPONENT  :>  EXPONENT  DIV  2: 

IF  EXPONENT  >  0  THEM  BASE  i-  SQR  (BASE) 
END  (*  WHILE  *) 

END:  (*  ELSE  *) 


RES  :•  U 
END:  (• 


INTEGER 


INTEGER 


') 


m 


. .. 


'*  V 

□ 

■  •  ■ ; 
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(•  R  ••  K  *) 


OPERATOR  **  (BASK l  REAL)  EXPONENT!  INTEGER)  R*S  t  RIAL} 
VAR  U«  RXAL) 

UK  SXPl  BOOLEAN! 


BEGIN 

ir  BAS*  »  0.0 

■now  ir  exponent  -  o 
THEN  BEGIN 


MUTtUl 

NRZTILN  ('ERROR!  ZERO  **  ZERO')! 

SVR  (0)  (*  RETURN  TO  OPERATING  8Y8TEM  *) 

END  (*  THEN  *) 

ELSE  IP  EXPONENT  <  0 
THEN  BEGIN 
NRITEUtl 

WRITE LH  ('ERROR!  ZERO  '*  NEGATIVE')! 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END  (*  THEN  *) 

ELSE  U  t-  0  (•  ZERO  »*  POSITIVE  *) 


ELSE  IP  BASE  -  1.0  THEN  U  I-  1.0 

ELSE  BEGIN  <*  BASE  O  0  OR  1.  USUAL  BRANCH  *) 


NEG  EXP  I-  EXPONENT  <  0) 

EXPONENT  I-  ABS  (EXPONENT)! 

IP  EXPONENT  -  1  THEN  U  X-  BABE 

ELBE  IP  EXPONENT  -  2  THEM  U  !-  8QR  (BASE) 

ELSE  BEGIN 
U  1“  1 1 

WHILE  EXPONENT  >  0  DO  (•  EXPONENT  -  0  falls  through  •) 

BEGIN 

IP  EXPONENT  MOO  2  -  1  THEM  U  »-  U  *  BASE! 

EXPONENT  !»  EXPONENT  OXV  2| 

IP  EXPONENT  >  0  THEN  BABE  >•>  8QR  (BASE) 

END  <•  WHILE  •} 

END!  (*  ELSE  •) 

IP  NEGJOP  THEN  U  !-  1.0  /  U 
END!  7*  ELSE*) 


RES  !-  U 

END|  (*  REAL  •*  INTEGER*) 


>S' 


(*  R  **  R  *) 


OPERATOR  **  (BASK,  EXPONENT!  REAL)  RES  I  REAL! 
(*  Raqulras:  R**X,  SQR,  EXP,  LN  *) 

VAR  U:  REAL) 

X:  INTEGER! 


BEGIN 

IF  RASE  -  0.0 
THEN  IF  EXPONENT  •  0.0 
THEN  BEGIN 
WRITELN! 

WRITE LN  ('ERROR:  ZERO  "  ZERO’)! 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM 

END  (*  THEN  *) 

ELSE  IF  EXPONENT  <0.0 
THEM  BEGIN 
WRITELN! 

WRITELN  ('ERROR:  ZERO  **  NEGATIVE* )» 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM 

END  (•  THEN  *) 

ELSE  0  0.0  (*  ZERO  **  P08ITIVK  *) 

ELSE  IF  BASE  -  1.0  THEN  U  :-  1.0 

wr aw  BEGIN  (*  BASE  <>  0  OR  1.  USUAL  BRANCH  *) 

X  :-  TRUNC  ( EXPONENT ) | 

U  :«  BASE  «*  X:  (*  USE  REAL  **  INTEGER  *) 

IF  EXPONENT  <>  X 

THEN  IF  BASE  <  0.0 
THEN  BEGIN 
WRITELNf 

WRITELN  ('ERROR:  NEGATIVE  *'  REAL*)! 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  *) 

END  (•  THEN  *) 


(*  USUAL  BRANCH 

ASSERT:  BASE  >  0,  -1  <  EXPONENT  - 

ELSE  U  :-  U  •  EXP  ((EXPONENT  -  X  )  •  LN  (BASE)) 
END I  (*  ELSE  *) 

RXS  :-  U 

END |  (•  REAL  •*  REAL  *) 


OPERATOR  **  (BASE:  INTEGER:  EXPONENT:  REAL)  RES  :  REAL:  (*  X 

(*  R*quir«a:  R»*R,  R**X,  SQR,  EXP,  LN  •) 

VAR  U:  REAL: 

BEGIN 
U  J“  BASE: 

RES  :“  O  »»  EXPONENT  (*  USE  REAL  *'  REAL  •) 

END!  (*  INTEGER  '*  REAL  •) 


•) 


*) 


<  1  *) 


*•  R  *) 
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OPERATOR  *•  (EASE I  TAYLOR 1  EXPONENT)  INTEGER)  RES  i  TAYLOR)  (•  T  ••  E  •) 

(•  Requires*  VRHOLL,  R**X,  T8QR,  T«T,  SCALP,  T_IDENTJBER0(T)  *) 

(*  DECISION  TABLE  Cor  special  cases) 


EXPONENT* 

0 

1 

2 

>  2 

<  0 

I  Box  1 

1 

1 

BASB.TC  - 

0 

I  Ondef . 

-  BASE 

8QR  (BASE) 

—  0 

1 

Ondef.  | 

I  Box  2 

|  BOX  3 

l 

BOX  4 

|  BOX  5 

1 

1 

i 

1 

1 

BA8E.TC [ 1 ]  - 

0 

|  1 

(  -  BASE 

i 

i 

§ 

e 

U 

1 

* 

Ondef.  | 

1 

I 

1 

i 

1 

1  Box  6 

1 

1 

I 

MJ  <> 

0 

j  1 

1 

|  -  BASE 

1 

1 

8QR  (BASE) 

1  1 
|  By  r«curr*nc«  | 

Box  6  is  the  usual  method.  *) 


VAR  D  I  TAYLOR) 

V  I  KVECTOR) 

RATIO  1  REAL) 

K,  J,  REV  <  DIMTYPE) 

BEGIN 

U. LENGTH  »-  8 AES. LENGTH) 

O.T  *■  BASE.Tr 

(•  ZERO  RE8DLT  AMD  WORE  VECTORS  •) 

O.TC  x-  Via* DLL i  V  )«  VRMDLL) 

IF  T  I DENT  ZERO  (BASE)  (*  Decision  table  Box  1 

•  ) 

THEM  IP  EXPONENT  <•  0 
THEN  BEGIN 
NRITELM) 

NRITELN  ('ERROR*  ZERO  **  ZERO')) 

SVR  (0)  (*  RETORN  TO  OPERATING  SYSTEM  *) 

END  (•  THEN  •) 

EijSS  (*  O.TC  Is  already  •“  0  *) 


ELSE 

IF 

EXP0R4NT  »  0 

(* 

Decision 

table 

Box 

2 

*) 

1 

■HE! 

1  O.TC(I)  i-  1.0 

ELSE 

IF 

EXPONENT  -  1 

(« 

Decision 

table 

Box 

3 

*) 

1 

Ml 

1  0  *-  BASE 

ELSE 

IF 

EXPONENT  "  2 

(* 

Decision 

table 

Box 

4 

*) 

TON  0  l-  TSQR  (BASE) 
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(*  Decision  table  Box  5  •) 


I 


ELSE  IF  BASE.TC [ 1 ]  -  0 
THEM  ir  EXPONENT  <  0 
THEN  BEGIN 
WRITELN) 

WRITE LN  (’ERROR!  ZERO  TAYLOR  VARIABLE  **  NEGATIVE' )> 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END  (•  THEN  •) 

ELSE  BEGIN  (*  Use  repeated  multiplications  of  series  *) 

U.TC[1]  s-  1.0» 

WHILE  EXPONENT  >  0  DO 
BEGIN 

IF  EXPONENT  MOD  2  -  1  THEN  U  t«  U  *  BASE; 

EXPONENT  l-  EXPONENT  DIV  2> 

IF  EXPONENT  >  0  THEN  BASE  t-  TSQR  (BASE) 

END  (*  WHILE  *) 

END  (*  ELSE  *) 

(*  Decision  table  Box  6.  Usual  branch  *) 
BASE.TC [11  •*  EXPONENT)  (*  REAL  **  INTEGER  •) 

:«  1.0  /  BASE . TC 1 1 1 1 
IF  U.  LENGTH  >-  2  THEN  BEGIN 

U.TC12]  l-  EXPONENT  *  U.TC[1]  *  BA8E.1C[2]  •  RATIO) 

FOR  K  1-  3  TO  U.  LENGTH  DO 
BEGIN 

FOR  J  !-  1  TO  K  -  1  DO 
BEGIN 

REV  S-  X  -  J  ♦  1) 

V[J]  i-  BASE.TC [REV]  •  (EXPONENT  *  (K  -  J)  -  J  +  1>) 

END)  (•  rOR  J  *) 

U.TC [K]  i-  RATIO  •  (EXPONENT  *  SCALP  (U.TC,  V,  0) 

END  (*  FOR  K  *1 
END  (•  IF  U. LENGTH  >“  2  •) 

END)  (•  ELSE  •) 

RES  J-  U 

END)  (•  TAYLOR  ••  INTEGER  •) 


ELSE  BEGIN 
U.TC[1] 
RATIO 


(•  T  ••  R  *) 

OPERATOR  **  (BASE)  TAYLOR)  EXPONENT)  REAL)  RES  !  TAYLOR) 

(•  Requires!  VRNULL,  T*»K,  R»*R,  R**K,  T*T,  TSQR,  TSQRT,  SCALP, 
T_IDENT_ZERO,  T_IDENT_CON8TANT  *) 

(•  DECISION  TABLE 

EXPONENT) 


BA8E.TC  -  0 
BASE.TCM]  <-  0 


BASE.TC [1]  >  0 

Boxes  4  or  6  are  the  usual  method.  *) 


for  special  cases: 

Not  Int. 


0 

1 

2  or  1/2 

Int.  >  2 

or  < 

Box  1 

1 

1 

BOX  2 

|  Box  3 

1 

|  Box  4 

|  BOX  5 

Undef . 

1 

1 

1 

-  BASE 

i 

|  SQR  or  SQRT 

l 

|  SQR  or  SQRT 

1 

|  —  0 
| 

1 

I  Undef 

1 

1 

1 

-  BASE 

|  T  **  K 

1 

I  Undof 

1 

I 

1 

1 

1 

•  BASE 

1 

1 

1 

1  SQR  or  SQRT 

1 

1 

1 

|  T  **  K 

1  Box  6 

1 

I  recur 

♦ - - - - - - ♦- - - - + - — - * 


VAR  a  t  TAYLOR: 

V  t  RVECTORl 

RATIO  i  REAL I 

X  i  INTEGER: 

X,  J,  REV  t  OIMTYPE: 

BEGIN 

D.UMn  l>  BASS. LENGTH: 

O.T  «-  BASE. T I 

{*  ZERO  RESULT  AMO  WORK  VECTORS  •) 

O.TC  l-  VRHULL:  V  l-  VRH VLLl 

IP  exponent  -  0.0  (*  Dacialon  tabla  Box  1  •) 

THIN  XT  T  IDEWT  ZERO  (BASE) 

THEM  BEGIN 
NRIRLM; 

WRXTELM  ('ERROR:  ZERO  *•  ZERO' ) t 

BVR  (0)  (•  RETORM  TO  OPERATXMG  SYSTEM  •) 

END  (*  THEM  *) 

ELSE  O.TCI1]  »-  1.0 

ELSE  It  EXPOMEWT  “  1.0  (*  Dacialon  Ubli  Box  2  *) 

THEM  O  i-  BABE 

ELSE  IP  EXPOMEWT  »  2.0  (*  Dacialon  tabla  Box  3  •) 

THEM  O  «•  TSQR  (BASE) 

ELSE  IP  EXPONENT  -  O.S  (•  Doclalon  tabla  Box  3  •) 

THEM  D  »-  TSQRT  (BASE)  (•  Error  1C  BASE.TC(I)  <-  0.0  *> 

XU(E  BEGIN 

I  I-  TRDMC  (EXPONENT) I  (*  Dacialon  tabla  Box  4  *) 

IP  (EXPONENT  ■  X)  AMD  (I  >  2) 

THEN  0  «•  BASE  **  I  (*  TAYLOR  ••  INTEGER  *) 

ELSE  IP  BABE.TC t 1]  <>  0.0 

THEM  BEGIN  (*  Dacialon  tabla  Box  5  •) 

NRITSLMl 

MUTE  ('ERROR:  ZERO  TAYLOR  VARIABLE  "  NEGATIVE'): 

WRITELM  ('  OR  NON-INTEGER’): 

8VR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END  (•  THEN  *) 

ELSE  BEGIN  (*  Dacialon  tabla  Box  6.  Usual  branch  ') 

O.TCdl  BASE. TC Ml  •*  EXPONENT:  (*  REAL  *•  REAL  •) 

RATIO  i-  1.0  /  BASE.TCCl] I 
IP  0. LENGTH  >-  2  THEM  BEGIN 

O.TC(2]  :-  EXPONENT  *  U.TC(t)  *  BASE. TC 12]  *  RATIO: 

POR  E  :-  3  TO  D. LENGTH  DO 
REGXM 

POR  J  1  TO  X  -  1  DO 
BEGIN 

REV  i-  X  -  J  ♦  1: 

V(J]  i-  BASE.TC [REV]  •  (EXPONENT  •  (K  -  J)  -  J  ♦  1) : 

END:  (*  POR  J  •) 

U.TC(K)  «•  RATIO  •  SCALP  (O.TC,  V,  0)  /  (K  -  1) 

END  (•  POR  K  •) 

END  (•  IP  U. LENGTH  >-  2  •) 

END  (•  ELSE  •) 

END:  (•  BLEB  •) 

RES  U 

END:  (*  TAYLOR  •*  REAL  •) 
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(*  R  *•  *  •) 


OPERATOR  **  (BASE;  REAL;  EXPONENT;  TAILOR)  US  i  TAILOR; 

{•  Raquiraa;  vrmull,  r**r,  r**e,  SQR,  lm,  scalp  •) 

TAR  J,  X,  INDEX;  DIHTIPt; 

RATIO,  LOG  BASRs  RIAL; 

X  i  EVBCTOR; 

0  i  TAILOR; 

BRGIM 

X  i-  VRMTJLL; 

0. LENGTH  r-  EXPONENT . LENGTH ; 

O.T  1-  EXPONENT. T; 

D.TC  t-  VRNDLL; 

tJ.TCM)  t-  BABB  *•  EXPONENT. TC[1) ; 

(*  RIAL  **  RIAL  -  may  gmarat*  error*  •) 

IP  MOT  T  I  DINT  CONSTANT  (XXPOMRNT)  THEM  BRGIM 
IP  BASS  <-“0.0 
TURN  BRGIM 
WRITE  LM  | 

WRITRLM  ( ' ERROR  Dl  URAL  ••  TAILOR;  BARR  <-  ZERO'); 

SPR  (0)  (•  RRTOBM  TO  OPRRATIMG  SYSTEM  •) 

RMD  {•  THEM  •) 

RLSR  BRGIM 

LOG  BABB  I-  LM  (BA SB); 

FObTk  i-  2  TO  U. LENGTH  DO 
BRGIM 

RATIO  I-  1.0  /  (X  -  1); 

POM  J  I-  2  TO  X  DO 

BRGIM 

INDEX  l«  X  -  J  ♦  1; 

It  INDEX)  ;*  EXPONENT.  TC  [J]  •  (J  -  1)  •  RATIO 
RMD;  (*  POE  J  •) 

U.TC [X]  I-  L0G_BA8R  •  SCALP  (D.TC,  X,  0) 

RMD  (•  FOR  E  ”) 

RMD  (•  RLSR  •) 

RMD;  (•  IP  •) 

RES  ;■  0 

END;  <*  REAL  *•  TAILOR  *) 


OPERATOR  ••  (BASRl  INTEGER;  EXPONENT;  TAILOR)  RES  1  TAILOR;  (•  X  •*  T  *) 
(*  Raqnlraa;  VRMULL,  R»*T,  R«#R,  R**K,  SQR,  LM,  SCALP  •) 

VAR  0;  REAL; 

BRGIM 
0  I-  BASE; 

RES  1-  0  *•  EXPONENT  <•  USE  REAL  **  TAILOR  *) 

END;  (*  INTEGER  ••  TAILOR  •) 


C.5.  Interval  Taylor  Power  Operators  and  function*. 


(•  IT_POW.LI8 

-  INTERVAL 

TAYLOR  POWERS 

<<<<<<<<<<<<<<<<<<<<<< 

Contents: 

ITSQR (IT) 

FUNCTION 

ITSQR  (T: 

ITAYLOR) 

:  ITAYLOR; 

ITSQRT ( IT ) 

FUNCTION 

ITSQRT  (T: 

ITAYLOR) 

:  ITAYLOR; 

ITEXP(IT) 

FUNCTION 

ITEXP  (T: 

ITAYLOR ) 

:  ITAYLOR; 

I  **  K 

OPERATOR 

•  * 

(BASE: 

INTERVAL; 

EXPONENT:  INTEGER) 

RES  : 

INTERVAL; 

K  •*  I 

OPERATOR 

•  # 

(BASE: 

INTEGER; 

EXPONENT:  INTERVAL) 

RES  : 

INTERVAL; 

I  ••  I 

OPERATOR 

** 

(BASE: 

INTERVAL; 

EXPONENT:  INTERVAL)  RES 

:  INTERVAL; 

IT  **  K 

OPERATOR 

** 

(BASE: 

ITAYLOR; 

EXPONENT:  INTEGER) 

RES  : 

ITAYLOR; 

IT  ••  I 

OPERATOR 

•  * 

(BASE: 

ITAYLOR; 

EXPONENT:  INTERVAL) 

RES  : 

ITAYLOR; 

K  *•  IT 

OPERATOR 

** 

(BASE: 

INTEGER; 

EXPONENT:  ITAYLOR) 

RES  : 

ITAYLOR; 

I  *•  IT 

OPERATOR 

(BASE: 

INTERVAL; 

EXPONENT:  ITAYLOR) 

WS  : 

ITAYLOR; 

IT  **  IT 

OPERATOR 

** 

(BASE: 

ITAYLOR; 

EXPONENT:  ITAYLOR) 

RES  : 

ITAYLOR; 

FUNCTION  ITSQR  (T:  ITAYLOR )  :  ITAYLOR;  (*  ITSQR  (IT)  •) 

(•  Requires:  IVRNULL,  ISC  ALP,  IT_IDENT_CONSTANT(IT),  ISQR  •) 

VAR  I,  J,  K,  HALF:  DIMTYPE ) 

X,  Y:  I  VECTOR; 

U  :  ITAYLOR; 

BEGIN 

X  :-  IVRNULL;  Y  :-  IVRNULL; 

U. LENGTH  :»  T. LENGTH; 

U.T  :“  T-T; 

U.TC  :»  IVRNULL; 

U.TC[1]  :-  ISQR  (T.TC [ 1 ] ) ; 

IF  NOT  IT_IDENT_CONSTANT  (T)  THEN  BEGIN 
X[1]  :-  T.TC [11; 

FOR  K  :-  2  TO  U. LENGTH  DO 
BEGIN 

X  [K]  :-  T.TC  [K] ; 

HALF  :-  K  DIV  2; 

FOR  J  :*  1  TO  HALF  DO 
BEGIN 

I  :-  K  -  J  +  1; 

Y  [J]  T.TC [I] 

END;  (*  FOR  J  •) 

U.TC (K]  2  *  I SCALP  (  X,  Y,  HALF); 

IF  K  MOD  2  -  1  THEN 
BEGIN 

HALF  :-  HALF  +  1; 

U.TC(K):-  U.TC (KJ  ♦  ISQR  (T.TC (HALF)  ) 

END  (•  IF  •) 

END  (•  FOR  K  •) 

END;  (*  IF  NOT  IT_I DENT_CONSTANT  (T)  *) 

ITSQR  :«  U 

END;  (*  FUNCTION  ITSQR  (ITAYLOR)  •) 
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FUNCTION  IT8QRT  (Tt  ITAYLOR)  <  ITAYLOR;  (*  ITSQRT ( IT )  * 

(*  Raquirass  IVRNULL;  ISQRT,  IT_IDXNT  CONSTANT (IT), 

I9QR,  I SCALP  *)  “ 

VAR  I,  K,  INDEX,  INDX2:  DIKTYPS; 

RATIO i  INTERVAL; 

X,  Y  <  I VECTOR; 

U  t  I  TAYLOR  I 

BEGIN 

IF  (IINF  (T.TCtl])  <  0.0)  OR 

( (I  INF  (T.TC(IJ)  -  0.0)  AND  (T. LENGTH  >-  2)) 

THEN  BEGIN 
WRITELNl 

NRITELN  ('ERROR;  SQUARE  ROOT  OF  ITAYLOR  VARIABLE  <-  ZERO.') I 
SVR  (0)  END;  (•  RETURN  TO  OPERATING  SYSTEM  •) 

(*  •) 

X  IVRNULL;  Y  t-  IVRNULL; 

U. LENGTH  !“  T. LENGTH; 

U.T  ;»  T.T; 

U.TC  ;»  IVRNULL; 

U .TC  1 1 )  »-  ISQRT  (T.TCC11);  (*  INTERVAL  *) 

IF  NOT  IT_IOENT_CON8TANT  (T)  THEN  BEGIN 
RATIO  «-  1  /  (2  •  U.TC[1] ); 

IF  U. LENGTH  >-  2  THEN  BEGIN 
U.TC(2]  l»  T.TC12]  •  RATIO; 

IF  U. LENGTH  >-  3  THEN  BEGIN 

U.TC13)  (T.TCI31  -  I8QR  (U.TC(2]>)  •  RATIO; 

FOR  X  I-  4  TO  U. LENGTH  DO 
BEGIN 

IF  K  MOD  2  •  0 
THEN  BEGIN 

INDEX  j-  K  DIV  2; 

INDX2  I-  INDEX  -  t; 

X (INDX2)  U.TC (INDEX] 

END  («  THEN  *) 

ELSE  BEGIN 

INDEX  (X  ♦  1)  DIV  2; 

U.TC (X)  i-  XSQR  (U.TC [INDEX] ) 

END;  (*  IF  *) 

FOR  I  1  TO  INDX2  DO 
BEGIN 

INDEX  K  -  I; 

Y [I]  I-  U.TC  (INDEX]; 

END;  (•  FOR  I  *) 

U.TC(K]  (U.TC (X]  -  2  *  ISCALP  (X,  Y,  INDX2 ) )  »  RATIO 
END  (•  FOR  K  •) 

END  (•  IF  U. LENGTH  >-  3  *) 

END  («  IF  U. LENGTH  >-  2  *) 

END;  (*  IF  NOT  IT_  IDENT_CONSTANT  (T)  •) 

ITSQRT  U 

END;  (*  FUNCTION  ITSQRT  (ITAYLOR)  *) 
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t*  imp  (it)  *) 


FUNCTION  mXP  (T;  ITATLOR)  t  ITATLOR; 

(*  tequlraa;  IVRRULL,  IT_IDENT_CONSTANT(  IT) ,  I  SCALP  •) 
VAR  J,  X,  INDEX i  DIMTYPI;  ”  ~ 

RATIO i  INTERVAL; 

X  l  IVECTOR) 

0  l  ITATLOR; 

BEGIN 

X  l-  IVRRULL; 

0. LENGTH  «■  T. LENGTH; 

U.T  l-  T.T; 

O.TC  I-  IVRRULL; 

o.tc[i]  i-  axp  (T.TCtU); 

IP  NOT  IT_IDENT_CONSTAirr  (T)  THEN  BEGIN 
FOR  I  !■  T  TO  oTlENGTH  DO 

BEGIN 

RATIO  l-  1  /  INTPT  (X  -  1))l 
FOR  J  l»  2  TO  K  DO 
BEGIN 

INDEX  t-  X  -  J  +  1; 

X[ INDEX]  i-  T.TCU1  •  (J  -  1)  *  RATIO 
END;  (*  FOR  J  *) 

O.TC(X)  I-  ISCALP  (O.TC,  X,  X) 

END  (•  FOR  X  *) 

END;  (•  IF  HOT  IT_I DENT JMN8T ANT  (T)  *> 

ITEXP  ;«  O 

END;  (•  FUNCTION  ITEXP  (ITATLOR)  *) 


(*  I  **  X  •) 

OPERATOR  ••  (BASE;  INTERVAL;  EXPONENT;  INTEGER)  RES  I  INTERVAL; 

(*  Require*;  X**X  •) 

VAR  Ut  INTERVAL; 

NEGJEXP;  BOOLEAN; 

BEGIN 


(*  BASE  —  INTEGER  and  EXPONENT  >-  0.  Os*  INTEGER  **  INTEGER  *) 
IF  (BASK  -  INTPT  (TRUNC  (ISOP  (BASE))))  AND 
(EXPONENT  >“  0) 

THEN  0  INTPT  (TRUNC  (ISOP  (BASE))  ••  EXPONENT) 

EL8E  IF  (0.0  IN  BASE)  AND  (EXPONENT  <-  0) 

THEN  BEGIN 
MRXTXLN; 

NRITELN  ('ERROR  IN  INTERVAL  *•  INTEGER;  ZERO  "  <-  EERO') I 

SVR  (0)  (•  RETORN  TO  OPERATING  STSTEN  •) 

END  (•  THEN  *) 


nn  begin 

me  to  i-  exponent  <  o» 
recount  «-  ua  (exponent)! 

IP  monRI  -  1  TUN  0  I-  BASK 

SLAB  ir  EXPONENT  -  2  TUN  0  t-  ISQR  (BABB) 

KLSB  BEGIN 

a  »-  nam  o.o)i 

WHILE  EXPONENT  >  0  DO  (•  EXPONENT 

BEGIN 

IP  EXPONENT  MOO  2  -  1  TUN  U  I-  U  •  BASE! 
EXPONENT  i-  EXPONENT  DIV  2| 

IP  EXPONENT  >  0  TUN  BABE  I-  ISQR  (BASE) 
END  (*  WHILE  •) 

END*  (*  ELBE  *) 

IP  MBS  EXP  TUN  0  i-  t  /  0 
END)  7*  ELSE*) 

RES  «•  0 

END I  (•  INTERVAL  ••  INTEGER  *> 


(*  uaual  branch  *) 


-  0  falls  through  •) 


(•  K  **  I  *> 


OPERATOR  *•  (BASE I  INTEGER!  EXPONENT!  INTERVAL)  RES  i  INTERVAL! 
(*  Raquirsai  K**K  •) 

VAR  Os  INTERVAL! 


BBGIN 

IP  BA8E  -  1  THEN  0  I-  INTPT  (1.0) 


—  1  *) 


(•  EXPONENT  —  INTEGER  »  0.  Oss  INTEGER  *'  INTEGER  •) 
ELBE  IP  (EXPONENT  -  INTPT  (TRDNC  ( I SUP  (EXPONENT))))  AND 
(I80P  (EXPONENT)  >-  0.0) 

THEN  O  I-  INTPT  (BABE  ••  TRDNC  ( IBOP  (EXPONENT))) 


ELSE  IP  BASE  <-  0  (•  EXPONENT  contains  at  lsast  on#  rsal  n unbar  *) 

THEN  BBGIN 
WRITELMl 

WRITBLN  ('ERROR  IN  INTEGER  *•  INTERVAL!  BABE  <-  SERO')» 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END  <*  TUN  •) 


ELBE  U  i-  IEXP  (EXPONENT  *  ILN  (INTPT  (BABE)))! 


(•  Usual  branch  *) 


END!  (*  INTEGER  ••  INTERVAL  •) 


)ir-'  -j/y  •  *. 


(*  I  **  I  •) 

OPERATOR  ••  (BASEl  INTERVAL;  EXPONENT i  INTERVAL)  REE  t  INTERVAL! 

(*  Requires i  X**K,  I**X,  K*»I  *> 

VAR  U:  INTERVAL; 

BEGIN 

IF  BASE  -  INTPT  (1.0)  THEN  0  l-  BASE  (*  BASE  —  1  *) 

(*  EXPONENT  —  INTEGER.  Dee  INTERVAL  **  INTEGER  •) 
ELBE  IP  EXPONENT  -  INTPT  (TRDNC  (ISUP  (EXPONENT) ) ) 

THEM  0  »-  BASE  **  TRUNC  (ISUP  (EXPONENT)) 

(•  BASE  —  INTEGER.  Use  INTEGER  **  INTERVAL  •) 
ELSE  IF  BASE  -  INTPT  (TRUNC  (ISUP  (BASE))) 

THEN  U  :-  TRUNC  (ISUP  (BASE))  »*  EXPONENT 

(*  BASE  end  EXPONENT  each  contain  at  least  one  real  niMber .  *) 

ELSE  IP  I INF  (BASE)  <-0.0 
THEN  BEGIN 
WRITELMt 

WRITEUt  ('ERROR  IN  INTERVAL  **  INTERVAL;  BASE  <-  ZERO'); 

SVR  (0)  (•  RETURN  TO  OPERATING  SY8TEM  •) 

END  (•  THEN  •) 

ELSE  U  I-  IEXP  (EXPONENT  •  I  Lit  (BASE));  (•  Usual  branch  *> 

RES  >«  U 

END;  (*  INTERVAL  '*  INTERVAL  •) 


(*  IT  *•  K  •) 

OPERATOR  ••  (BASEl  ITAYLORl  EXPONENT!  INTEGER)  RES  I  ITAYLOR; 

(•  Requires:  IVRNULL,  I*»K,  IT8QR,  I*IT,  IT*IT,  I SCALP, 

IT  IDENT  ZERO(IT)  •) 


(*  DECISION  table  (or  special  caaesi 


EXPONENT: 

0 

1 

2 

>2  <0 

1 

| 

BOX  1 

1 

BASE.TC  -  0 

! 

Undef. 

•  BASE 

SQR  (BASE) 

■»  0  Undef.  I 

! 

1 

Box  2 

I  Box  3 

i 

i 

i 

i 

I 

Box  4 

1 

BOX  5  I 

0  IN  BASE.TCI1] 

J 

1 

1 

1 

I  -  BASK 

1 

SQR  (BABE) 

1 

1 

1 

By  nult.  Undef.  | 

1 

1 

i 

| 

i 

i 

1 

Box  6  I 

BASE.TC [1]  <>  0 

1 

+« 

1 

1 

1  -  BASK 

- 

1 

1 

SQR  (BABB) 

1 

1 

1 

By  recurrence  | 

Box  6  is  the  usual  Method.  •) 


VAX  D  s 

V  « 

RATIO  i 
K,  J,  REV 


ITAYLOR; 

I VICTOR; 
INTERVAL; 

i  D INTYPE; 


BIG  IK 

U.LENOTH  i-  BASE. LENGTH; 

U.T  i-  BASB.T; 

(•  ZBRO  RESULT  AMD  WORK  VECTORS  *) 
0.1C  i-  IVRNULL;  V  t-  IVRNULL; 


•) 


ir  IT_IDEHT_ZERO  (BASE)  (*  Decision  table  BOX  1 

THEM  IF  EXPONENT  <-  0 
THEM  BEGIN 
WRITHE; 

WRITEIM  ('ERROR;  ZERO  **  ZERO'); 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  *) 

END  (*  THEN  •) 

ELSE  (*  U.TC  is  already  «”  0  *) 


if  EXPONENT  •  0  (*  Decieion  table  Box  2  *) 

THEN  O.TCll]  i-  INTPT  (t.0) 

ELSE  IP  EXPONENT  -  1 
ran  U  I-  BASE 

ELBE  IF  EXPONENT  «  2 

ran  U  I-  ITSQR  (BASE) 

ELSE  IF  0.0  IN  BASI.TC  ( 1] 
ran  IF  EXPONENT  <  0 

ran  begin 
WRITELN; 

WRITBLM  ('ERROR;  ZERO  ITAYLOR  VARIABLE  "  NEGATIVE' ); 

SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  *) 

end  (•  ran  *) 

ELSE  BSCIH  (*  Use  repeated  Multiplications  of  series  *) 

U.TC[1)  «-  INTPT  (1.0); 

WHILE  EXPONENT  >  0  DO 
BEGIN 

IF  EXPONENT  MOD  2  -  1  TITO  U  i«  U  •  BASE; 

EXPONENT  «■  EXPONENT  DIV  2; 

IF  EXPONENT  >  0  THEN  BASE  t-  ITSQR  (BASE) 

END  (*  WHILE  •) 

END  <*  E LSI  •> 

(*  Decision  table  Box  6.  Uaual  branch  *) 
;-  BABE. TC (11  ••  EXPONENT;  (*  INTERVAL  ••  INTEGER  •) 

1-  1  /  BASE.TCIH; 

IF  U. LENGTH  >-  2  ran  BEGIN 

U.TCC21  ;-  (EXPONENT  *  RATIO)  •  U.TC[1)  •  BASI.TC(2]| 


ELSE  BEGIN 

O.TCfl] 

RATIO 


(*  Daciaton  tabla  Box  3  *) 
{•  Daciaton  tabla  Box  4  •) 
(*  Daciaton  table  Box  5  *1 
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FOR  K  t-  3  TO  U. LENGTH  DO 
BIG  XU 

FOR  J  i*  1  TO  It  *  1  DO 

BIGIN 

RXV  :-  K  -  J  ♦  1» 

V(J)  I-  BASI.TC [RBV]  •  (EXPONENT  •  (X  -  J)  -  J  ♦  1); 
BHD |  (*  FOR  J  •) 

U.TC{X]  :-  RATIO  *  X SCALP  (O.TC,  V,  X)  /  (X  -  1) 

BID  (•  FOR  X  •) 

BHD  (*  IF  U. LENGTH  >-  2  •) 

END I  (•  ELSE  *) 

RES  :-  U 

END;  (•  ITAYLOR  ••  INTIGER  *) 


(•  IT  **  X  •) 

OPERATOR  *'  (BASE;  ITAYLOR;  EXPONENT:  INTERVAL)  RES  i  ITAYLOR; 

(*  Raqulraa:  IVRNULL,  I**X,  ITSQR,  I*IT,  IT*IT,  IT**X,  1**1, 

X**I,  I  SCALP,  IT_IDENT_ZBRO(  IT)  ,  XT_XDXMT_COHSTANT(  IT)  •) 

VAR  D  I  ITAYLOR  I 

V  I  I VECTOR; 

RATIO  :  INTERVAL; 

X,  J,  REV  :  DIMTYPEl 

BEGIN 

D. LENGTH  I-  BASE. LENGTH; 

U.T  BASE.T; 

(*  ZERO  RESULT  AND  NORX  VECTORS  *) 

U.TC  :-  IVRNULL;  V  :■  IVRNULL; 

IF  IT_IDENT_ZERO  (BASE) 

THEN  IF~ I  INF  (EXPONENT)  <»  0.0 
THIN  BEGIN 
NRITELN; 

NRITELN  ('ERROR:  ZERO  **  ZERO'); 

SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  •) 

END  <•  THEN  •) 

ELSE  (•  U.TC  ia  already  —  0  •) 

(•  EXPONENT  —  INTEGER.  Uaa  ITAYLOR  **  INTEGER  •) 
ELSE  IF  EXPONENT  -  XNTPT  (TRUNC  (I8UP  (EXPONENT))) 

THEN  U  :-  BASE  **  TRUNC  (ISUP  (EXPONENT)) 

ELSE  IF  I INF  (BA8E.TC(1] )  <-  0.0 
THEN  BEGIN 

NRITELN; 

NRITELN  ('ERROR  IN  ITAYLOR  ••  INTERVAL:  BASE  <-  ZERO.'); 

SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  *) 

END  (•  THEN  •) 


(*  tac uralon.  Muni  branch  •) 

(«  i mini  **  interval  *) 


ELSE  BEGIN 

O.TCtl]  I-  BA8E.TCI1)  **  EXPONENT* 

IF  MOT  IT_I DENT_CONSTANT  (BASE)  TEEM 

begin” 

RATIO  t-  1  /  BASE.TCCUl 

IF  0.  LENGTH  >-  2  THEM  BEGIN 

U.TC[2]  l-  (EXPONENT  *  RATIO)  *  O.TCll]  *  BASE.TC[2]f 

FOR  X  I-  3  TO  0. LENGTH  DO 
BEGIN 

FOR  J  I-  1  TO  X  -  1  DO 
BEGIN 

RET  I-  X  -  J  ♦  1» 

V(J)  i-  BASE.TCtREV]  *  (EXPONENT  *  (X  -  J)  -  J  ♦  1)» 
END i  (•  FOR  J  *) 

O.TC(K)  I-  RATIO  *  I SCALP  (O.TC,  V,  0)  /  (X  -  1) 

END  (•  FOR  K  •) 

END  <*  IF  U. LENGTH  >-  2  *) 

END  (•  IF  NOT  IT_IDENT_CONSTANT  (BASE)  *) 

ENDt  (*  ELSE  •)  ”  “ 

RES  t*  U 

END)  (*  ITAYLOR  **  INTERVAL  *) 


(•  X  •*  IT  *) 

OPERATOR  «•  (BASEi  INTEGER i  EXPONEMTl  ITAYLOR)  RES  I  ITAYLORi 
(•  Raqulraai  IVRNOLL,  X**I,  X**X,  UK. 

IT_IDENT_COM8TANT(IT),  I  SCALP  *) 

VAR  J,  X,  INDEX i  DIMTYPXi 

RATIO,  L0G_BA8Xi  INTERVAL  j 
X  i  I  VECTOR  I 

0  I  ITAYLORI 

BEGIN 

X  I-  IVRNOLL | 

O. LENGTH  i-  EXPONENT. LENGTH | 

O.T  I-  EXPONENT. T| 

O.TC  I-  IVRNOLL  I 

O.TC(I)  t“  BASE  ••  EXPONENT. TC [11 1 

(•  INTEGER  **  INTERVAL  -  aay  vanarnt*  arrora  *) 

IF  NOT  IT_IDENT_CON8TANT  (EXPONENT)  THEM  BEGIN 
IF  BASE  <»  0 
THEN  BEGIN 
MRITXLMl 

NRITXLN  ('ERROR  IN  INTEGER  •*  ITAYLORi  BASE  <»  ZERO' ) l 
8VR  (0)  (•  RETORN  TO  OPERATING  SYSTEM  •) 

END  (•  THEN  •) 

ELSE  IF  BASE  -  1  THEN  U.TC(I)  I-  INTPT  (1.0) 


•7»- 


pi 


ELSE  BEGIN 

LOG_BA8E  1-  ILM  (INTPT  (BASE) ); 

FOR- K  2  TO  0. LENGTH  DO 

BEGIN 

RATIO  «-  1  /  INTPT  <K  -  1); 

FOR  J  :«  2  TO  K  DO 
BEGIN 

INDEX  i-  K  -  J  ♦  1; 

X  [INDEX]  :-  EXPONENT.  TO  (J)  *  (J  -  1)  *  RATIO 
END ;  (*  FOR  J  *) 

U.TCIK]  i-  L0G_BASE  *  I SCALP  (D.TC,  X,  K) 

END  (*  FOR  X  7) 

SID  (•  ELSE  •) 

END;  (•  IF  *) 

RES  I-  V 

END;  (*  INTEGER  **  ITATLOR  *) 


I  **  IT 


OPERATOR  **  (BASE:  INTERVAL;  EXPONENT:  ITATLOR)  RES  I  ITAYLOR; 
(*  Requires:  IVRNULL,  X**I,  K**K,  ILN,  I**I,  I»*K, 
IT_IDENT_CONSTANT(  IT)  ,  I  SCALP  •> 

VAR  J,  K,  INDEX:  D INTYPE l 

RATIO.  LOGBASE :  INTERVAL) 

X  :  IVECTOR; 

U  :  ITATLOR; 


BEGIN 

X  :• 

U. LENGTH  :• 

O.T  :■ 

O.TC  :■ 

U.TC[1]  :« 


IVRNULL; 

EXPONENT  .  LENGTH ; 

EXPONENT. T; 

IVRNULL; 

BASE  **  EXPONENT . TO [ 1 ) ) 

(•  INTERVAL  **  INTERVAL  -  uy  generate  errors  *) 

(*  Note  that  there  is  neither  accuracy  or  speed  advantage  to 
be  gained  by  treating  BASE  INTEGER  using  K**IT 
because  the  BASE  is  first  converted  to  an  interval 
there.  •) 


IF  NOT  IT_IDENT_CONSTANT  (EXPONENT)  THEN  BEGIN 
IF  IINF  (BASE)  <-  0.0 
THEN  BEGIN 
MRITELN; 

WRITELN  ('ERROR  IN  INTERVAL  **  ITATLOR)  BASE  <-  ZERO'); 
SVR  (0)  (•  RETURN  TO  OPERATING  SYSTEM  *) 

END  (•  THEN  •) 


ELSE  IF  (IINF  (BASE)  -  ISUP  (BASS))  AND 
(IINF  (BASE)  -  1.0) 

THEN  U.TC [ 1 J  :•  INTPT  (1.0) 


1.0  *) 


I 
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ELSE  BEGIN 

LOG  BASE  tm  ILN  (BASE); 

POR*K  i-  2  TO  0. LENGTH  DO 
BEGIN 

RATIO  l-  1  /  INTPT  (E  -  1>» 

TOR  J  i-  2  TO  K  DO 
BEGIN 

INDEX  «-  K  -  J  +  1; 

X [INDEX]  I-  EXPONENT. TC[J]  *  (J  -  1)  *  RATIO 
END;  (*  TOR  J  •) 

U.TC [K]  i-  LOG_BASS  *  I SCALP  (O.TC,  X,  X) 

END  (•  TOR  X  *’) 

END  (*  ELSE  •) 

END;  {*  IT  *> 

RES  t-  U 

END;  <*  INTERVAL  «*  I TAYLOR  *) 


OPERATOR  ••  (BASE,  EXPONENT)  ITAYLOR)  RES  I  ITAYLOR;  (*  IT  **  IT  *) 

(*  Raqulr«a>  IT**I,  I**IT,  1**1,  X**I,  I**X,  X**X, 

I  TUI,  ITXXP,  IT8QR,  I*IT,  IT*IT, 

IVRNULL,  I  SCALP,  IT_ZDENT_CONSTANT  (IT)  *) 

VAR  Ut  ITAYLOR; 

BEGIN 

IP  IT  IDEWTJMNSTAMT  (EXPONENT) 

tun  0  7-  BASE  *«  EXPONENT. TC [  1  ]  (•  Os*  ITAYLOR  **  INTERVAL  *) 

ELSE  IT  n  I DENT  CONSTANT  (BASE) 

THEN  u" )-  BASE.TCdl  ••  EXPONENT  (*  Uaa  INTERVAL  •*  ITAYLOR  *) 

ELSE  IP  IINP  (BASE.TCtD)  <-  0.0 
THEN  BEGIN 
WRITE LN) 

WRITELN  (’ERROR)  ITAYLOR  BASE  <-  ZERO'); 

8VR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END  (*  THEN  *) 

ELSE  BEGIN 

IP  BASE.T  <>  EXPONENT. T 
THEN  BEGIN 
WRITELN; 

WRITE  (’ERROR)  ••  OP  ITAYLOR  VARIABLES  WITH'); 

WRITELN  ('  UNEQUAL  SCALE  r ACTORS' ); 

SVR  (0)  (*  RETURN  TO  OPERATING  SYSTEM  •) 

END;  (*  IT  •) 

U  >-  EXPONENT  •  ITLN  (BASE); 

0  i-  ITEXP  (U); 

END;  (•  ELSE  •) 

RES  i-  U 

END;  (*  ITAYLOR  **  ITAYLOR  •) 

(*  mo  OP  IT  POW.LIB  »»»»»>»»»»»»»»»»»»»>»»  *) 
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C.6.  Reel  tod  Intsrvsl  Tsvlor  Function* 


(•  FIT  ms. LIB 


MU  AMD  INTXKVAL  TAYLOR  FUMCTIONS 


«««««« 


Contents t 

TSQR(T)  FUNCTION  TSQR  (Ti  TAYLOR)  t  TAYLOR  I 
TSQRT(T)  FUNCTION  TSQRT  (Ti  TAYLOR)  t  TAYLOR) 

TKXF(T)  FUNCTION  TSXP  (Tt  TAYLOR)  i  TAYLOR) 

TLM(T)  FUNCTION  TLM  (Tl  TAYLOR)  t  TAYLOR) 

T_SIH_COe  PROCKDURE  T_SIN_COS  (Tl  TAYLOR)  VAR  8,  Cl  TAYLOR )l 
TSIM(T)  FUNCTION  TSIN  (Tl  TAYLOR)  I  TAYLOR) 

TC08(T)  FUMCTIOM  TC06  (Tl  TAYLOR)  t  TAYLOR) 

TRUNGX(T)  FUMCTIOM  TRUNGB  (Ti  TAYLOR)  t  TAYLOR) 

TARCTAM(T)  FUMCTIOM  TARCTAN  (Tl  TAYLOR)  t  TAYLOR) 

TTAM(T)  FUMCTIOM  TTAM  (Tt  TAYLOR)  t  TAYLOR) 

TDIFF(T)  FUNCTION  TOIFF  (Tt  TAYLOR)  I  TAYLOR) 

TIMTGRL(T)  FUMCTIOM  TINTS RL  (Tl  TAYLOR)  I  TAYLOR) 

ITSQR(IT)  FUNCTION  IT8QR  (Tl  ITAYLOR)  I  ITAYLOR) 

I TSQRT (IT)  FUMCTIOM  I TSQRT  (Ti  ITAYLOR)  I  ITAYLOR) 

ITEXP(IT)  FUMCTIOM  ITBXP  (Tl  ITAYLOR)  I  ITAYLOR) 

ITLM(IT)  FUMCTIOM  ITIM  (Tl  ITAYLOR)  t  ITAYLOR) 

IT_8IN_COS  PROCBDUR1  IT_8IM_Coe  (Tl  ITAYLOR)  VAR  8,  Cl  ITAYLOR)) 
ITSIN(Yt)  FUMCTIOM  ITs7w  'Tl  ITAYLOR)  I  ITAYLOR) 

ITCOS(IT)  FUMCTIOM  ITCOS  (Tl  ITAYLOR)  I  ITAYLOR) 

ITRUHGK(IT)  FUMCTIOM  ITRUMGI  (T)  ITAYLOR)  I  ITAYLOR) 

ITARCTANf IT)  FUMCTIOM  1TAMCTAM  (Ti  ITAYLOR)  i  ITAYLOR) 

ITTAM(IT)  FUMCTIOM  I  TTAM  (Tl  ITAYLOR)  I  ITAYLOR) 

ITDIFF(T)  FUMCTIOM  ITDIFF  (Tl  ITAYLOR)  I  ITAYLOR) 

ITIMTGRL(T)  FUMCTIOM  ITIMTGRL  (Tl  ITAYLOR)  I  ITAYLOR) 


FUMCTIOM  TSQR  (Tt  TAYLOR)  I  TAYLOR) 

(*  Requires ■  VRMULL,  SCALP,  SQR 

VAR  I,  J,  K,  HALF  I  DIMTYFI) 

X,  Yi  RVBCTOR) 

U  i  TAYLOR) 


(•  TMQR(T)  *) 


BIGIN 

X  i-  VRMULL)  Y  i-  VRMULL) 
0.LINST8  i-  T.LXNGTH) 

U.T  l“  T.T) 

U.TC  i-  VRMULL) 

U.TC(1)  i-  SQR  (T.TC t 1 1 ) ) 
I[1]  i-  T.TC(IJ) 
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I 


FOR  X  i-  2  TO  0. LENGTH  DO 

usni 

X(Xl  «»  T.TC(K]  i 
BALT  l-  X  DJV  2} 

FOR  J  i-  1  TO  HALF  DO 
BEGIM 

t  »•  x  -  a  ♦  i» 

YtJ)  i-  T.TCU1 
BNP;  (*  FOR  J  *) 

O.TCIX]  I-  2.0  •  SCALP  <  X,  Y,  0)» 

IF  X  MOO  2  -  1  THEM 
BEGIN 

HALF  «-  HALF  «■  1} 

O.TClXJi-  D.TCIX]  +  8 (JR  (T.TClHALF)) 
BHD  (*  IF  *) 

BHDi  (*  FOR  X  •) 

T8QR  >-  0 

END/  (*  FUNCTION  TSQR  (TAYLOR)  *) 


FUNCTION  TSQRT  (Ti  TAYLOR)  I  TAYLOR I 

(*  Requires i  VRNOLL,  SORT,  SQR,  SCALP  *) 
VAR  I,  X,  INDEX ,  INDX2:  DIKTYPEl 
RATIO}  REAL; 

X,  Y  t  RVXCTORl 
0  I  TAYLOR} 


(*  TSQRT(T)  *) 


BEGIN 

IF  (T.TCI1J  <  0.0)  OR 

((T.TC(1)  -  0.0)  AMD  (T.LmOni  >•  2)) 

THEN  BEGIN 

writsuii 

mtlTBLN  ('ERROR}  SQUARE  ROOT  or  TAYLOR  VARIABLE  <-  ZERO.')) 


SVR  (0)  END | 
(•  ELSE  •> 


<•  RETURN  TO  OFERATXNS  SYSTEM  *) 


0. LENGTH  1-  T. LENGTH/ 

U.T  «■  T.Tl 

X  i-  VRMOLLl  Y  I"  VRMOLLl 
U.TC  }■  VRMOLLl 


O.TC(I)  I-  SQRT  (T.TCtlJ)l 


(*  REAL  *) 


IF  0. LENGTH  >-  2  THEM  BEGIM 
RATIO  ««  1.0  /  (2.0  •  O.TCtl] )» 

O.TC(2]  T>TC(2]  *  RATIO} 

IP  0. LENGTH  >-  3  THEM  BEGIN 

O.TC131  i-  (T.TC13]  -  BOR  (O.TC121))  •  RATIO} 


.  aJ'*  ,-Vi'  : .  (Jf  { 


FOR  X  I-  4  TO  0. LI NOTH  DO 
BIGIN 

IP  X  WOO  2-0 
THB)  BEGIN 

INDEX  l-  X  DIV  2; 

INDX2  »-  INDEX  -  1> 

Xf INDX2]  «-  D.TCflNDEX] 

END  (*  THEN  *) 

ELSE  BEGIN 

INDEX  !-  (X  ♦  1)  DIV  2; 

O.TCtX]  SQR  (O.TC [INDEX] ) 

SKS;  (*  IF  *) 

FO*  I  »-  1  TO  INDX2  DO 
BEGIN 

INDEX  !-  X  -  II 
Y(X]  «■  O.TC [ INDEX] ; 

END;  (*  FOR  I  •) 

O.TC[X]  i-  (T.TC [X]  -  O.TC [X]  -  2.0  *  SCAUP  (X,  Y,  0))  •  RATIO 
END  (*  FOR  X  *) 

END  (•  IF  U. LENGTH  >-  3  «) 

END;  (•  IF  0. LENGTH  >-  2  ») 

TSQRT  t-  0 

END;  («  FUNCTION  TSQRT  (TAYLOR)  *) 


FUNCTION  TEXP  (Ti  TAYLOR)  t  TAYLOR; 
(*  toqulraai  VRHULL,  SCALP  *) 
VAR  J,  X,  INDEX.  DIMTYPE; 

RATIO.  REAL; 

X  i  RVECTOR; 

0  .  TAYLOR; 


X  I-  VRNULL; 

0. LENGTH  T. LENGTH; 

O.T  .«  T.T; 

O.TC  1-  VRNULL; 

O.TCtl]  !-  EXP  (T.TCM1  ); 


(*  TEXP(T)  *) 


FOR  X  »-  2  TO  0. LENGTH  DO 
BEGIN 

RATIO  1.0  /  (X  -  1); 

FOR  J  2  TO  X  DO 
BEGIN 

INDEX  X  -  J  +  1; 

X[ INDEX]  I-  T.TC [J]  •  (J  -  1)  *  RATIO 
END;  (*  FOR  J  •) 

O.TCtX]  SCALP  (O.TC,  X,  0) 

END;  (*  FOR  X  *) 

TEXP  »-  O 

END;  (*  FUNCTION  TEXP  (TAYLOR)  •) 


(•  runt)  •> 


i  f 

! 


FUNCTION  TLN  (Ti  TAYLOR)  I  TAYLOR) 

(*  Rtqulrtai  VRMULL,  LN,  SCALP  •) 

VAR  J,  X,  moot  I  DXMTOIl 
RATIO,  RATI 2 1  RIAL) 

X  t  KVXCTORi 

0  i  TAYLOR) 

BEGIN 

IP  T.TCtl]  <■  0.0 

TRIM  BKGIM 
WRITXLM) 

WRITBLM  ClRRORi  LM  OP  TAYLOR  VARIABLE  <-  XERO.')) 

SVR  (0)  BHD)  (•  RETURN  TO  OPERATING  SYSTXM  *) 

(*  ELSE  *) 

X  »■  VRMULL; 

U. LENGTH  I-  T. LENGTH) 

O.T  t-  T.T) 

O.TC  »■  VRMULL; 

O.TCtU  l«  IN  (T.TCtl] )f 
RATIO  i-  1.0  /  T.TCtl] ) 

IF  0. LENGTH  >«  2  THEM  BEGIN 
O.TC[2]  i-  T.TC(2]  *  RATIO) 

FOR  X  i-  3  TO  0. LENGTH  DO 
BEGIN 

RATI2  t-  1.0  /  (X  -  1)| 

FOR  J  i-  2  TO  X-1  DO 
BEGIN 

INDEX  I-  X  -  J  +  1/ 

X[ INDEX]  l-  T.TCtJ]  *  (INDEX  -  1)  *  RATI2 
END)  (•  FOR  J  *> 

U.TC(X]  «-  RATIO  *  (T.TCfXJ  -  SCALP  (O.TC,  X,  0)) 

END  (•  FOR  X  •) 

END)  («  IF  U. LENGTH  >-  2  •) 

TLM  i-  O 

END)  (•  FUNCTION  TLM  •) 


PROCEDURE  T_8IN_C08  (Tl  TAYLOR)  VAR  S,  Cl  TAYLOR))  (*  TJSINJC08  •) 

(•  Aaqulrast  VRMULL,  SIN,  COS,  SCALP  •)  ”  ” 

VAR  J,  X,  INDEX i  DIMTYPE) 

RATIO)  REAL) 

X  i  KVECTOR) 

BEGIN 

X  i-  VRMULL) 

8. LENGTH  i-  T. LENGTH) 

S.T  »“  T.T) 

8.TC  i*  VRMULL) 

S.TCtl]  i-  SIN  (T.TCtl] )| 

C. LENGTH  t-  T. LENGTH) 

C.T  »”  T.T) 

C.TC  »-  VRMULL) 

C.TCt 1)  1-  COS  (T.TCtl] )) 
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FOR  X  i-  2  TO  T.LXNOTH  00 

nsni 

RATIO  «•  1.0  /  (K  -  1)1 
FOR  J  i-  2  TO  R  DO 
BEGIN 

INDEX  l-  X  -  J  +  1 1 

X[ INDEX)  I-  T.TCUJ  •  (J  -  1)  *  RATIO | 
END)  (•  FOR  J  *) 

S.TCCE]  l-  SCALP  (C.TC,  X,  0); 

C.TC(X)  l-  -  SCALP  (S.TC,  X,  0) 

END)  <*  FOR  X  •) 

ENDt  (*  PROCEDURE  T  SIR  COR  •> 


FUNCTION  TSIN  (Tl  TAYLOR)  t  TAYLOR)  t*  TSZN(T)  *> 

(•  Requires i  VRNULL,  T_8IN_COB,  SIN.  COR.  SCALP  •> 

VAR  8,  Cl  TAYLOR)  ~  “* 

BEGIN 

TJBINJCOR  <T,  8,  C); 

TSIN  7-  8 

END)  <*  FUNCTION  TSIN  (TAYLOR)  *) 


FUNCTION  TCOS  (Tl  TAYLOR)  I  TAYLOR)  (*  TCOR(T)  •) 

(*  Requires i  VRNULL,  T_8IN_C06,  SIN,  COR,  SCALP  •) 

VAR  8,  Cl  TAYLOR;  ~ 

BEGIN 

T_SIH_COS  (T,  S,  C)> 

TC08  i-  C 

END)  (*  FUNCTION  TCOS  (TAYLOR)  •) 


FUNCTION  TRUNGR  (Tl  TAYLOR)  I  TAYLOR;  (*  TRDNQE(T)  •) 

(•  Requires!  X/T.  TSQR  *) 

VAR  Ui  TAYLOR; 

BNGIN 

U  i-  TSQR  (T); 

O.TCH)  i-  U.TCHJ  ♦  1.0; 

O  i-  1  /  O; 

TRUDGE  ;■  U 

END;  (•  FUNCTION  TRUNGR  •) 


FUNCTION  TARCTRM  (Tl  TAYLOR)  I  TAYLOR;  {•  TARCTAN(T)  *> 

(•  Requires i  VRNULL,  TRUNGR,  K/T,  TSQR,  ARCTAN,  SCALP  •) 

VAR  J,  X,  INDEX;  DIMTYPE; 

RATIOl  REAL; 

X  i  RVBCTOR; 

U,  V  i  TAYLOR; 


x  i-  twoui 

O.LXHOTH  »“  T.LXHGTHf 

0.T  »■  T.Tl 

O.TClU  «“  AHCTAH  <T.TC[1])» 


FOR  X  t-  2  TO  O.LXHOTH  DO 
BXGIN 

RATIO  t-  1.0  /  (X  -  Dl 
FOR  J  i-  2  TO  X  DO 
MSXM 

XMDXX  i-  X  -  J  ♦  1» 

X [ INDEX]  i-  T.TCIJ]  *  (J  -  1)  *  RATIO! 
XHDl  (•  FOR  J  •) 

U.TCtK]  t-  SCALP  (V.TC.  X,  0) 

XHDl  <*  FOR  X  *) 

TARCTAN  «-  0 

XHDl  (*  FOHCTIOR  TARCTAH  •) 


fuhction  ttah  <t«  tatlor)  »  taylor; 

(•  MqolrMi  VRHOLL,  T/T,  T_SIH_C06,  SCALP  *> 

VAR  8,  Cl  TAYLOR I 

BXOXH 

T_«IH_COS  (T,  8,  C)i 
TTAH  I"  S  /  C 

XHDl  (•  FOHCTIOR  TTAH  •) 


(*  TTAH(T)  *> 


FOHCTIOR  TDXFF  (Ti  TAYLOR)  I  TAYLORt 
VAR  X,  INDEX l  DXMTYPXi 
RATIO!  RXALi 
0  i  TAYLOR; 


(•  TDXPF(T)  •) 


BXOXH 

O.T  !■>  T.T» 

IF  T.LXHOm  -  1  THXH 
BXOXH 

O.LXHOTH  «-  1; 

O.TC[1)  »*  0.0 

XHD  (•  THXH  *) 

KL8X  BXOXH 

O.LXHOTH  l-  T.LXHOTH  -  1| 

RATIO  i-  1.0  /  T.T; 

FOR  X  l-  1  TO  O.LXHOTH  DO 
BXOXH 

INDEX  i-  X  ♦  If 

O.TC(X)  i-  T.TCtIHDEX)  •  RATIO  •  X 
XHD  («  FOR  X  •) 

XHDl  (•  SLOE  *) 

TDIFF  «■  O 

XHD;  (*  FOHCTIOR  TDXFF  •) 


f 

tp 

•V 


FUNCTION  TINTGRL  (Tt  TAYLOR)  >  TAYLOR) 

VAR  K,  INDEX)  D INTYPE) 

0  I  TAYLOR) 

BEGIN 

IP  T. LENGTH  <  DIM  THEN  U. LENGTH  I-  T. LENGTH 
ELSE  0. LENGTH  1»  DIN) 

U.T  t«  T.T» 

U.TCll]  t-  0.0) 

FOR  X  i-  2  TO  0. LENGTH  DO 

BEGIN 

INDEX  K  -  1) 

U.TC[K]  I-  T.TC [INDEX]  *  T.T  /  INDEX 
END)  (•  FOR  X  •) 

TINTGRL  1-  0 

END)  (•  FUNCTION  TINTGRL  •) 


FUNCTION  IT  SCR  (Tt  ITAYLOR)  i  ITAYLOR) 

(•  MqttlrMt  IVRNULL,  ISQR,  I  SCALP  *) 

VAR  I,  J,  X,  HALF:  OIHTYFB) 

X,  Yl  I  VECTOR) 

U  I  ITAYLOR) 

BEGIN 

U. LENGTH  I-  T. LENGTH) 

U.T  1“  T.T) 

X  I-  IVRNULL)  Y  )<*  IVRNULL; 

U.TC  I*  IVRNULL) 

U.TC[1]  ISQR  (T.TC[1])) 

X(1J  t-  T.TC(U) 

FOR  X  )■  2  TO  U. LENGTH  DO 
BEGIN 

X[X]  I-  T.TC [X]  ) 

HALF  i-  X  DIV  2) 

FOR  J  1-  1  TO  HALF  DO 
BEGIN 

I  !-  X  -  J  ♦  1) 

Y[J]  t-  T.TC [I] 

END)  (•  FOR  J  *) 

U.TC(X]  t-  2  *  I SCALP  (  X,  T,  HALF)) 

IF  X  HOD  2  -  1  THEN 
BEGIN 

HALF  t-  HALF  ♦  1) 

U.TClXJt-  U.TC [X]  ♦  ISQR  (T.TC [HALF]  ) 
END  (*  IF  *) 

END)  (*  FOR  X  *) 

ITSQR  »-  U 

END)  (•  FUNCTION  ITSQR  (ITAYLOR)  *) 


(•  TINTGRL(T)  *) 


(*  ITSQR (IT)  *) 


V. 


(•  XTSQRT(XT) 


POM CTIOM  ITSQRT  (Tl  ITATLOR )  <  ITATLOR  I 

(•  tafttlrui  IVRNULL)  ISQRT,  IT_IDEirr_COW8TAWT(  IT) , 

IBQR,  ISCALP  *)  ” 

VAR  X,  X,  INDEX,  XMDX2t  DXMTYPE; 

RATIO i  INTERVAL; 

X,  Y  f  I VECTOR I 
0  i  ITATLOR; 

BEGIN 

xr  (iinp  (t.tcH]  <  o.o)  or 

UIINP  (T.TCMl  -  0.0)  AND  (T. LENGTH  >-  2)) 

THIN  BEGIN 
WRITS LH; 

WRITE LN  ( ' ERROR 1  SQUARE  ROOT  or  ITATLOR  VARIABLE  <-  ZERO .')» 
SVR  (0)  Ml  (*  RETURN  TO  OPERATING  SYSTEM  *> 

(*  ELSE  *) 

X  «•  XVRNOLL)  Y  l-  IVRNULL; 

U. LENGTH  r-  T. LENGTH; 

U.T  l-  T.T; 

U.TC  I-  XVRNULL; 

U.TCtl]  tm  ISQRT  (T.TCtlJ)l  (*  INTERVAL  *) 

IP  NOT  IT  IDEWTJWHBTANT  (T)  THEN  BEGIN 
RATIO  «-  7  /  (2~*  U.TCtl]  )l 

IP  U.UWOTB  >»  2  THEN  BEGIN 
U.TCI2]  i-  T.TC (2)  *  RATIO; 

IP  U. LENGTH  >-  3  THEN  BEGIN 

U.TCI3)  i-  (T.TCI3J  -  IBQR  (0.TCI2J))  *  RATIO; 

POR  K  i-  4  TO  0. LENGTH  DO 

BEGIN 

IP  K  MOO  2*0 
THEM  BEGIN 

INDEX  t-  X  DXV  2; 

INOX 2  i*  INDEX  -  1; 

X [INDX2]  i*  U.TC [INDEX] 

END  {•  THEN  *) 

ELSE  BEGIN 

INDEX  i*  (X  ♦  1)  DIV  2; 

U.TC(X]  IBQR  (U.TC [INDEX] ) 

END;  (*  IP  *) 

POR  I  i*  1  TO  XNDX2  DO 
BEGIN 

INDEX  i-  X  -  I; 

Y(X]  l*  U.TC [INDEX] ; 

END;  (•  POR  X  •) 

U.TCtX]  l-  (U.TC [X]  -  2  •  I8CALP  (X,  Y,  INDX2) )  •  RATIO 
END  (*  POR  X  •) 

END  (*  IP  U. LENGTH  >-  3  •) 

END  (*  IP  U. LENGTH  >-  2  *) 

END;  (*  IP  NOT  XT_X DENTJCON8TANT  (T)  *) 

1TSQRT  i*  U 

END;  (•  FUNCTION  IT8QRT  (ITATLOR)  *) 


1 V*/.  Ml.  mm 


(•  ITEXP(IT)  •> 


FUNCTION  ITEXP  (Tt  ITAYLOR)  I  ITAYLOR; 
(*  Aaquiraai  IVRNULL,  I  SCALP  •) 
VAN  J,  K,  INDEX i  DIMTYPE; 

RATIO!  INTERVAL | 

X  l  IVECTORj 
0  i  ITAYLOR; 


I-  IVRNULLf 
«-  T. LENGTH  I 
»-  T.T / 
i-  IVFNULLl 


U.TCdJ 


IXXP  (T.TCHl  )/ 


FOR  X  i"  2  TO  D. LENGTH  DO 
BEGIN 

RATIO  !-  1  /  INTPT  (E  -  1)); 

FOR  J  l-  2  TO  X  DO 
BEGIN 

INDEX  i-  X  -  J  +  t; 

X[ INDEX]  «-  T.TC [J]  •  (J  -  1)  *  RATIO 
END;  (*  FOR  J  *) 

U.TCfK]  I-  I  SC  ALP  (U.TC,  X,  X) 

END;  (*  FOR  X  *) 

ITEXP  1-  U 

END;  (•  FUNCTION  ITEXP  (ITAYLOR)  •) 


FUNCTION  ITUI  (Tt  ITAYLOR)  I  ITAYLOR; 

(•  Rsquiraat  IVRNULL,  ILN,  ISCALP  •) 
VAR  J,  X,  INDEX  :  D INTYPE; 

RATIO,  RATI 2 >  INTERVAL; 

X  ;  I VECTOR; 

U  S  ITAYLOR; 


(•  ITUI  (IT)  *) 


BEGIN 

IF  I  INF  (T.TCtU)  <-  0.0 

THEN  BEGIN 
NRITELN; 

NRITELN  (’ERROR)  LN  OF  ITAYLOR  VARIABLE  <-  ZERO.’); 
SVR  (0)  END;  (*  RETURN  TO  OPERATING  8 

(•  ELSE  •) 


X  i-  IVRNULL; 

0. LENGTH  i-  T. LENGTH; 

U.T  i"  T.T; 

U.TC  i-  IVRNULL; 

U.TCtl]  ILN  (T.TC[1] ); 
RATIO  t-  1  /  T.TC[1] ; 


IF  U. LENGTH  >-  2  THEN  BEGIN 
U.TC(2]  i-  T.TC(2]  •  RATIO; 


■ 


1) 


TOR  X  j-  3  TO  D.LSHOTH  00 

non 

RATI 2  i-  1  /  XMTFT  (X  -  1>| 

FOR  J  1-  2  TO  K-1  DO 
MGIN 

non  i-  x  -  j  ♦  ii 

XtZMOn)  »-  T.TC[J)  •  RATI 2  •  (IMDBX  - 
m >1  (•  FOR  J  •) 

U.TCIXJ  I-  RATIO  •  (T.TC(K)  -  ISCALP  (O.TC,  X,  X)  > 
no  (*  FOR  X  *) 

SMDl  (•  IF  O.LXWTH  >-  2  •) 

ITU  i-  0 

KMDl  (•  FOMCTIOH  ZTLM  •) 


FROCXDORX  IT  an  COR  (T<  ItATLORl  VAR  S,  Cl  ITAYLOR) |  (• 

(•  RaquTraaT  ZVIMOLL,  I SIRS,  I COS,  ISCALP  *) 

VAR  J,  X,  non  ■  DZMTTFRl 
RATZOi  INTERVAL 1 
X  i  IVXCTORi 

awn 

X  >-  XVRMOLLf 

8.UH0TH  l-  T.LMQTHl 
8.T  I-  T.Tl 

8.TC  i-  IVRMULLl 

S.TCtn  i*  zszm  (T.IC(tJ); 

C.LKMGTH  tm  T.UBIOTRi 
C.T  1“  T.Tl 

C.TC  «-  IVRHOLLl 
C.TCI1]  I-  zcoa  (T.TC[1))» 

FOR  X  i-  2  TO  T. LENGTH  00 
SWIM 

RATIO  l-  1  /  XMTFT  (X  -  1); 

FOR  J  i-  2  TO  X  DO 

awn 

zmin  i-  x  -  a  *  ii 

XIZRDRX]  t-  T.TC(J]  *  RATIO  *  (J  -  1)i 
■HD I  (•  FOR  J  •) 

S.TC(K)  t-  I SCALP  (C.TC,  X,  K)| 

C.TCIX]  i-  -  ISCALP  (S.TC,  X,  X) 

■MO I  (•  FOR  X  •) 

mo i  (•  procsdorr  it  an  DOS  •) 


FOMCTIOH  nan  (Ti  itatlor)  i  itaylor  (  -  (• 

(*  Raquiraai  IVRROLL,  IT_aiM_COa,  ISIMX,  I COS,  ISCALP  •) 
VAR  a.  Cl  ITATLOR r 

awn 

ir_an_co«  (t,  a,  c)? 
mn  »•  a 

BMDl  (•  FOMCTIOH  ITS  III  (ITATLOR)  •) 


it  an  cos  •) 


ZTSn(ZT)  *) 


FUNCTION  CTC06  (Tt  ITAYLOR)  I  ITAYLORl  (*  VKOeilt)  *) 

(*  Requires i  IVRNULL,  IT_8IN_C08,  ISINE,  ICOS,  I  SCALP  •) 

VAR  8,  Cl  ITAYLOR;  ”  ” 

BIG  IK 

ITJSIN_C06  (T,  8,  C)l 
ITCOS  7-  C 

END;  I*  FUNCTION  ITCOS  ( ITAYLOR)  *) 


FUNCTION  ITRUNGE  (Tl  ITAYLOR)  i  ITAYLOR;  (*  ITRDNCK(IT)  *) 

(•  Requires!  K/IT,  ITSQR  •) 

VAR  Ot  ITAYLOR) 

BEGIN 

U  i-  ITSQR  (T) l 
U.TCll]  i-  U.TCll]  *  1) 

0  |«  1  /  Ol 
TRUNSE  1'  0 

END;  (*  FUNCTION  ITRUNGE  *) 


FUNCTION  ITARCTAN  (Tl  ITAYLOR)  1  ITAYLOR;  <*  ITARCTAN(IT)  *) 

(*  Requires i  IVRNULL,  ITRUNGE,  K/IT,  ITSQR,  IARCTAN,  I SCALP  *) 

VAR  J,  K,  INDEX  I  DIKTYPB; 

RATIOl  INTERVAL; 

X  I  IVECTOR; 

U,  V  i  ITAYLOR; 

BEGIN 

X  I-  IVRNULL; 

D. LENGTH  I-  T. LENGTH; 

U.T  1-  T.T; 

U.TCll]  1-  IARCTAN  (T.TCIU); 

(*  Generates  error  if  AB8  (T.TCIU)  >1.0  *) 

V  i-  ITRUNGE  (T); 

FOR  K  I-  2  TO  U. LENGTH  DO 
BEGIN 

RATIO  I-  1  /  INTFT  (K  -  1); 

FOR  J  i-  2  TO  X  DO 
BEGIN 

INDEX  i-  K  -  J  ♦  1; 

X( INDEX)  i-  T.TCWJ  •  (J  -  1)  •  RATIO; 

END;  (•  FOR  J  *) 

U.TCtK]  I-  I SCALP  (V.TC,  X,  R) 

END;  (•  FOR  K  •) 

ITARCTAN  l-  0 

END;  (*  FUNCTION  ITARCTAN  •) 


FUNCTION  ITT AN  (Tl  ITAYLOR)  i  ITAYLOR; 

(*  Requires  I  IVRNULL,  IT/IT,  IT_BIR_COS,  I SCALP  •) 
VAR  8,  Ct  ITAYLOR;  ”  ~ 

BEGIN 

IT_8IN_C08  (T,  8,  C); 

I1TAN  7»  8  /  C 

END;  (•  FUNCTION  ITT AN  •) 


(*  ITTAN(IT)  •) 
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<*  itdiff(it)  *) 


r 


1  function  itdiff  (ti  itaylor)  <  itaylor) 

VAR  X,  INDEX l  DIMTYPE) 

RATIO t  INTERVAL) 

0  i  ITAYLOR) 

BEGIN 

O.T  «-  T.T) 

"  IF  T. LENGTH  -  1  THEN 

BEGIN 

U.  LENGTH  l-  1| 

U.TC[1]  «-  0.0 

END  (*  THEN  •) 

j  ELSE  BEGIN 

O.LERGTH  >•  T. LENGTH  -  1) 

RATIO  i-  1  /  INTPT  (T.T) I 

FOR  (■•ItO  0. LENGTH  DO 
BEGIN 

INDEX  1-  X  +  1) 

U.TC(K]  i-  T < TC [ INDEX]  *  RATIO  *  X 
i  END  (*  FOR  X  *) 

j  END)  <*  ELSE  •) 

!  ITDIFF  i-  D 

l  END;  (•  FUNCTION  ITDIFF  •) 


FUNCTION  ITINTGRL  (Ti  ITAYLOR)  >  ITAYLOR)  (•  ITINTGRL(IT)  *) 

VAR  X,  INDEX)  D INTYPE J 
RATIO)  INTERVAL) 

U  >  ITAYLOR) 

BEGIN 

IF  T. LENGTH  <  DIN  ISBN  U. LENGTH  I-  T. LENGTH  ♦  1 
ELSE  U. LENGTH  I-  DIM) 

U.T  1-  T.T) 

U.TCtl]  »-  INTPT  (0.0)) 

RATIO  !-  INTPT  (T.T)) 


U 


FOR  X  l-  2  TO  U. LENGTH  DO 
BEGIN 

INDEX  i"  X  -  1) 

U.TC(K)  i-  T.TC [INDEX]  *  RATIO  /  INDEX 
END)  (•  FOR  X  •) 

ITINTGRL  i-  U 

END)  (•  FUNCTION  ITINTGRL  •) 


(•  END  OF  RIT  FN8.LIB  >»»»»»>>»»»»»»»»»»»»»»>  *) 


1 
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C.7.  Utility  ronctlona  and  **2SS£lE3t • 


I*  UTIL. LIB  -  UTILITY  PXOdDOm  I  IUHCTIOHB  «««««««« 


Contantai 

(unction  VRRULL  s  RVBCTORt 

roncrzon  ivrrull  «  ivbctor# 
ranctzon  t  zonnr  mo  (t«  nnw)  »  boolbar# 
roncrzon  t~i dbrt jxhbtart  (Ti  thrum)  i  boolbar; 
roncrzon  njnwnTsBM  (n  nutn)  >  nnuMf 
roncrzon  itjtDtnr“ corstart  (n  znnoi)  >  boolbar# 
procboorb  mtijRimu  (iwti  ammii 
pmocbppre  m*irt~iwnxvhLjmiMB  (mi  znnoi)i 
procboorb  read  Trtbkval  *bribb  (van  n  itaylor)  i 
pbocbpprb  mtznjnouM  Tr<  taylor)# 
roncrzon  naorr~ (r  i  itaylor)  «  taylori 


*) 


roncrion  vrroll  ■  kvrctor; 

tar  Zt  DDRYPSI  Ot  RVBCrORl 
BBOIR 

rOR  I  f  1  TO  DIM  DO  all I  »»  0.0# 
VRRULL  i-  O 

bhd#  (*  roncrzon  vrroll  *) 


roncrzon  ivrrull  <  i vector; 

VAR  Zi  DIMTYPR#  Ol  ZVSCTORl 
BBOIR 

cor  i  i-  i  to  dim  do  o(i]  i-  nrrPT  <o.o>» 

I VRRULL  I«  U 

bmd;  (•  roncrzon  ivrrull  *) 


roncrion  t_idbrt_mro  (Ti  taylor)  «  boolbah  » 

(•  trob  If  tin  aarlaa  la  ldantlcallr  aqnal  to  asro  •) 
VAR  0<  BOOLBAB# 

Kt  IBTBOBR# 

BBOIR 
0  1“  TROB# 

X  f»  I# 

RRZLB  (X  <»  T.LBROTH)  ARD  (0)  DO 

BBOIR 

ir  T.TCtX]  <>  0.0  TOR  0  I-  rALSBl 

X  i»  X  ♦  1 

BHD#  (•  BULB  •) 

r  iDtnr  mo  •>  o 

rib#  (•  roncrzon  r  ioobt  bbro  *) 


in 


FUNCTION  T  I  DOT  CONSTANT  (Tt  TATUM )  i  BOOUMl 

(*”  tiros- if  tho  Nriw  Is  identically  equal  to 
VAR  Ot  BOOUMl 
Ki  MTBGSRl 

■sen 

o  i*  tnobi 

X  t"  2| 

mau  (X  <-  T.LSNOTH)  AMO  (0)  00 

non 

IF  T.TC(K)  <>  0.0  THBM  0  I-  F ALAS I 

X  »•  X  ♦  1 

XMDl  {•  mu  *) 

T  IDXMT  CONSTANT  t-  D 

WXDl  (•  FUNCTION  TJCDSNT  CONSTANT  *> 


FUNCTION  CTwIDtNT_HB»0  (Tl  ITAVUMt)  I  MOOT  NAM  I 

<*  'rum  Tf  tho  MrlN  is  Identically  oqoal  to  ooro  *) 

VAR  Ul  BOOLXMl 

xi  intbgeri 

ssen 

0  <*  TRUE; 

X  «•  1» 

mu  (X  <•  T.LSNOTH )  AMD  (0)  00 

non 

IF  T.VC(K]  <>  INTFT  (0.0)  TMW  D  >•  FALSBi 

X  •-  X  ♦  I 

XMDi  (•  miu  *) 

IT  I DENT  ZERO  I-  0 

ESDI  («”  FUNCTION  IT_ID«Wr_*El»  (ITATL0R)  *) 


FUNCTION  IT_'  JtWT_CONETAET  (Ti  ITATLOR)  I  BOOUMl 
(•  “true  Tf  tho  oorioo  is  identically  equal 
to  «  constant  lntorral  *) 

VAR  0>  BOOUMl 
Xt  MTEQERi 

BEGIN 
U  t"  TRUE I 
X  i*  2| 

WHILE  (X  <-  T. LENGTH)  AND  (U)  DO 

SBGn 

IF  T.TCIX)  <>  INTFT  (0.0)  TEEN  0  I-  FALSBl 

X  *•  X  ♦  1 

END |  (*  WHILE  •) 

ITJtOBWTjOOHSTAMT  «-  0 

END|  (•”  FUNCTION  XT_I DENT  ^CONSTANT  (ITAVUMt)  •) 


FNOCEOORE  WRITE  INTERVAL  (IWTi  INTERVAL) | 
BEGIN  ~ 

writs  (•[•»  nr. inf.  *,  nr. sob.  •)•)! 

ENDi  (•  FNOCSDORS  WRITS  INTERVAL  •> 


proccdorc  vritc_irtcrval_ccaim  (inti  ntnoi) i 
vu  it  ozNnvsi  ^ 
amir 

ro*  x  t-  i  to  srr.lbroth  do  rmir 
man  (Xt4,  •  •)» 

WRXn_XHTCRVAI.  (SRR.TC[X])| 

writer  no 

bnd»  (•  wocwwi  man  irtrrvai  stun  *) 


proccdobc  ms  nmnu  nun  (vu  pi  itaylor); 

VAR  It  DURTRl 

RR8XM 

RBU)  (DATA,  r.T)» 

nmu  ('KALB  FACTOR i 1 ,  r .T)» 

FOR  I  «-  1  TO  P.LBRQTR  DO  IRCAD  (DATA,  P.TClI]  )| 

nunu  ( ' scrim  for  p«')» 

RlOn  IRTCKVALJMRICS  (P)  t 

RHDt  ”(•  PRodmORl  RRAD  IRTCRVAL  SERIES  •) 


PROCCDORC  HRin  IKRXBS  (T:  TATIiOR)  I 
VAR  It  DDRTPRj 
BRCIH 

POR  I  l-  1  TO  T.LBROTR  DO 

MRXTCUI  (I»4,  1  T.TCtll) 

RMDl  <•  PROCRDORR  MUT1  SCRIM  •> 


PtMCTIOM  TMIDPT  (P  i  ITAYLOR)  t  TAYLOR; 

VAR  It  DXMTYPVi  Ot  TAYLOR) 

BSQIH 

O.T  1“  P.Tl 
D.LCRGTH  I-  P.LCHOTH; 

POR  I  1  TO  O.LCRQTH  00 

0.TC(I]  t*  (ZIHP  <P.TC(I])  ♦  XSOP  (P.TC(I1  ) )  *  0.5t 
ITHZDPT  t-  O 

CRD;  (•  P0HCTZ0H  IWIDPT  *) 


(*  crd  or  orxb.hn 


»»»»»»»»»»»»»»»»»»»»»  •) 
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APMBDIX  0 


DMT  NMMl  for  KSCQJtOLV  and  IMQ_SOLV 


fhii  Appendix  in  addressed  to  tMitri  with  neennn  to  *ascal-K  «ho  wish  to 
an  the  programs  W*g_SOX<V  and  imq_solv  to  explore  the  solutions  of  initial- 
value  profelsas  for  ordinary  differential  equations* 

These  pro press  allow  direct  user  intervention  in  the  selection  of  the 
integration  step  sine  at  each  step.  As  such,  they  are  quite  useful  for  hands-on 
exploration,  hut  they  are  not  intended  to  serve  as  general  purpose  solvers* 

the  illustrations  given  here  are  for  the  Pascal-SC  conpller  for  the  Xilog 
xlo  operating  systsn  and  its  text  editor*  Users  with  the  Xilog  no-1 
aicroooagtuter  can  obtain  the  software  described  in  this  report  by  sending  a 
hard-sectored,  single-sided  I*  floppy  disk  to  the  second  author.  The  directions 
given  here  will  apply  with  only  Minor  uodificationa  to  the  use  of  fnscsl-K 
oonpllera  for  other  ay stuns. 

Step  1.  Create  the  source  progran. 


■OXT  Hyprob.f 

on  AMtQJMX.V.1  or  IW*g_SOLV.S 

CWWI  nun  «»  to  the  expression  for  the  equation  to  be  solved 
■space  the  source  oode  for  the  tag*  (ITSQft)  operator  by  the  operators 
and  functions  needed  to  evaluate  the  expression  in  your  equation i  they 
are  obtained  free  -.LIS  as  needed.  Many  operators  call  other 
operators,  which  suet  also  be  included  in  your  source  oode.  The  oalls 


to  tbo  buio  Uttml  operators  oro  already  UelaM  la  the  heading 
of  XDBQJMLV. 


the  proeeea  of  aaeanhllng  all  the  neoeeeary  operators  la  the  oorreot  order  la 
rather  tedious  aad  error-pronei  be  patient  and  careful,  a  good  library  eaaager 
would  be  an  aaeet. 

Step  2.  Oouplle. 

For  H»Q_aoC.V,  use i  CMP  Nyprob  HNLXB 

For  1MQJKS.V,  ueet  CMP  Nyprob  X_LX» 

Step  3.  If  the  compiler  deteota  errors,  ter  example,  an  omitted 
subroutine,  then  edit  Nyprob. 8  again  and  reooapl le. 

Step  4.  Zf  you  wish  to  reed  Input  data  frou  a  file  Instead  of  the  console, 
then  thla  file  oust  be  created,  the  Initial  Taluee  of  X  and  T  are  reed  fren  an 
external  file  which  can  be  naned  as  the  near  chooses.  Nyprog.DAT  will  be  need 
here  ter  Illustration.  Output  can  also  be  directed  to  an  external  file  instead 
of  the  console,  for  exaaple,  to  Myprog.OOT. 

Step  5.  Run. 

For  Input/out pot  using  the  terminal,  use 


XQP  Myprog  SOON  SCON  SOON 
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ABSTRACT  fCAMkw  m  199mm  BP  If  mmiM  Nf  IAMB  Sr  Mat  amaSot) 

Taylor  series  have  a  long  history  of  usefulness  in  numerical  analysis, 
especially  for  the  numerical  solution  of  the  initial  value  problem  for  systems 
of  ordinary  diffsrential  equations.  Since  recurrence  relations  for 
coefficisnts  of  Taylor  sarles  are  veil  known,  it  is  possible  to  automate  the 
computation  of  arithmetic  operations  and  various  standard  functions  with 
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20.  ABSTRACT  -  cont’d. 

argvjMnts  which  are  themselves  series.  If  the  language  used  for  scientific 
confutation  supports  user  defined  operators  and  data  types,  then  the 
facilities  built  into  the  language  compiler  itself  can  be  used  to  generate 
machine  code  for  the  evaluation  of  Taylor  coefficients.  Examples  of  such 
languages  are  Pascal-SC,  Algol  68,  and  ADA  (a  trademark  of  the  U.  S. 

Department  of  Defense) .  Pascal-SC  (Pascal  for  Scientific  Computation)  offers 
the  user  highly  accurate  floating-point  and  interval  arithmetic,  the  latter 
being  useful  for  automatic  confutation  of  guaranteed  error  bounds.  In  this 
language,  series  with  real  coefficients  are  introduced  as  type  TAYLOR,  and  the 
corresponding  series  with  interval  coefficients  as  type  ITAYLOR.  Source  code 
is  given  for  the  operators  +,  -,  *,  /,  **  and  the  functions  SQR,  SQRT,  EXP, 
SIN,  COS,  ARCTAN,  and  LN  with  arguments  of  these  types  and  some  other  useful 
functions  and  procedures.  Integer,  real,  and  interval  constants  are  also 
allowed  in  TAYLOR  or  ITAYLOR  expressions.  Suggestions  for  the  implementation 
of  additional  operators  or  functions  are  given.  An  application  of  Taylor 
series  and  the  methods  of  interval  analysis  to  the  solution  of  the  initial 
value  problem  for  ordinary  differential  equations  is  made  using  types  TAYLOR 
and  ITAYLOR.  An  analysis  of  the  stability  of  this  method  is  made,  which  shows 
that  the  recurrence  relations  for  generation  of  the  Taylor  series  for  the 
solution  exhibit  a  mild  instability  which  has  no  significant  effect  on  the 
values  of  the  solution  computed  by  analytic  continuation. 


